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B-1B  Avionics/ Automatic  Test  Equipment 
Maintenance  Queueing  Analysis 


I.  Background  and  Problem  Statement 


Background 

Modern  USAF  aircraft  have  complex  avionics  components 
which  require  -frequent  unscheduled  maintenance.  As  part  of 
their  maintenance  concept,  USAF  aircraft  are  delivered  to 
their  bases  with  complex  and  expensive  automatic  test 
equipment  (ATE).  Avionics  components,  called  line 
replaceable  units  (LRUs),  are  removed  from  the  aircraft  when 
a  malfunction  is  detected  and  are  taken  to  a  nearby  repair 
shop  where  the  ATE  is  located.  Maintenance  technicians  then 
use  the  ATE  for  fault  detection,  isolation,  and  ultimately 
repair  of  the  avionics  LRU.  The  ATE  consists  of  specialized 
test  stations.  Each  test  station  is  devoted  to  a  different 
grouping  of  the  avionics  LRUs.  For  example,  one  test 
station  might  be  used  for  the  repair  of  radio  frequency  <RF> 
LRUs,  while  another  test  station  eight  he  used  for  the 
repair  of  digital  computer  LRUs.  Since  there  are  many 
Aircraft  at  each  base,  and  since  each  aircraft  has  many 
complex  LRUs,  there  may  be  numerous  LRUs  sent  to  the  repair 
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thU  r«t»irch  effort  is  a  queueing  analysis  of  B-1B 
avionics  maintenance.  In  this  analysis,  the  queueing 
"customer's*  ars  ths  avionics  components  of  th»  B-1B,  and  ths 
queueing  "servers"  are  automatic  test  equipment  stations. 

To  solve  this  queueing  si tuition,  both  analytical  and 
simulation  techniques  were  considered.  In  looking  at  any 
queueing  situation,  1  believe  it  is  very  worthwhile  to  try 
both  simulation  and  analytical  techniques,  and  to  keep  an 
open  mind  about  which  technique  is  best  until  both  are 
developed. 

I  am  in  obligation  to  my  thesis  advisor,  Lt  Col  Jim 
Bexfield,  for  his  guidance  and  expertise  in  queueing  theory. 
My  thesis  reader,  Maj  Ken  Feldman,  provided  helpful  advice 
from  a  real-world  management  point  of  view.  I  would  also 
like  to  thank  Lt  Col  Dick  Diehl  and  Lt  Col  Tom  Clark  for 
getting  me  started  on  this  project.  Finally,  1  would  like 
to  thank  the  many  people  at  the  B-1B  System  Program  Office 
and  at  Headquarters  Strategic  Air  Command  who  helped  me  keep 
the  project  going. 


Lance  M.  Roark 
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The  purpose  of  this  research  effort  was  to  develop  a 
technique  to  determine  B-1B  automatic  test  equipment  (ATE) 
station  quantities  required  to  support  the  B-1B  avionics 
Component*  at  base  level.  As  part  of  this  effort,  both 
simulation  and  analytical  solutions  were  developed.  A 
detailed  and  complex  simulation  model  was  developed  in  the 
Q-GERT  simulation  language.  In  addition*  an  analytical 
model  was  developed  based  on  the  theory  of  open  queueing 
networks  and  other  queueing  technique*.  However*  the 
analytical  model  required  many  crude  and  simplifying 
assumptions*  and  the  analytical  results  were  not  entirely 
satisfactory.  The  Q-GERT  simulation  model  was  selected  as 
the  best  choice  for  the  remainder  of  the  research  effort. 

Two  techniques  were  developed  to  determine  test  station 
quantities  based  on  the  model  output.  The  first  technique' 
was  to  buy  sufficient  test  stations  to  achieve  a  four  day 
maximum  base  repair  cycle  time  for  the  avionics  components. 
The  second  technique  was  to  conduct  a  cost -benefit  analysis 
by  comparing  the  costs  of  additional  test  stations  (and  the 
benefits  of  shorter  repair  cycle  times)  to  the  benefits  of 
fewer  test  stations  (and  the  costs  of  longer  repair  cycle 


times).  Considerable  sensitivity  analysis  was  pe formed 
the  simulation  model,  and  the  research  effort  concludes 
a  range  of  management  options  for  consideration  by  the  E 
System  Program  Office. 


shop  each  day.  This  may  generate  a  significant  workload  on 
each  of  the  test  stations. 

The  workload  of  avionics  on  ATE  test  stations  is  of 
particular  interest  to  the  B-IB  System  Program  Office  (SPO) . 
The  ATE  which  will  be  deployed  with  the  B-IB  bomber  consists 
of  lour  specialized  test  station  types  which  will  be  used  to 
support  over  100  avionics  LRUs.  Unless  the  B-iB  SPO  has  the 
ability  to  quantify  and  predict  this  workload,  the  SPO  will 
be  uncertain  as  to  the  best  quantity  of  test  stations  of 
each  type  to  procure  for  each  B-IB  base.  Not  buying 
sufficient  test  stations  would  degrade  B—1B  avionics 
readinessl  buying  too  many  test  stations  would  be  needlessly 
expensive.  A  technique  to  estimate  the  avionics  workload  on 
the  ATE  would  not  only  be  useful  for  decision-making  about 
the  best  quantity  of  test  stations  to  procure,  but  could 
also  be  u*ed  to  justify  ATE  funding  requirements.  Due  to 
the  tremendous  lead  times  involved  in  the  procurement 
process,  such  decision-making  must  tak#  place  very  early  in 
the  life  of  the  B-IB  program. 

Stategjeol  of  etSfelBB 

The  B-IB  SPO  needs  a  technique,  and  ultimately  a  model, 
which  can  be  used  to  analyze  ATE  test  station  requirements. 
The  B-IB  SPO  will  continually  update  its  information,  and 
therefore  needs  its  own  model  and  user’s  manual .  This  model 


would  bo  used  to  examine  the  tradeoffs  of  cost  versus 
avionics  readiness.  This  model  would  also  be  useful  for 
"what  if"  trade  studies  (for  example,  impact  of  different 
operational  or  maintenance  concepts). 

Qbigctiye  q±  tbt  Rg*e«c;ch 

The  overall  objective  of  the  research  effort  was  to 
provide  the  B-IB  SPO  with  a  computerized  model  which 
provides  the  capability  to  assess  the  avionics  maintenance., 
workload  on  ATE  test  stations.  Further,  the  research  effort 
needed  to  develop  criteria  and  procedures  for  selecting  the 
best  quantities  of  stations  once  the  workload  had  been 
measured.  This  was  accomplished  by  the  sequential 
attainment  of  the  following  subobjectives: 

1.  It  was  essential  to  develop  a  detailed  and  accurate 
description  (model)  of  B-iB  avionics  maintenance.  As  part 
of  this  description,  it  was  necessary  to  consider  all 
possible  factors  that  could  have  a  bearing  on  the  avionics 
maintenance  workload  on  ATE  test  stations.  This  description 
was  obtained  from  a  review  of  various  B-IB  logistics  and 
maintenance  planning  documents  and  also  from  personal 
interviews  with  personnel  from  the  B-IB  SPO  and  HQ  SAC. 

This  conceptual  description  of  avionics  maintenance  was  to 
become  the  framework  of  all  subsequent  model  development. 


2.  Sine*  the  rtMvch  effort  involved  *  practical 
application,  it  was  necessary  to  obtain,  collect,  and  review 
relevant  data  needed  for  the  model  development.  B-1B 
operational  data  <such  as  number  of  aircraft:  per  bape* 
flying  hour*  per  aircraft  per  month,  etc.)  were  obtained 
from  the  B-1B  SPO.  Reliability  and  maintainability 
estimates  for  each  avionics  LRU  were  obtained  from  the  B-1B 
associate  contractors.  In  addition,  operational  reliability 
and  maintainability  data  were  colloeted  from  actual 
experience  with  the  F-16  fighter  and  the  B-S2  Offensive 
Avionics  System  <OAS)  update.  The  data  were  examined  by 
goodness-of-f it  tests  and  other  techniques  to  determine  the 
most  realistic  probability  distributions  *hieh  were 
eventually  incorporated  during  the  model  development. 

3.  A  very  detailed  and  complex  sibulatlon  model  of  B-1B 
avionics  maintenance  was  developed.  This  model  was 
developed  in  the  O-SCRT  simulation  language,  the  model 
simulates  the  fide  Of  avionics  LRUs  from  LRU  failure  to  LRU 
repair  (and  return  to  base  supply).  Once  the  model  was 
developed,  the  next  step  was  to  design  the  simulation 
experiment.  As  part  of  this  experimental  design,  variance 
reduction  techniques  were  used.  The  model  was  also  used-  for 
sensitivity  analysis  since  much  of  the  data  inputs  were 
preliminary  contractor  estimates. 


4. 


A  second  model,  based  on  queueing  theory,  was 


developed.  The  analytical  modal  Mas  based  on  the  theory  Of 
queueing  networks  and  other  advanced  techniques.  The 
analytical  model*  not  surprisingly,  requires  many  more 
simplifying  assumptions  then  does  the  simulation  model.  The 
queueing  theory  approach  includes  writing  the  balance 
equations  and  then  solving  these  equations  by  a  numerical 
technique.  The  research  effort  also  examined  the  tradeoffs 
between  simulation  and  analytical  techniques  for  this 
particular  application  and  concluded  that  the  simulation 
model  was  best  suited  for  the  B-1B  SPO. 

5.  The  research  effort  also  included  two  tradeoff  studies 
on  avionics  maintenance.  First,  it  was  necessary  to  conduct 
a  cost  "-benefit  analysis  on  procurement  of  additional  test 
At at ions  versus  procurement  of  additional  LRU  spares.  This 
was  part  of  the  effort  in  determining  the  best  quantities  of 
stations  to  procure.  Second ,  the  impact  of  a  hypothetical 
6— IB  deployment  for  a  major  conventional  war  was  considered. 


Liitcsfcuc*  S«K£b  sod  B*yi*w 

Since  this  proposed  thesis  research  involved  a  queueing 
analysis  of  B-18  avionics,  the  focus  of  the  literature 
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search  has  been  on  the  mathematical  techniques  of  queueing 
theory.  In  queueing  theory,  "customers"  stand  in  line 
(called  a  queue)  awaiting  service  by  some  type  of  service 
facility.  The  time  between  arrivals  of  the  customer  and 
also  the  time  to  serve  a  customer  are  usually  described  by 
some  probability  distribution.  In  this  thesis  effort,  the 
customers  represent  avionics  LRUs,  and  the  act  of  service 
represents  maintenance  of  a  LRU  on  the  ATE  test  station. 
Queueing  theory  will  be  used  to  determine  the  minimum  number 
of  test  stations  (by  type)  which  can  accommodate  the 
maintenance  workload  at  each  base. 

With  any  queueing  system,  the  system  must  be  designed 
so  that  the  mean  system  capacity  to  process  customers  can 
handle  the  average  workload  of  customers.  If  this  is  not 
the  case,  a  "traffic  jam"  or  bottleneck  will  result  and  the 
queue  site  will  grow  indefinitely.  But  even  if  the  mean 
system  capacity  can  handle  the  average  workload,  a  "traffic 
jam"  can  still  occur  because  the  actual  volume  of  workload 
fluctuates  (in  a  statistical  sense).  Because  of  these 
fluctuations,  the  queue  size  may  become  too  big  or  the  time 
a  customer  spends  in  the  queueing  system  may  be  unacceptably 
high  (1:1).  The  main  purpose  o#  queueing  analysis  is  to 
precisely  measure  the  queue  size  and  the  time  a  customer 
spends  in  the  system. 

Queueing  theory,  in  its  most  ambitious  form,  can  be 
used  to  analyze  a  queueing  network  or  syatem.  A  queueing 
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network  is  often  represented  as  a  set  of  queues,  and  the 
network  may  process  many  different  types  or  classes  of 
customers.  The  individual  queues  are  usually  described  in 
terms  of  numbers  of  servers,  queueing  (scheduling) 
disciplines,  and  probability  distributions  for  the  service 
and  arrival  times  at  each  of  the  queues.  Some  queueing 
networks  can  be  represented  as  Markov  processes.  A  Markov 
process  represents  a  system  as  a  generally  discrete  set  of 
disjoint  states.  In  the  case  of  queueing  theory,  the  state 
of  the  system  or  network  is  usually  the  number  of  customers 
awaiting  service  at  each  of  the  queues.  The  set  of  states 
and  the  transitions  between  the  states  define  a  set  of 
linear  equations.  By  adding  an  equation  mhich  states  that, 
the  state  probabilities  must  Sdd  to  one,  it  is  possible,  at 
least  in  principle,  to  derive  the  equilibrium  probability  of 
each  state  (2:25).  However,  more  complicated  situations 
cannot  be  solved  in  practice,  and  it  is  necessary  to  resort 
to  computer  simulation  to  obtain  any  results  for  the 
queueing  analysis. 

If  the  state  probabilities  can  be  obtained,  it  is 
usually  desirable  to  condense  this  information  into 
performance  measures  of  the  queueing  system.  These 
performance  measures  can  then  be  used  to  assess  the  adequacy 
of  the  queueing  system.  In  the  case  of  simulation  models, 
the  performance  measures  are  typically  estimated  directly. 

In  any  event,  the  most  common  measures  of  performance  for  a 


queueing  system  are  the  average  number  of  customers  in  the 
system,  the  average  number  of  customers  waiting  in  queue, 
the  average  amount  of  time  a  customer  spends  in  the  system, 
and  finally  the  average  amount  of  time  a  customer  waits  in 
queue  (3:261).  In  some  cases,  it  is  passible  to  compute 
these  performance  measurements  directly  without  explicitly 
computing  the  system  state  probabilities. 

The  most  common  queueing  model  is  a  single  server 
exponential  model.  "Single  server"  means  that  there  is  only 
one  server  to  process  customers I  customers  stand  in  line 
before  this  server  and  wait  for  their  turn  to  be  served. 
"Exponential"  means  that  both  the  time  between  customer 
arrivals  and  the  time  to  service  each  customer  have 
exponential  probability  distributions.  This  is  a  relatively 
simple  and  papular  model,  and  it  can  be  found  in  many  texts 
such  as  Ross.  Ross  has  described  the  mathematics  of  the 
state  probabilities  and  the  system  performance  measurements 
in  full  detail  (3:264-267). 

There  are  many  extensions  to  the  single  server 
exponential  model.  Ross  has  developed  the  solution  to  the 
case  where  the  queue  for  the  single  server  has  a  finite 
capacity  (3:270-272).  Cooper  presents  the  case  where  the 
arrival  of  customers  follows  a  "quasi-random  process".  In 
such  a  process,  the  requests  for  service  are  generated  by  a 
finite  number  of  sources  where  each  source  follows  an 
exponential  distribution  (4:102-116).  This  contrasts  with 


the  simple  single  server  exponential  model  described 
earlier,  because  in  the  simple  model  the  customers  arrive 
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•from  an  infinitely  large  population  of  sources.  Finally, 
Kleinrock  describes  the  most  general  exponential  model.  It 
allows  for  the  queueing  system  to  have  more  than  one  server, 
a  finite  capacity  queue,  and  a  finite  population  of  sources 
(5:108-110).  Almost  all  other  exponential  models  are 
special  or  limiting  cases  of  this  general  model. 

Cooper  also  develops  a  model  which  allows  for  different 
types  of  customers.  The  customers  still  arrive  and  are 
served  according  to  an  exponential  distribution,  but  the 
different  types  of  customers  are  allowed  to  have  different 
mean  service  times.  The  equations  of  state  probabi 1 i ties 
cannot  be  solved  directly  for  such  a  model,  but  they  can  be 
approximated  by  a  finite  state  space  and  solved  iteratively 
by  a  numerical  technique  such  as  Gauss-Seidel  iteration 
<4: 123, 158-165) . 

There  are  also  many  applications  where  the  service  time 
is  not  exponential,  but  can  be  approximated  by  a  sum  or 
weighted  average  of  exponential  distributions.  In 
mathematical  terms,  the  service  time  has  an  Erlang  or  a 
hyperexponential  probability  distribution.  Such  an 
approximation  is  valid  for  many  situations  since  the  analyst 
can  control  the  variance  of  the  model  service  time  relative 
to  the  mean  of  the  model  service  time  by  the  selection  of  an 
appropriate  sum  or  weighted  average.  Cooper  develops  a 


solution  to  such  a  model;  the  solution  is  called  the  method 
of  phases  <4:171-175).  This  technique  allows  the 
transformation  of  a  general  service  time  model  into 
essentially  an  exponential  service  time  model.  In  this 
technique,  the  states  represent  not  the  number  of  customers 
in  the  system,  but  rather  the  number  of  exponential  phases 
yet  to  be  completed  in  the  system.  This  technique  may  also 
be  applied  to  cases  where  the  time  between  customer  arrivals 
is  not  exponential. 

It  is  possible  to  analyze  multiple  queue  systems  called 
networks  in  which  a  customer  requires  service  at  more  than 
one  queueing  station.  Customers  enter  the  system  at  various 
points,  queue  for  service,  and  then  proceed  to  some  other 
point  in  the  network  to  receive  additional  service. 

Kleinrock  shows  that  if  each  queue  has  exponential  service 
times,  and  if  the  time  between  arrivals  (from  outside  the 
network)  at  each  queue  is  also  exponential,  and  if  the 
customers  move  from  one  queue  to  another  according  to  a 
specified  probability  distribution,  then  each  queue  can  be 
separately  analyzed  and  solved  by  the  simple  exponential 
model  (5:147-161).  A  special  case  of  this  is  the  closed 
network,  which  has  a  fixed  and  finite  number  of  customers. 
Customers  circulate  throughout  the  system  but  never  leave 
it,  and  no  new  customers  may  enter  the  system  (6:14-15). 

If  the  analysis  of  the  network  can  be  broken  down  to 
analysis  of  the  individual  queues,  the  solution  of  the 
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network  is  relatively  inexpensive.  However,  there  are  many 
applications  which  cannot  be  solved  by  this  technique  since 
they  do  not  meet  the  restrictive  assumptions.  Sauer  and 
Chandy  discuss  many  of  these  applications  which  include 
simultaneous  resource  possession,  priority  queue  disciplines 
(where  the  customers  are  not  served  first-in,  first-out), 
and  other  variations  of  the  queueing  network.  Sauer  and 
Chandy  discuss  techniques  which  can  be  used  to  obtain 
numerical  approximations  to  these  applications  <2:26-30). 

An  even  more  general  model  has  been  developed  by 
Basket,  et  al ,  which  allows  for  many  different  classes  of 
customers  in  the  network.  The  model  also  allows  for  several 
different  service  disciplines,  and  for  a  broad  class  of 
service  time  distributions.  The  model  was  developed  to 
describe  a  computer  system  as  a  network  of  processors  and  a 
collection  of  customers  (jobs  and  tasks).  Different  service 
centers  may  have  different  queueing  disciplines.  Different 
customers  may  have  different  routes  through  the  network; 
customers  may  also  change  from  one  class  to  another  when 
changing  service  centers.  The  solution  calculates  the 
marginal  distribution  of  the  number  of  customers  at  each 
queue  (7:248-260).  Cox,  as  an  AFIT  thesis,  has  written  a 
FORTRAN  computer  program  which  solves  the  model  of  Basket, 
et  al .  His  computer  program  uses  special  algorithms  which 
are  computationally  efficient  for  large  networks 
(8:3-3,3-26) . 
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Traditionally,  the  most  common  technique  used  to  model 


avionics  maintenance  is  computer  simulation.  Clark  and 
Allen  have  developed  a  Q-GERT  simulation  program  which  is  a 
model  of  F-lll  avionics  and  automatic  test  equipment 
maintenance.  The  model  considers  not  only  the  -failures  and 
repairs  o-f  avionics  LRUs,  but  also  the  -failures  and  repairs 
o-f  the  ATE  itself.  The  model  also  considers  the  impact  of 
scarce  resources  such  as  maintenance  technicians  or  spare 
parts  (12:394-397).  The  F-16  prime  contractor,  General 
Dynamics,  has  also  developed  a  simulation  model  of  the 
avionics  maintenance  workload  on  ATE.  The  General  Dynamics 
model,  called  the  Station  Loading  Model,  is  a  FORTRAN 
computer  program  which  uses  Monte  Carlo  (random  number) 
techniques  to  simulate  the  arrivals  and  repair  of  avionics 
LRUs.  The  model  is  very  detailed  in  the  modeling  of  the 
maintenance  of  the  LRUs  and  has  separate  inputs  for  set  up, 
ATE  warm  up,  primary  test  time,  and  other  maintenance  task 
times.  General  Dynamics  has  also  written  a  Station  Loading 
Model  User’s  Manual  which  explains  the  inputs,  logic,  and 
outputs  of  the  model  (13). 

Three  AFIT/LS  students  (Bryson,  Husby,  and  Webb)  in 
1982  wrote  a  joint  thesis  on  a  simulation  model  of  F-1& 
avionics  maintenance.  The  F-16  research  effort  was  similar 
in  purpose  to  the  current  B-1B  research  effort,  and 
similarly  involved  the  development  of  a  Q-GERT  simulation 
model.  This  previous  research  effort  was  of  particular 


/  // 


$>* 


n .rir-7  vv  vvvvy.vv.T 


W 


probability  distributions  may  be  used  in  a  simulation  model. 
However,  Cooper  points  out  many  disadvantages  to  a 
simulation  approach.  First,  there  is  always  some 
uncertainty  as  to  when  equilibrium  (steady  state)  is 
reached.  Second,  simulation  introduces  a  sampling  error 
since  the  model  results  are  based  on  only  a  -finite  number  o-f 
events.  Finally,  an  error  is  introduced  by  the  random 
number  generator.  Cooper  concludes  that  simulation  should 
only  be  used  when  mathematical  analysis  is  not  -feasible. 

Queueing  theory  is  a  useful  analytical  technique  which 
can  be  used  to  study  many  situations  involving  customers 
(avionics  LRUs)  waiting  for  service  (maintenance).  The 
actual  queueing  model  developed  in  this  research  effort 
involved  the  combination  of  three  techniques.  Specifically, 
the  solution  technique  to  be  described  later  uses 
Kleinrock's  discussion  of  queueing  networks,  Cooper's 
discussion  of  the  method  of  phases  approximation,  and 
Cooper's  numerical  approach  for  problems  with  classes  of 
customers  with  different  mean  service  times.  However,  the 
actual  analytical  queueing  model  still  required  many 
simplifying  assumptions,  and  the  results  obtained  were  not 
entirely  satisfactory.  The  simulation  model  could 
accomodate  a  more  realistic  (and  complex)  approach.  The 
simulation  model  developed  for  this  research  effort  used 
many  features  taken  liberally  from  the  Clark  and  Allen  model 
and  the  General  Dynamics  Station  Loading  Model. 
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II. 


B-1B  Avi.oni.cg  Mai.ntenan£e 


B-1B  Program  Descri.eti.gn 

The  B-1B  is  a  new  multi-role  bomber  which  will 
eventually  replace  the  B-52  as  the  penetrating  bomber 
element  of  the  strategic  TRIAD.  The  B-1B  also  has 
collateral  missions  as  a  conventional  bomber  or  as  a  cru 
missile  launch  platform.  The  B-1B  could  potentially  be 
for  naval  and  theatre  conventional  warfare,  or  for  theat 
nuclear  warfare.  Current  plans  call  for  the  production 
deployment  of  a  force  of  100  B-1B  aircraft  assigned  to  f 
main  operating  bases  (MOBs) .  Delivery  of  26  aircraft  to 
#1  will  start  in  June  1985,  delivery  of  16  aircraft  to  Mi 
#2  will  start  in  September  1986,  delivery  of  32  aircraft 
MOB  #3  will  start  in  January  1987,  and  delivery  of  16 
aircraft  to  MOB  #4  will  start  in  October  1987.  Each  B-li 
squadron  consists  of  16  Primary  Aircraft  Authorizations 
(PAA) .  MOB  #1  also  has  an  additonal  10  aircraft  for  a 
Combat  Crew  Training  Squadron  (CCTS) .  MOB  #3  will  be  th 
base  for  two  squadrons.  It  should  be  noted  that  the  PAA 
total  assigned  to  the  MOBs  consists  only  of  90  aircraft; 
this  allows  for  a  reserve  of  10  aircraft  to  allow  for 
aircraft  attrition  or  airframe  maintenance  at  a  depot 


The  B-1B  has  an  integrated  avionics  system  totaling 
over  424  installed  line  replaceable  units  (LRUs)  of  which 
there  are  approximately  212  repairable  LRUs.  Tentatively, 
109  LRUs  have  been  designated  for  base  level  repair  on  the 
B-1B  automatic  test  equipment  (ATE).  Other  repairable  LRUs 
are  designated  for  base  level  repair  on  other  support 
equipment  or  for  depot  level  repair.  The  B-1B  avionics 
consist  of  offensive  avionics,  defensive  avionics,  and 
miscellaneous  avionics  associated  with  other  systems.  The 
B— IB  offensive  avionics  uses  updated  B-52  Offensive  Avionics 
System  (OAS)  equipment  as  well  as  a  terrain  following  radar 
and  a  new  inertial  navigation  system.  The  B-1B  defensive 
avionics  consists  of  a  complex  electronic  countermeasures 
(ECM)  system  known  as  the  ALQ-161.  A  smart  jamming 
enhancement  and  a  tail  warning  function  have  also  been 
included  on  the  B-1B.  Another  important  avionics  system  is 
the  Central  Integrated  Test  System  (CITS).  CITS  is  an 
automated  system  that  performs  fault  isolation  and 
verification  of  system  performance  in  real  time.  It 
provides  information  both  to  the  crew  during  flight  and  on 
the  ground  for  maintenance  debriefing  and  troubleshooting. 

It  is  on  the  CITS  data,  for  the  most  part,  that  the 
organzational  level  maintenance  actions  are  based  (16:1-1). 


QC9§0iS.*ti.9D*l  Level.  Maintenance 

Organizational  level  maintenance  consists  of  those 
tasks  normally  performed  on-aircraft  (on  the  flight  line)  by 
SAC  maintenance  technicians.  For  the  avionics  LRUs,  this 
maintenance  consists  of  debriefing,  fault  isolation,  removal 
and  replacement  of  the  failed  LRU,  and  clean-up  tasks  such 
as  documentation  (16:11-2).  The  crew  debriefing  is  at  the 
end  of  each  flight.  Typically,  the  B-1B  air  crew 
communicates  discrepancies  to  the  maintenance  team  chief  on 
the  ground  while  the  aircraft  is  returning  to  base.  The 
maintenance  chief  then  calls  the  appropriate  maintenance 
personnel  to  the  crew  debriefing  room.  When  the  crew 
arrives,  each  crew  member  is  matched  with  his  maintenance 
counterpart  to  discuss  the  malfunctions  and  to  agree  on  the 
description  of  the  discrepancy  which  will  then  be  written  in 
various  maintenance  documents.  Maintenance  will  then  be 
performed  to  remedy  each  discrepancy  (17).  For  avionics 
LRUs,  it  is  usually  first  necessary  for  the  maintenance 
technicians  to  isolate  the  cause  of  the  discrepancy  to  a 
single  LRU.  For  example,  a  reported  radar  system 
discrepancy  may  be  found  to  be  due  to  a  failed  radar 
antenna.  The  CITS  information  will  be  used  by  the 
technicians  for  the  fault  isolation.  Whenever  the  CITS 
information  is  ambiguous  or  incomplete,  the  maintenance 
technicians  would  have  to  engage  in  manual  troubleshooting 


to  find  the  failed  LRUs.  Once  the  fault  is  found,  the 
maintenance  technicians  remove  the  failed  LRU  from  the 
aircraft  and  replace  it  with  a  good  LRU  <obtained  from  base 
supply)  of  the  same  type.  The  technicians  then  perform  a 
system  verification  to  ensure  that  the  discrepancy  has  been 
resolved.  Once  the  technicians  complete  the  maintenance, 
they  then  fill  out  all  maintenance  documentation  and  also 
take  the  failed  LRU  to  a  production  control  point  of  the 
Avionics  Maintenance  Squadron  <AMS>  at  intermediate  level. 

Organizational  level  maintenance  is  performed  by 
maintenance  technicians  from  the  Organizational  Maintenance 
Squadron  (OMS) .  The  OMS  technicians  not  only  remove  and 
replace  LRUs,  but  also  service  the  aircraft  and  conduct 
various  pre-flight  and  post -flight  inspections.  Each 
aircraft  has  a  dedicated  maintenance  team  headed  by  a  team 
chief  (usually  a  Master  Sergeant).  The  OMS  technicians  are 
not  avionics  specialists,  and  have  no  intermediate  level 
duties.  For  this  reason,  they  might  occassional ly  call  for 
special  assistance  from  the  AMS  technicians  if  required 
(21) . 

Intermediate  Level  MglQtftQSQSS 

The  intermediate  level  shops,  which  are  located  at  base 


level,  house  the  specialist  activities  that  maintain  and 


support  the  aircraft  systems.  One  of  the  most  important 
organizations  at  intermediate  level  is  the  AMS,  which  is 
responsible  for  the  repair  of  the  failed  avionics  LRUs.  The 
repair  of  LRUs  designated  for  maintenance  on  the  B-1B  ATE  is 
an  AMS  responsibility,  although  the  AMS  has  other  duties  as 
wel  1 . 

The  scheduling  of  avionics  LRU  repair  is  accomplished 
approximately  on  a  first-come,  first-seree  basis.  However, 
there  might  be  an  adjustment  tm~ allow  for  a  higher  priority 
for  LRUs  with  a  low  «pare  asset  posture.  Once  scheduled,  an 
AMS  technician  repairs  the  LRU  on  the  ATE  test  station  of 
t+ie  appropriate  type.  In  most  cases,  the  technician  is  able 
to  successfully  repair  the  LRU  at  intermediate  level!  such  a 
maintenance  action  is  designated  as  Repairable  This  Station 
(RTS).  For  an  RTS  maintenance  action,  the  repair  takes 
place  in  a  complicated  sequence  of  events.  First,  the 
technician  must  set  up  the  LRU  on  the  ATE  test  station.  The 
LRU  is  physically  connected  to  the  station  via  an  interface 
test  adapter  (ITA).  The  technician  then  runs  a  performance 
test  until  the  fault  in  the  LRU  is  found.  The  fault  is 
usually  in  one  of  the  printed  circuit  boards  of  the  LRU. 

The  printed  circuit  boards  are  known  as  shop  replaceable 
units  (SRUs) .  Once  the  bad  SRU  is  identified,  the 
technician  goes  to  supply  to  get  a  good  replacement  SRU.  A 
small  number  of  high  failure  rate  SRUs  are  kept  in  a  forward 
supply  point!  other  SRUs  are  obtained  from  base  supply.  In 


any  event,  the  technician  removes  the  bad  SRU  from  the  LRU 
and  replaces  it  with  the  good  SRU.  The  technician  then  runs 
one  complete  performance  test  to  verify  the  success  of  the 
repair.  The  repaired  LRU  is  then  removed  from  the  station 
and  taken  to  base  supply.  The  technician  also  fills  out  all 
required  maintenance  documentation.  Repairs  may  not  always 
be  this  simple,  however.  The  performance  test  may  not 
always  be  able  to  fault  isolate  to  a  single  SRU.  Rather, 
the  test  might  only  fault  isolate  to  a  group  pf  SRUs,  or  may 
indicate  ambiguous  results.  In  such  a  case,  the  technician 
resorts  to  manual  troubleshooting  to  precisely  locate  the 
fault.  Another  complication  is  that  not  all  maintenance 
actions  lead  to  actual  repairs.  One  possibility  is  that  the 
LRU  has  failed  in  a  way  which  is  beyond  the  capability  of 
the  AMS  to  repair.  For  example,  if  the  failure  has  occurred 
in  a  chassis  of  the  LRU,  and  not  a  removeable  SRU,  the  LRU. 
would  then  be  sent  to  depot  for  repair.  Such  a  maintenance 
action  is  called  Not  Repairable  This  Station  (NRTS) .  The 
LRU  is  then  sent  to  depot  (for  greater  facilities  and  higher 
skill  level  technicians)  for  repair.  Another  possibility  is 
that  the  AMS  technician  cannot  find  any  problem  with  the 
LRU.  This  might  occur  if  the  LRU  was  unnecessarily  removed, 
or  if  there  is  some  incompatibility  between  the  CITS  and  the 
ATE.  This  type  of  mainteance  action  is  called  Retest  Okay 
(RTOK) .  The  LRU  is  taken  to  base  supply  without  any  repair 
being  required. 


AMS  technicians  that  repair  LRUs  on  ATE  can  come  from 
one  of  four  branches.  Each  AMS  technician  has  an  Air  Force 
Specialty  Code  (AFSC)  which  describes  his  particular 
specialty.  There  is  a  dedicated  AFSC  for  offensive  avionics 
LRUs,  for  defensive  avionics  LRUs,  for  communication  and 
navigation  LRUs,  and  for  automatic  flight  control  and 
instrument  LRUs.  There  is  no  cross-uti 1 ization  of  personnel 
between  branches. 

All 

Not  only  do  LRUs  fail  and  require  maintenance,  but  so 
do  the  ATE  test  stations  themselves.  The  test  stations, 
being  automatic  test  equipment,  are  complex  electronic 
devices.  When  a  station  fails,  a  technician  must  be  called 
to  repair  the  station.  The  technician  actually  uses  the 
station  itself  as  part  of  the  repair  process.  The 
technician  must  first  isolate  the  station  fault  to  a  test 
replaceable  unit  <TRU) .  A  TRU  is  the  station  equivalent  of 
an  LRU.  The  technician  must  then  further  isolate  the  fault 
to  a  bad  SRUi  the  technician  removes  and  replaces  the  failed 
SRU  in  a  manner  similar  to  LRU  repair.  Station  maintenance 
is  accomplished  by  technicians  from  the  Precision 
Measurement  Equipment  Laboratory  (PMEL)  branch  of  the  AMS. 


III. 


D$t a  Collection  and  Gathering 
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S=l§  Sgur ce  Selection 


The  B-1B  prime  contractor,  Rockwell  International, 
selected  the  B-1B  ATE  vendor  by  a  competitive  source 
selection.  This  was  a  source  selection  conducted  by 
Rockwell  and  not  by  the  Air  Force.  As  part  of  the  source 
selection  process,  Rockwell  International  required  each 
competing  potential  vendor  to  prepare  a  station  loading 
analysis  (which  is  logistics  terminology  tor  an  avionics/ATE 
queueing  analysis).  To  allow  each  vendor  the  opportunity  to 
per -Form  such  an  analysis,  Rockwell  provided  each  vendor  with 
detailed  reliability  and  maintainability  data  on  all  of  the 
avionics  LRUs  designated  for  repair  on  the  ATE.  This  data, 
with  two  exceptions  that  will  be  explained  later,  became  the 
baseline  for  the  analysis  conducted  in  the  next  two 
sections.  This  reliability  and  maintainabi 1 i ty  data,  which 
was  included  in  the  Rockwell  Request  for  Proposal  (RFP) ,  is 
included  in  the  computer  listing  included  in  Appendix  A.  An 
example  of  this  data  is  shown  in  Table  I. 
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Table  I.  Sample  Reliability/Maintainability  Data 


QTY  AR 

wyc  nomenclature  assqq  a/c  rate 


RETEST 
*/.  OK 


58AAA  Avionics  Processor  BM 


124.0  20.0  16.0 


6.671 


26.0 


The  term  WUC  refers  to  a  -five  character  work  unit  code 
designation  for  each  LRU.  This  particular  LRU  is  called  an 
avionics  processor.  The  LRU  is  provided  by  Boeing ,  one  of 
the  three  B-1B  associate  avionics  contractors.  There  are 
eight  avionics  processors  installed  on  every  B-1B  aircraft. 
The  term  AR  RATE  refers  to  the  arrival  rate  of  this 
particular  LRU  into  the  AMS.  The  arrival  rates  included  in 
the  Rockwell  RFP  were  estimated  by  the  associate  contractors 
before  March  of  1983.  These  values  were  superseded  by  more 
recent  contractor  estimates  that  were  available  during  this 
research  effort,  as  explained  in  the  next  subsection.  The 
term  RETEST  7.  OK  means  for  this  LRU  that  25*/.  of  the 
maintenance  actions  will  be  RTOKs.  The  term  TP  refers  to 
performance  test  time.  The  avionics  processor  requires  124 
minutes  to  run  one  complete  performance  test.  This  is  the 
second  term  in  the  original  RFP  which  has  been  adjusted. 
Rockwell  International  estimated  the  LRU  performance  test 
times  for  a  hypothetical,  generic  ATE  design.  Rockwell 
estimates  presumably  reflect  the  level  of  contemporary  ATE 
technology.  Immediately  after  contract  award,  however,  the 


winning  ATE  vendor  (General  Dynamics  Electronics)  suggested 
that  their  ATE  design,  due  to  expanded  ATE  computer 
capabilities,  could  achieve  a  dramatic  6>QV.  reduction  in  the 
LRU  performance  test  times.  This  means,  for  example,  that 
GDE  estimates  that  the  avionics  processor  performance  test 
time  would  only  be  49.6  minutes.  The  GDE  estimates  were 
adopted  as  the  baseline  for  this  research  effort.  The 
impact  of  this  reduction  is  described  in  the  discussion  on 
sensitivity  analysis  included  in  Section  VI.  The  term  TS  is 
the  set-up  time,  also  in  minutes,  required  to  initially 
set-up  the  LRU  on  the  ATE  test  station.  The  term  TD  refers 
to  the  time  in  minutes  for  the  technician  to  tear — down  the 
LRU  from  the  ATE  after  the  maintenance  action  is  completed. 
Finally,  the  term  TR  refers  to  the  time  to  repair  the  LRU 
(remove  and  replace  an  SRU)  also  in  minutes. 

In  the  Rockwell  instructions  provided  in  the  RFP,  the 
furnished  maintainability  data  were  expected  to  be  used  to 
compute  the  mean  task  times  for  both  LRU  RTS  actions  and  LRU 
RTOK  actions.  These  means  would  be  computed  as  follows: 


(MEAN  REPAIR  TIME, 

RTS) 

=  TS 
+  TR 

+ 

+ 

(0. 

TP 

5)  TP 
+  TD 

(1) 

(MEAN  TASK  TIME, 

RTOK) 

=  TS 

+ 

TP 

+  TD 

(2) 

Note  that  the  time  for  fault  isolation  is  assumed,  on  the 
average,  to  be  one  half  the  time  of  one  complete  performance 
test.  In  other  words,  the  time  of  fault  isolation  is 


assumed  to  be  distributed  uniformly  over  the  first 
performance  test.  Interviews  with  SAC  maintenance  and 
logistics  personnel  (21)  revealed  many  deficiencies  in  the 
Rockwell  approach.  First,  it  failed  to  account  for  the 
travel  time  to  obtain  an  SRU  spare  for  a  repair  (the 
technician  who  obtains  the  SRU  spare  is  the  same  individual 
accompl ishing  the  repair).  Second,  the  mean  repair  time 
assumes  that  the  fault  isolation  is  always  successful,  and 
does  not  allow  for  the  possibility  of  extended 
troubleshooting.  Third,  the  equation  for  the  mean  task  time 
of  a  RTOK  maintenance  action  excludes  the  possibility  of  a 
second  performance  test.  In  practice,  however,  the 
technician  may  feel  the  need  for  a  second  performance  test 
if  he  feels  that  the  first  test  was  not  completely  correct 
and  unambiguous.  Fourth,  and  finally,  the  Rockwell 
equations  do  not  consider  the  possibility  of  a  NRTS 
maintenance  action.  Refinements  to  these  equations  are 
proposed  in  the  next  subsection. 

Reliability  and  Maintainability  Estimates 

LRU  reliability  estimates  from  each  of  the  B-1B 
associate  avionics  contractors  are  regularly  furnished  to 
the  B-1B  SPO.  The  most  recent  estimates  (as  of  August  1983) 
formed  the  baseline  for  this  research  effort.  Reliability 


in  a  logistics  sense  (as  opposed  to  an  engineering  sense)  is 
measured  in  Mean  Time  Between  Demand  (MTBD) .  Demand  refers 
to  a  demand  on  supply  (which  occurs  when  an  LRU  is  removed 
from  the  aircraft  for  corrective  maintenance) .  This 
includes  not  only  true  (inherent)  failures,  but  also 
includes  induced  failures  and  RTDKs.  The  time  in  MTBD  is 
measured  in  aircraft  flying  hours  (and  not  LRU  operating 
hours).  The  estimates,  by  LRU,  are  included  in  the  computer 
listing  included  in  Appendix  A.  In  addition  to  the 
contractor  estimates,  MTBD  estimates  were  also  obtained  from 
HQ  SAC  maintenance  and  logistics  personnel  (21).  The  SAC 
estimates  were  typically,  but  not  always,  significantly  more 
pessimistic  than  the  contractor  estimates.  The  SAC 
estimates  were  based  on  the  current  experience  of  like  LRUs 
on  the  B-52,  FB-111,  and  other  systems. 

LRU  maintainability  estimates  were  obtained  from  the 
Rockwell  International  ATE  RFP  as  explained  in  the  last 
subsection.  However,  the  lower  (more  optimistic)  0DE 
performance  test  times  were  used  as  the  baseline  of  this 
research  effort.  In  addition,  the  algorithms  to  compute  the 
mean  task  times  were  modified  as  follows: 

(MEAN  REPAIR  TIME,  RTS)  =  TS  +  (0.5) TP  +  (.75>TR 

+  ( . 25) (2) TR  +  TP  +  TD  (3) 

(MEAN  TASK  TIME,  RTOK)  -  TS  +  TP  +  (0.3) TP  +  TD  (4) 

(MEAN  TASK  TIME,  NRTS)  -  TS  +  (0.5) TP  +  TD  (5) 


These  modifications  were  developed  based  on  discussions  with 
HQ  SAC  maintenance  and  logistics  personnel  (21).  For  the 
LRU  repairs,  the  fault  isolation  to  a  single  SRU  by  the 
performance  test  is  assumed  to  be  successful  only  75%  of  the 
time.  For  the  other  25%,  the  repair  times  are  doubled  to 
account  for  manual  fault  isolation.  For  the  LRU  RTOKs,  it 
is  assumed  that  a  second  performance  test  is  run  30%  of  the 
time.  The  actual  percentage  would  depend  on  the  degree  of 
confidence  that  the  technician  has  in  the  ATE  software 
(performance  test).  A  third  mean  task  time  has  been  added 
for  NRTS  maintenance  actions.  It  is  assumed  that  the 
technician  can  successfully  fault  isolate,  but  that  no 
repair  at  base  level  is  possible.  Therefore,  the  LRU  is 
removed  from  the  station  to  be  sent  to  a  depot  for  repair. 

In  addition  to  the  calculation  of  the  above  mean  task 
time  estimates,  the  models  of  the  next  two  sections  also 
include  an  allowance  for  time  to  obtain  a  SRU  spare.  This 
only  applies  to  repairs  (RTS  actions),  and  not  to  RTOK  or 
NRTS  actions.  The  travel  time  to  the  forward  supply  point 
is  estimated  at  15  minutes  round  trip,  and  the  time  to 
obtain  a  SRU  spare  from  base  supply  is  estimated  at  one 
hour.  It  is  also  estimated  that  a  SRU  spare  will  be 
available  on  the  forward  supply  point  25%  of  the  time,  and 
will  be  available  at  base  supply  the  remaining  75%  (21). 


§!ri§  i_onal_  Data 

As  described  earlier,  the  building  block  of  the  B-1B 
fleet  is  the  16  PAA  squadron.  A  typical  weekly  flying 
schedule  for  a  16  PAA  squadron  was  obtained  from  HQ  SAC 
personnel  <15>.  In  any  given  week,  only  8  aircraft  will 
actually  fly  daily  training  missions.  The  other  8  aircraft 
will  have  alert  obligations  or  will  be  undergoing  scheduled 
(phase)  inspections.  The  S  aircraft  that  will  fly  missions 
will  typically  fly  a  total  of  21  sorties  in  a  week,  leading 
to  107  flight  hours  in  a  week.  A  significant  aspect  of  the 
flying  schedule  is  that  the  sorties  are  not  spread  evenly 
during  the  days  of  the  week.  Typically,  3  sorties  will  be 
flown  on  Monday,  6  sorties  on  Tuesday  and  Wednesday,  4 
sorties  on  Thursday,  and  2  sorties  will  be  flown  on  Friday. 
In  addition,  for  any  day  of  the  week,  the  sorties  are  not 
spread  evenly  over  the  day.  Sorties  will  typically  be 
launches  in  the  morning  or  in  the  late  evening.  What  this 
means,  from  the  point  of  view  of  a  modeling  strategy,  is 
that  the  arrival  rate  of  the  avionics  LRUs  will  not  be 
constant  over  the  maintenance  day  or  during  different  days 
of  the  week. 

In  addition  to  B-1B  flying  hour  data,  information  was 
also  obtained  (21)  on  maintenance  technician  manning  levels 
These  values  represent  maintenance  manpower  authorizations, 


and  are  shown  in  Table  II.  Values  are  shown  for  the  16 
squadron  and  the  10  PAA  CCTS. 

Table  II.  Maintenance  Manpower  Authorizations 


(Actual 

(Assumed 

Di 

Author i z  at i ons ) 

Labor /Shi 

Qa£®9Q£¥ 

16  PAA 

(CCTS) 

10  PAA 

16  PAA 

(C 

iQ 

Organizational  Maintenance 

Squadron  (OMS) 

128 

78 

38 

Communication  Navigation 

Branch,  AMS 

22 

10 

7 

Auto  Flight  Control / 

Instrument  Branch,  AMS 

14 

6 

4 

Offensive  Avionics 

Branch,  AMS 

55 

18 

16 

Defensive  Avionics 

Branch,  AMS 

35 

14 

io 

Precision  Measurement 
Equipment  Laboratory 
Branch,  AMS 

18 

0 

5 

These  manpower  authorizations  represent  total  slots 
available  to  each  squadron.  As  a  first  approximation,  t 
maintenance  policy  will  call  for  two  8-hour  shifts  of 
maintenance  per  work  day  (Monday  through  Friday).  Manpo 
authorizations  will  then  be  divided  evenly  between  the  t 
shifts.  As  a  better  approximation,  the  technicians  may 
perform  maintenance  during  the  "graveyard"  shift  on  an 
occasional  basis  as  workload  requirements  dictate.  In 
addition,  a  certain  number  of  OMS  technicians  would  have 
be  available  during  the  "graveyard"  shift  to  launch  and 


service  aircraft.  For  modeling  purposes,  the  authorizations 
were  divided  evenly  between  the  two  shifts.  Moreover,  only 
607.  of  the  technicians  were  assumed  to  be  available  for 
direct  labor.  The  other  40 7.  would  account  for  illness, 
leave,  or  duties  outside  the  scope  of  this  research  effort. 

Goodness-gf -Fi.t  Analysis 

The  baseline  inputs  used  in  this  research  effort  were 
reliability  and  maintainability  estimates  obtained  from  the 
B— IB  associate  contractors.  However,  the  estimates 
represent  mean  values,  and  do  not  specify  the  distributional 
nature  of  the  reliability  and  maintainability  random 
variables.  To  select  the  best  distributional  assumptions 
for  subsequent  model  development,  operational  data  was 
obtained  from  the  F-16  program  and  the  B-52  OAS  program  for 
goodness-of-f i t  analysis. 

The  first  factor  considered  was  avionics  reliability. 
Reliability  in  this  context  means  Mean  Time  Between  Demand 
(MTBD)  and  not  the  usual  MTBF.  It  would  be  highly  desirable 
to  perform  goodness-of-f it  tests  on  the  LRU  inter — arrival 
times  to  obtain  the  best  distributional  approximation. 
However,  actual  operational  maintenance  data  is  not 
collected  in  this  manner.  Maintenance  technicians  record 
the  time  of  the  removal  of  the  avionics  LRU  from  the 
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aircraft;  it  is  not  possible  to  reconstruct  the  actual 
inter-arri val  times  of  the  LRU  failures.  For  this  reason, 
it  was  only  possible  to  test  if  the  LRU  arrival  process 
could  be  approximated  by  a  Poisson  process.  Another  issue 
in  testing  the  distribution  of  the  LRU  arrivals  is  the 
measure  of  time.  Time  could  either  be  measured  in  flying 
hours  or  in  calendar  time.  These  two  approaches  would  not 
be  equivalent  since  the  flying  hours  per  base  per  month  will 
not  be  perfectly  constant;  the  variation  in  the  flying  hours 
would  add  to  the  variability  of  the  LRU  arrival  process. 
Rather  than  derive  a  distribution  for  the  flying  hours  and 
derive  a  distribution  for  the  arrival  process  (in  flying 
hours),  it  was  simpler  to  fit  a  distribution  to  the  LRU 
a -rival  process  measured  in  calendar  time.  In  the  models 
discussed  in  the  next  two  sections,  flying  hours  per  base 
per  month  are  treated  as  a  fixed  constant  and  not  as  a 
random  variable.  The  variability  in  the  flying  hours  is 
reflected  in  the  distribution  of  the  arrival  process.  In 
any  event,  the  key  question  was  if  the  LRU  arrival  process 
could  be  approximated  by  a  Poisson  process. 

Operational  maintenance  data  on  19  major  F-16  avionics 
LRUs  were  obtained  from  the  F-16  SPO  Centralized  Data  System 
(CDS).  All  of  these  19  LRUs  are  repaired  on  the  F-16  ATE. 
These  LRUs  were  examined  on  an  aggregate  basis,  and  also  two 
LRUs  were  examined  on  an  individual  basis.  A  chi-square 
goodness-of-f it  was  used  to  test  the  null  hypothesis  that 


the  LRU  arrival  process  was  a  Poisson  process.  Failing  to 
reject  the  null  hypothesis  does  not  prove  that  the  LRU 
arrival  process  is  truly  Poisson;  it  merely  suggests  that  it 
should  be  a  reasonable  approximation.  In  examining  the  F-16 
reliability  figures,  two  concerns  become  quickly  evident. 
First,  in  general,  the  F-16  avionics  reliability  had 
significantly  improved  over  time  since  the  initial  F-16 
deployment.  To  overcome  this  concern,  the  data  used  in  the 
chi-square  goodness-of-f i t  test  reflected  only  the  last  six 
months  of  F-16  experience,  thereby  capturing  a  mature, 
steady-state  situation.  It  should  be  pointed  out, 
therefore,  that  the  models  described  in  the  next  two 
sections  should  only  be  used  to  estimate  the  workload  during 
mature,  steady-state  experience.  Second,  the  F-16 
reliability  varied  considerably  from  one  base  to  another. 
This  is  in  part  due  to  configuration  differences,  and  in 
part  due  to  differences  in  the  quality  of  maintenance 
documentation.  To  estimate  the  variability  at  a  single 
base,  the  data  used  in  the  test  was  generated  by  a  single 
base,  Nellis  Air  Force  Base.  For  the  data  for  a  mature  six 
month  period  from  a  single  F-16  base,  it  was  concluded  that 
the  Poisson  process  would  be  a  reasonable  approximation  to 
the  LRU  arrival  process.  This  was  true  for  the  arrivals  of 
the  two  individual  LRUs  tested,  and  it  was  true  for  the 
aggregate  arrivals  of  the  19  LRUs  considered.  Another  point 
that  needs  to  be  mentioned  is  that  the  time  interval 


selected  (for  each  data  point  in  the  test)  was  one  week. 
This  would  then,  of  course,  smooth  over  any  possible 
differences  between  days  of  the  week  or  different  times  of 
the  day.  The  statistical  results  are  presented  in  Table 
III. 

Table  III.  Results  of  Chi-Square  Goodness-of -Fi t  Tests 
Null  Hypothesis:  LRU  Arrival  Process  is  Poisson  Process 
(Measured  on  a  Weekly  Basis) 

Chi-Square  Degrees 

Jest  Statistic  of  Freedom  Significance 

F-16  Inertial 

Navigation  Unit  2.32  3  >.250 

F-16  Low  Power 

RF  Unit  1.63  3  >.250 

F-16  Major 

Avionics  LRUs  2.56  4  >.250 

The  second  factor  considered  was  avionics  repair  time 
(or  other  task  time)  at  intermediate  level.  Again, 
operational  data  was  obtained  from  the  F-16  program.  For 
reasons  discussed  earlier,  the  data  was  obtained  for  a 
mature  six  month  period  from  a  single  base.  It  was  not 
possible  to  get  separate  data  on  set-up  times,  test  times, 
repair  times,  or  tear-down  times;  maintenance  documentation 
is  not  that  detailed.  It  was  only  possible  to  get  task 
times  for  the  overall  maintenance  action.  It  was  possible, 
however,  to  separately  analyze  repairs,  RTOKs,  and  NRTS 
actions.  As  will  be  seen,  the  mean  and  variance  are  quite 
different  for  different  types  of  maintenance  actions. 


The  task  time  Mas  -first  collected  for  a  single  F-16  LRU 


(the  low  power  RF  unit).  52  observations  were  obtained  for 
repairs  (RTS  actions),  and  19  observations  were  obtained  for 
RTOKs.  There  were  only  3  NRTS  actions  during  the  six  month 
period  which  is  not  sufficient  to  do  a  meaningful  test.  For 
the  repairs  and  the  RTOKs,  the  following  distributions  were 
tested  for  the  possible  distribution  of  the  maintenance  task 
times:  exponential,  Erlang,  and  lognormal.  An  Erlang  was 

used  instead  of  a  (general >  gamma  because  an  Erlang  is 
easier  to  simulate,  and  because  an  Erlang  distribution  can 
be  used  in  the  method  of  phases  technique  of  analytical 
queueing  theory.  It  would  have  been  necessary  to  resort  to 
a  (general)  gamma  only  if  all  Erlang  distributions  were  poor 
approximations.  The  results  for  the  RTS  maintenance  actions 
are  shown  in  Table  IV. 


Table  IV.  Results  of  Chi-Square  Goodness-of-Fi t  Tests 
Null  Hypothesis:  LRU  Maintenance  Task  Times  Have 
Specified  Distribution 
(Case  for  RTS  Maintenance  Actions) 


Specified 

Distribution 

Exponential 
Erlang  (k«2) 
Erlang  (k*3> 
Lognormal 


Chi-Square 

Statistic 

21.10 

4.40 

4.68 

11.21 


Degrees 
of  Freedom 

7 

6 

6 

6 


Significance 

<.005 

>.250 

>.250 

<.100 


For  an  Erlang  distribution,  the  ratio  of  the  mean  to  the 
standard  deviation  is  known  to  be  the  reciprocal  of  the 
square  root  of  the  k  parameter.  By  taking  the  sample  mean 


divided  by  the  sample  standard  deviation,  it  was  possible  to 
obtain  a  point  estimate  o-f  that  k  parameter.  The  test  was 
run  with  k=2  and  k=3  since  2  and  3  were  the  closest  integer 
values  to  that  point  estimate  of  2.6. 

Similar  results  were  obtained  for  the  RTOK  maintenance 
actions.  The  Erlang  (this  time  with  k=5)  was  judged  to  be  a 
reasonable  approximation.  This  time,  the  lognormal  was  also 
a  reasonable  approximation.  Of  course,  with  only  19  data 
points,  it  is  difficult  to  establish  a  meaningful  and 
powerful  test. 

At  this  point  two  assumptions  were  made.  First,  the 
Erlang  distribution,  with  the  right  k  value,  could  be  used 
to  approximate  the  maintenance  task  times.  Second,  since 
there  was  not  sufficient  data  to  perform  tests  on  the  NRTS 
maintenance  actions,  it  was  assumed  that  the  distribution  of 
a  NRTS  action  would  be  the  same  as  for  a  RTOK  action.  The 
next  step  was  to  look  at  a  second  F-16  LRU  (the  inertial 
navigation  unit)  and  also  a  B-52  OAS  LRU  (the  signal  data 
converter).  The  purpose  was  to  make  point  estimates  of  the 
coefficient  of  variation  (the  mean  divided  by  the  standard 
deviation)  to  see  if  there  were  any  patterns  in  k-values. 
This  data  is  presented  in  Table  V. 


Table  V. 


Estimates  of  Coefficient  of  Variation 


Estimate 


k- value 


F-16  Low  Power  RF  Unit 

-  Repair 

-  RTOK 


.622 

.436 


F-16  Inert  Nav  Unit 

-  Repair 

-  RTOK 


.745 

.473 


B-52  Signal  Data  Converter 

-  Repair 

-  RTOK 


.585 

.531 


The  B-52  OAS  data  was  not  large  enough,  unfortunately,  to 


conduct  a  meaningful  goodness-of-f i t  test.  There  were  only 


7  repairs  and  only  17  RTOK  actions.  It  should  also  be  noted 


that  the  actual  k-value  of  an  Erlang  distribution  must,  of 


course,  be  an  integer.  At  this  point  in  the  research 


effort,  it  was  assumed  that  the  RTS  actions  could  be 


approximated  by  an  Erlang  distribution  with  a  k-value  of  2, 


and  that  RTOK  and  NRTS  actions  could  be  approximated  by  an 


Erlang  distribution  with  k=5. 


IV.  Simulation  Model 


Model_  Overview 

A  detailed  Q-GERT  simulation  program  was  developed 
which  can  be  used  to  measure  the  workload  of  avionics 
maintenance  on  intermediate  level  automatic  test  equipment 
(ATE).  The  model  simulates  the  flow  of  avionics  line 
replaceable  units  (LRUs)  from  LRU  failure  to  LRU  repair  and 
return  to  base  supply.  This  flow  is  shown  in  Figure  1.  The 
model  also  simulates  the  failures  and  maintenance  of  the  ATE 
test  stations  themselves. 

The  model  begins  with  a  simulation  of  avionics  LRUs 
which  fail  on  B-1B  aircraft  while  in  flight.  Sorties  are 
assumed  to  be  launched  early  in  the  morning  and  at  night,  so 
the  arrival  rate  of  avionics  LRUs  is  not  constant  throughout 
the  day.  The  arrival  rate  is  also  different  for  different 
days  of  the  week  due  to  different  flying  schedules.  As  an 
LRU  fails,  the  model  assigns  various  attributes 
(character i sti cs)  to  that  LRU  which  describes  its  subsequent 
repair.  These  attributes  include  the  test  station  type 
requirement,  the  intermediate  level  technician  requirement, 
the  type  of  maintenance  action  required,  and  the  hours  that 
will  be  required  for  the  LRU  repair. 


37 


Figure  1 


Flow  of  Avionics  LRU/ATE  Repair 


Once  an  avionics  LRU  fails,  it  requires  maintenance 
from  the  Organi 2 ati onal  Maintenance  Squadron  (OMS) 
technicians.  This  maintenance  consists  of  debriefing, 
removal  of  the  failed  LRU  from  the  aircract,  and 
transportation  of  the  LRU  to  the  avionics  shop.  The  model 
treats  OMS  technicians  as  constrained  resources,  and 
therefore  failed  LRUs  must  wait  (or  queue)  for  an  available 
technician  until  maintenance  may  be  performed. 

Once  at  the  avionics  shop,  the  LRUs  queue  for 
intermediate  level  maintenance  at  the  Avionics  Maintenance 
Squadron  (AMS).  The  LRUs  require  both  a  technician  and  a 
test  station  before  maintenance  can  commence.  Each  of  the 
station  types  and  each  of  the  AMS  technician  types  are 
treated  as  constrained  resources  in  the  model. 

Not  only  do  avionics  LRUs  fail  and  require  maintenance 
but  the  ATE  test  stations  themselves  fail  and  require 
repair.  The  model  generates  ATE  failures,  and  for  each 
failure,  assigns  attributes  which  describe  the  subsequent 
repair.  The  model  assumes  that  ATE  repair  takes  precedence 
over  LRU  repair  since  an  ATE  test  station  must  be  in  good 
working  order  in  order  to  be  used  during  the  LRU  repair. 

The  key  outputs  of  the  model  are  resource  utilization 
and  LRU  base  repair  cycle  time.  Resource  utilization  for 
all  station  types  and  all  technician  types  is  included.  Th 
base  repair  cycle  time  is  the  time  from  LRU  failure  until 
return  (of  the  LRU)  to  base  supply.  It  is  this  second 


output  -feature  that  will  be  used  to  determine  test  station 
requirements.  The  Q-GERT  network  diagram  and  the  computer 
listing  -for  the  model  are  both  included  in  Appendix  A. 

Detailed  Model,  Description 

The  model  generates  avionics  LRU  failures  according  to 
the  B-1B  flying  hour  profile  at  each  base.  LRU  failures  are 
generated  as  a  Poisson  process  for  a  four  hour  period  twice 
each  day.  B-1B  aircraft  are  assumed  to  be  launched  in  the 
morning  from  0800  until  1200,  and  in  the  late  evening  from 
2000  until  2400.  The  B-1B  sortie  duration  is  assumed  to 
last  5.0  hours,  so  the  avionics  LRU  failures  are  generated 
from  1100  until  1500  and  from  0100  until  0500  for  each  day. 
Nodal  modification  is  used  to  start  and  stop  the  Poisson 
process  for  each  four  hour  period.  The  model  also  keeps 
track  of  days  of  the  week  since  the  flying  hour  profile  (and 
therefore  LRU  arrival  rate)  is  different  for  different  days 
of  the  week . 

In  order  to  generate  LRU  failures  according  to  a 
Poisson  process,  it  is  necessary  to  simulate  exponential 
interarrival  times  over  a  fixed  period  (in  this  case,  from 
1100  until  1500  and  from  0100  until  0500).  This  is 
accomplished  in  the  model  by  a  FORTRAN  user  function,  where 
the  numerical  value  assigned  to  the  user  function  is  the 


exponential  interarrival  times  between  successive  LRU 


•failures.  The  FORTRAN  user  function  also  defines  attribute 
values  for  each  LRU  failure  which  are  used  in  the  later 
portion  of  the  model.  Detailed  task  listings  of  the  FORTRAN 
initialization  subroutine  and  the  user  function  are  shown  in 
Figure  2  and  Figure  3.  The  i ni t i al i zat 1  on  subroutine  simply 
defines  all  variables  which  are  needed  in  the  user  function. 
Variable  information  is  transmitted  to  the  user  function 
through  use  of  FORTRAN  COMMON  statements.  The  FORTRAN  user 
function  then  computes  the  mean  interarrival  time  of 
avionics  LRU  failures  for  the  correct  day  of  the  week  and 
the  correct  time  of  day.  The  actual  LRU  failure 
interarrival  time  is  then  simulated  according  to  an 
exponential  distribution.  The  user  function  then  uses 
several  more  random  number  draws  to  determine  the 
characteri sties  of  the  LRU  maintenance  which  will  take  place 
later  in  the  model.  These  characteristics  are  stored  as 
attributes  of  the  failed  LRU. 

Once  the  avionics  LRUs  have  failed,  they  require 
maintenance  at  organizational  level  (i.e.,  at  the  flight 
line).  This  maintenance  is  performed  by  two  QMS 
technicians.  The  avionics  LRUs  must  wait  in  queue  until 
technicians  are  available.  The  maintenance  consists  of 
three  distinct  tasks  performed  in  series.  First,  the 
maintenance  crew  must  attend  a  debriefing  (with  the  aircrew) 


which  takes  a  constant  time  of  1.0  hour.  Second,  the 


Subroutine  UI 


Initialize  variables  by  data  statements 

Add  LRU  arrival  rates  to  obtain  aggregate  avionics  Mean 
Time  Between  Demand  (MTBD  in  -flying  hours) 

Compute  LRU  arrival  rate  as  percentage  of  aggregate 
avionics  arrival  rate  (reciprocal  of  MTBD) 

Add  test  station  failure  rates  to  obtain  aggregate  ATE 
Mean  Time  Between  Maintenance  (in  operating  hours) 

Compute  test  station  failure  rate  as  percentage  of 
aggregate  ATE  failure  rate 


Figure  2.  FORTRAN  Subroutine  UI 


Function  UF 


Get  day  of  the  week 

Convert  aggregate  avionics  MTBD  from  flying  hours  tc 
real  time 

-  Depends  on  day  of  the  week 

-  Depends  on  time  of  day  (day  or  night) 

Simulate  uniform  (0,1) 

-  Convert  to  exponential  time  between  LRU  failures 
Simulate  uniform  (0,1) 

-  Compare  uniform  number  to  LRU  failure  rate 
percentages  to  determine  which  LRU  has  failed 

Get  data  on  failed  LRU  and  place  attribute  values  ir 
transact i on 

-  Attribute  for  test  station  requirement 

-  Attribute  for  technician  (AFSC)  requirement 

-  Attribute  for  LRU  maintenance  (to  distinguish  frc 
ATE  maintenance) 

Simulate  uniform  (0,1) 

-  Compare  uniform  number  to  LRU  RTS,  RTOK  and  NRTS 
percentages  to  determine  type  of  maintenance  acti 

If  maintenance  action  is  RTS  (repair): 

Assign  maintenance  action  attribute  for  RTS 
Compute  mean  repair  time 

-  Simulate  two  uniform  (0,1) 

-  Convert  to  Erlang  repair  time  with  k=2 

-  Divide  repair  time  into  two  attributes 

—  Attribute  for  setup  and  fault  isolation  time 
—  Attribute  for  time  to  remove  and  replace  SRU 
teardown  LRU 

-  Simulate  uniform  (0,1) 

—  Compare  uniform  number  to  SRU  spare  part 
availability  percentage  to  determine  SRU 
availability  for  repair 
—  Set  attribute  for  SRU  spare  availability 


Figure  3.  FORTRAN  User  Function  UF 


If  maintenance  action  is  RTOK: 

-  Assign  maintenance  action  attribute  for  RTOK 

-  Compute  mean  task  time  for  RTOK 

-  Simulate  five  uniform  (0,1) 

-  Convert  to  Erlang  task  time  with  k=5 
Divide  task  time  into  two  attributes 
—  Attribute  for  setup  and  test  time 

—  Attribute  for  second  test  (if  required)  and 
teardown  time 

-  Set  attribute  for  SRU  spare  part  availability  (SRU 
spare  not  needed  for  RTOK) 

If  maintenance  action  is  NRTS: 

-  Assign  maintenance  action  attribute  for  NRTS 

-  Compute  mean  task  time  for  NRTS 

-  Simulate  five  uniform  (0,1) 

-  Convert  to  Erlang  task  time  with  k=5 

-  Divide  task  time  into  two  attributes 

—  Attribute  for  setup  and  fault  isolation  time 
—  Attribute  for  teardown  time 

-  Set  attribute  for  SRU  spare  part  availability  (SRU 
spare  not  needed  for  NRTS) 


Figure  3-  FORTRAN  User  function  UF  (Continued) 


technicians  actually  remove  the  failed  LRU  from  the 
aircraft;  the  removal  time  is  simulated  as  an  Erlang 
distribution  with  a  mean  of  two  hours  and  a  k-value  of  four. 
The  removal  time  was  assumed  to  be  Erlang.  The  mean  and 
k-value  were  estimated  from  B-52  OAS  data.  Finally,  the 
maintenance  technicians  take  the  failed  LRU  to  the  avionics 
shop.  This  is  assumed  to  take  a  constant  time  of  30 
minutes.  The  technicians  are  then  free  to  perform 
organization  level  maintenance  on  other  LRUs. 

Not  only  do  the  avionics  LRUs  fail  and  require 
maintenance,  but  so  do  the  ATE  test  stations.  This  is 
accomplished  in  the  model  by  a  FORTRAN  user  subroutine  which 
determines  the  ATE  failure  interarrival  times  and  assigns 
attributes  to  each  ATE  failure.  The  user  subroutine  is  very 
similar  to  the  user  function.  A  detailed  task  listing  of 
the  FORTRAN  user  subroutine  is  shown  in  Figure  4. 

Once  the  failed  LRUs  (or  in  an  abstract  sense,  the 
failed  ATE  test  stations)  arrive  at  the  shop,  they  must  wait 
in  queue  until  resources  are  available  so  that  the 
intermediate  level  (avionics  shop)  maintenance  can  be 
performed.  First,  the  failed  LRU  (or  station)  must  wait 
until  a  station  of  the  right  type  is  available.  The  model 
allows  for  four  station  types  with  a  user  input  quantity  of 
stations  of  each  type.  The  four  types  of  stations  for  the 
B-1B  are  the  Digital  station,  the  Digital  Analog  Video 
station,  the  Radio  Frequency  station,  and  the  Radar 


Subroutine  US 


Simulate  uni-form  (0,1) 

-  Convert  to  exponential  time  between  ATE  failures 
Simulate  uniform  (0,1) 

-  Compare  uniform  number  to  ATE  failure  rate  percentage 
to  determine  which  station  type  has  failed 

Place  attribute  values  for  station  repair  transaction 

-  Attribute  for  station  type 

-  Attribute  for  technician  (PMEL)  requirement 

-  Attribute  for  ATE  maintenance  (to  distinguish  from 
LRU  maintenance) 

-  Attribute  for  type  of  maintenance  action  (station 
repair  is  only  RTS) 

Simulate  two  uniform  (0,1) 

-  Convert  to  Erlang  repair  time  with  k=2 

-  Divide  repair  times  into  two  attributes 

—  Attribute  for  setup  and  fault  isolation  time 
—  Attribute  for  time  to  remove  and  replace  SRU 

Simulate  one  uniform  (0,1) 

-  Compare  uniform  number  to  SRU  spare  part  avi lability 
percentage  to  determine  SRU  availability  for  repair 

-  Set  attribute  for  SRU  spare  part  availability 


Figure  4.  FORTRAN  User  Subroutine  US 


Electronic  Warfare  station.  The  FORTRAN  user  function 
determines  by  a  random  number  the  actual  LRU  that  failed. 
Each  LRU  has  a  known  station  type  requirement  and  that 
station  type  requirement  is  set  as  an  attribute  value  by  the 
user  function.  The  LRU  is  then  routed  to  the  correct 
station  queue  by  conditional  branching.  Similar  branching 
occurs  for  station  maintenance.  Second,  the  LRU  (or 
station)  must  also  wait  until  a  technician  of  the  right  type 
is  available.  The  model  allows  for  the  five  AMS  maintenance 
technician  types  described  in  Section  II.  This  includes  the 
four  types  for  LRU  maintenance,  and  also  the  PMEL  technician 
type  for  ATE  station  maintenance.  The  quantity  of 
maintenance  technicians  (direct  labor  per  shift)  was  shown 
in  Table  II  in  Section  III.  The  FORTRAN  user  function 
determines  by  a  random  number  the  actual  LRU  that  failed. 
Each  LRU  has  a  known  technician  type  requirement  and  that 
technician  type  requirement  is  set  as  an  attribute  value  by 
the  user  function.  The  LRU  is  then  routed  to  the  correct 
technician  queue  by  conditional  branching.  All  station 
maintenance  is  routed  to  the  PMEL  technician  queue.  For 
example,  if  the  failed  LRU  is  from  the  B-1B  ECM  system,  then 
the  failed  LRU  will  require  (say)  a  Radar /EW  test  station 
and  a  defensive  avionics  technician. 

Once  both  (station  and  technician)  resources  are 
available,  the  technician  performs  the  first  part  of  the 
maintenance  action.  For  ATE  maintenance  and  for  LRU  repair, 


this  includes  the  time  for  setup  and  fault  isolation.  For 
LRU  RTOK  actions,  this  includes  the  time  to  set  up  and  run 
one  performance  test.  For  LRU  NRTS  actions,  this  includes 
the  setup  and  fault  isolation  time.  For  all  cases,  the  time 
to  accomplish  the  first  part  of  the  maintenance  action  is 
simulated  by  the  FORTRAN  user  function  (or  subroutine)  and 
stored  as  an  attribute. 

For  LRU  and  station  repairs,  a  SRU  spare  is  required 
before  any  further  maintenance  may  continue  on  the  LRU  or 
station.  In  the  FORTRAN  user  function  or  subroutine,  there 
is  a  907.  chance  that  the  correct  replacement  SRU  spare  will 
be  available  at  base  level.  If  the  spare  part  is  not 
available  for  an  LRU  repair,  the  LRU  is  taken  off  of  the 
test  station  and  placed  in  awaiting  parts  (AWAP)  status. 

The  station  and  technician  resources  are  then  freed  to  do 
other  work.  If  the  spare  part  is  not  available  for  station 
maintenance,  then  only  the  technician  resource  is  freed  to 
do  other  maintenance.  The  station  stays  in  down-fo» — parts 
status  until  the  correct  spare  part  can  be  ordered  and 
shipped.  For  both  LRUs  and  stations,  the  order  and  ship 
time  for  SRU  spares  is  assumed  to  be  a  constant  8  days.  On 
the  other  hand,  if  the  correct  spare  part  is  available  at 
base  level,  then  the  technician  obtains  the  SRU  spare  and 
proceeds  with  the  maintenance.  The  time  to  obtain  the  SRU 
spare  is  determined  by  probabilistic  branching.  There  is  a 
2 5%  chance  that  the  SRU  spare  will  be  available  at  the 
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forward  supply  point.  This  is  assumed  to  take  a  constant 
time  of  15  minutes.  There  is  a  75/1  chance  that  the  SRU 
spare  will  be  obtained  from  base  supply.  This  is  assumed  to 
take  a  constant  time  of  one  hour.  Of  course,  for  LRU  NRTS 
and  RTOK  actions,  no  SRU  spare  is  required,  and  there  is  a 
separate  branch  in  the  model  for  these  cases. 

The  maintenance  technician  then  peforms  the  second  part 
of  the  maintenance  action.  For  station  and  LRU  repair,  this 
includes  the  time  to  remove  and  replace  the  failed  SRU  and 
the  time  for  teardown.  For  LRU  RTOK  actions,  this  includes 
the  time  for  a  second  performance  test  and  the  time  for 
teardown.  There  is  a  30%  chance  that  a  second  performance 
test  will  be  required.  For  NRTS  actions,  this  includes  the 
time  to  teardown  the  LRU  only  since  no  repair  is  possible. 
For  all  cases,  the  time  to  accomplish  the  second  part  of  the 
maintenance  action  is  simulated  by  the  FORTRAN  user  function 
(or  subroutine)  and  stored  as  an  attribute.  Once  the  second 
part  of  the  maintenance  is  complete,  the  station  and 
technician  are  then  freed  to  do  other  work.  After 
maintenance  is  complete,  LRUs  are  then  taken  to  base  supply. 

Numerous  statistic  nodes  have  been  included  in  the 
model  to  measure  the  performance  of  the  maintenance  queueing 
system.  Interval  statistics  measure  the  base  repair  cycle 
time  which  is  the  time  from  LRU  failure  to  return  to  base 
supply.  Statistics  are  collected  separately  for  each 
station  type  (for  example,  the  Digital  station).  LRU 


maintenance  statistics  are  also  kept  separate  -from  ATE 
maintenance  statistics.  The  LRU  NRTS  actions  are  also 
removed  -from  the  statistics  since  no  repair  was  actually 
accomplished;  the  model  figure  of  merit  is  base  repair  cycle 
time. 

Finally,  the  model  keeps  track  of  all  changes 
associated  with  the  shift  changes.  Maintenance  technicians 
and  stations  are  available  from  0800  until  2400  each  working 
day  (i.e.,  two  maintenance  shifts  per  day).  Resource  alter 
nodes  are  used  to  control  the  beginning  and  end  of  each 
shift.  Since  Q-GERT  alter  nodes  are  nonpreempti ve,  however, 
this  means  that  a  maintenance  task  must  be  completed  before 
the  technician  is  allowed  to  go  home.  This  is  why  the 
maintenance  time  was  divided  into  two  parts.  It  is  possible 
to  have  a  technician  complete  the  first  part  of  a 
maintenance  action,  go  home  at  the  end  of  a  shift,  and 
complete  the  second  part  of  the  maintenance  action  the 
following  day.  Dividing  the  maintenance  time  into  two  parts 
therefore  minimizes  the  actual  amount  of  "overtime"  which 
occurs  in  the  simulation  model.  Actual  experience  with  the 
model  indicates  that  the  amount  of  "overtime"  varies  from 
zero  to  two  hours  per  day  (both  shifts  combined)  depending 
on  the  station  workload.  This  is  considered  reasonable  and 
realistic  since  the  SAC  maintenance  policy  does  allow  for 
occasional  "graveyard"  maintenance  to  respond  to  maintenance 
workload  requirements.  The  resource  alter  nodes  also  model 


scheduled  maintenance  on  the  ATE  test  stations.  The 
stations  are  assumed  to  require  a  daily  confidence  test  at 
the  beginning  of  the  first  maintenance  shift.  This  is 
assumed  to  take  a  constant  time  of  15  minutes  for  each 
station. 


B®§l9Q  °f  Simulation  Experiment 

The  purpose  of  the  simulation  experiment  is  to 
determine  if  varying  test  station  quantities  can  influence 
LRU  base  repair  cycle  time.  Therefore,  base  repair  cycle 
time  is  the  dependent  variable  of  interest,  and  test  station 
quantity  is  the  factor  to  be  varied.  Treatments  considered 
will  be  quantities  of  one,  two,  or  three  test  stations  (of 
each  station  type)  at  each  base.  To  reduce  variance  in  the 
simulation  experiment,  the  method  of  common  random  numbers 
is  used.  This  means  that,  for  each  block,  the  base  repair 
cycle  time  is  measured  against  the  three  treatments  with  the 
same  number  of  LRU  failures,  the  same  type  of  maintenance 
actions,  and  the  same  LRU  repair  times.  This,  of  course, 
requires  blocking  to  be  used  in  the  experimental  design 
since  the  observations  within  a  block  are  now  related.  In 
this  experimental  design,  there  are  three  treatments  and  ten 
blocks.  ThiB  experiment  has  18  degrees  of  freedom  in  the 
error  term.  A  crude  rule  of  thumb  is  that  the  degrees  of 
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■freedom  -for  the  error  term  should  be  at  least  10  (22:164). 
This  experiment  was  conducted  a  total  of  12  times  (one  -for 
each  of  -four  station  types  and  one  for  each  of  three 
aircraft  quantities).  A  sample  of  the  simulation  output 
(for  the  16  aircraft  base)  is  shown  in  Table  VI. 

For  each  of  the  12  experiments,  the  factor  of  test 
station  quantity  was  found  to  have  a  statistically 
significant  influence  (with«»<=  0.5)  on  LRU  base  repair  cycle 
time.  Since  the  variance  for  different  treatments  was 
clearly  not  constant  (see  Table  VI),  the  (normal  based) 

ANOVA  was  not  appropriate  to  analyze  the  experimental 
results.  For  this  reason,  the  experimental  results  were 
analyzed  using  Friedman’s  test,  which  is  the  nonpar ametr i c 
equivalent  of  one  way  ANOVA  with  blocking  (complete 
randomized  block  design).  Multiple  pairwise  comparisons  in 
each  of  the  1 2  experiments  were  also  almost  always 
statistically  significant.  In  fact,  32  out  of  36  pairwise 
comparisons  were  found  to  be  statistically  significant.  The 
detailed  data  on  the  experimental  design  is  included  in 
Appendix  B. 

Statistical  si gni f i cance,  however,  may  not  be  the 
critical  issue  in  determining  test  station  quantities. 

These  quantities  should  be  selected  to  minimize  the  overall 
cost  of  avionics  support.  This  cost  should  include  both  the 
cost  of  service  (cost  of  test  stations)  and  the  cost  of 
waiting  (cost  of  avionics  LRU  spares).  Test  station 
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quantities  should  be  selected  to  achieve  the  optimum  balance 
between  cost  of  service  and  cost  of  waiting  and  achieve  the 
lowest  overall  cost.  This  tradeoff  will  be  analyzed  in 
Section  VI. 

Test  for  Transient  Conditions 

Any  simulation  experiment  must  be  concerned  with  the 
effect  of  transient  conditions.  In  the  case  of  a  queueing 
model,  if  the  queues  are  initially  empty,  it  may  take  a 
certain  period  of  time  for  the  queue  length  to  grow  to  its 
steady  state  size.  The  simulation  experiment  described 
earlier  used  simulation  periods  of  6240  hours,  corresponding 
to  six  months  of  avionics  maintenance.  To  test  if  this 
period  was  sufficiently  long  to  achieve  steady  state 
conditions,  the  experiment  was  run  again  for  a  simulation 
period  of  12400  hours  for  all  four  station  types  with  five 
replications  each.  The  Wilcoxon  Signed  Ranks  Test  was  used 
to  test  if  there  was  any  statistically  significant 
difference  in  the  average  base  repair  cycle  time  between  the 
two  simulation  periods.  A  nonparametric  test  was  again 
selected  since  variance  was  not  constant  for  different  test 
station  utilizations.  For  any  10,  there  is  no 

statistically  significant  difference  in  b^se  repair  cycle 
time  between  the  two  simulation  periods.  The  conclusion 


reached  is  that  te  simulation  period  of  6240  hours  is 
sufficent  to  achieve  steady  state  results. 

Sensi.ti.yi.ty  Analyses 

Considerable  sensitivity  analysis  was  conducted  on  the 
various  model  inputs.  This  is  essential  since  the  baseline 
analysis  was  based  on  preliminary  contractor  estimates. 

Major  elements  investigated  include  LRU  reliability,  LRU 
test  times,  ATE  reliability,  ATE  repair  times,  and  other 
■factors  currently  used  in  the  model.  It  appears  that  the 
test  station  quantities  may  vary  considerably  with  only 
modest  and  quite  credible  changes  to  most  o-f  the  major 
elements.  For  this  reason,  it  is  too  early  in  the  B-1B 
program  to  precisely  determine  ATE  test  station  quantities. 
Rather,  it  is  only  possible  to  determine  a  reasonable  range 
of  quantities.  The  range  of  station  quantities  is  presented 
in  Section  VI. 

Verification  and  Validation 

Verification  means  ensuring  that  the  model  behaves 
exactly  as  it  is  intended.  To  assist  in  the  verification 
process,  the  simulation  model  is  the  synthesis  of  5  smaller 


models  which  were  developed  earlier.  Each  o-f  these  smaller 
models  corresponds  to  a  major  portion  of  the  final 
simulation  model.  Spec i f i cal  1 y ,  a  smaller  model  was 
developed  for  (1)  B-1B  flying  hour  profile  and  LRU  failure 
generation,  (2)  resource  allocation  (both  test  stations  and 
technicians),  (3)  ATE  failure  generation,  (4)  detailed  test 
and  repair  procedures,  and  (5)  FORTRAN  user  functions  and 
subroutines  used  in  the  model.  Each  of  these  smaller  models 
was  programmed  with  very  detailed  output  and  each  was 
sufficiently  simple  to  allow  manual  verification. 

Validation  is  a  process  to  ensure  that  the  model 
real i sti cal  1 y  portrays  the  real  world.  The  most  important 
element  of  the  continual  process  of  validation  was  the 
coordination  of  all  major  ground  rules  and  assumptions  with 
personnel  from  HQ  SAC  that  have  had  actual  experience  with 
avionics  maintenance  in  the  B-52  and  FB-111  programs.  These 
individuals  provided  a  significant  amount  of  feedback  and 
constructive  criticism.  In  addition,  an  effort  is  now  being 
conducted  to  run  this  model  on  the  B-52  DAS  and  its  ATE. 
Model  output  will  then  be  compared  to  the  actual  operational 
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One  of  the  most  important  theorems  at  queueing  net 
theory  is  known  as  Jackson’s  result;  the  theorem  applie 
a  so-called  open  network  ot  queues.  An  open  network  is 
multiple  queue  system  in  which  a  customer  typically  req 
service  at  more  thar  one  queueing  station.  It  is  assum> 
that  customers  enter  the  network  at  one  ot  the  queueing 
stations,  queue  tor  service,  and  then  depart  the  networ 
else  proceed  to  another  queueing  station  in  the  network 
additional  service.  It  is  also  assumed  that  each  queue 
station  has  s  parallel  servers  with  exponential  service 
times.  It  is  also  assumed  that  the  time  between  arriva 
(trom  customers  outside  the  network)  at  each  queueing 
station  is  exponential.  Finally,  it.  is  assumed  that 
customers  at  a  given  queueing  station  either  depart  the 
network  or  move  to  another  queueing  station  according  t< 
specified  probability  distribution  which  is  not  state 
dependent.  Jackson’s  result  states  that  for  such  a  neti 
the  state  pr obab i 1 i t i es  and  performance  measures  for  ea< 
queueing  station  may  be  obtained  by  applying  the  simple 
exponential  (M/M/s)  model  to  each  queueing  station 
spp ar  at  el y  ( 5 :  146-161) . 
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A  second  powerful  technique  is  mown  as  the  methods  of 
phases.  This  technique  applies  to  service  times  (or 
interrarival  times)  which  are  Erlang  or  hyperexponential . 

In  many  cases,  it  is  still  possible  to  obtain  a  closed-form 
expression  for  the  state  probabilities  by  transform  methods 
The  solution  technique  is  to  redefine  the  state  space;  the 
state  of  the  system  does  not;  refer  to  the  number  of 
customers  in  the  system  but  rather  to  the  number  of 
exponential  phases  yet  to  be  completed  in  the  system.  The 
key  to  the  solution  is  to  express  the  Erlang  distribution  a 
the  sum  of  exponential  distributions,  or  to  express  the 
hyperexponential  distribution  as  the  weighted  average  of 
exponential  distributions.  In  either  case,  it  is  then 
possible  to  construct  classical  exponential  rate  diagrams 
for  the  queueing  system.  A  closed-form  expression  may  then 
often  be  obtained  by  use  of  transforms  (4:171-175). 

A  third  technique  of  queueing  theory  involves  multiple 
classes  of  customers  with  different  mean  service  times.  It 
is  still  required  that  the  customers  arrive  and  are  served 
according  to  an  exponential  distribution.  The  state  space 
for  this  type  of  problem  is  a  set  of  vectors;  each  entry  in 
the  vector  represents  the  number  of  customers  of  each  class 
Depending  on  the  situation,  there  may  not  be  a  closed-form 
expression  for  the  state  probabilities.  It  is  possible, 
however,  to  write  the  balance  equations  for  the  state  space 
as  a  linear  system  of  equations  and  unknowns.  The  unknowns 


are  the  state  probabilities  -for  the  system.  The  solution 
technique  is  to  approximate  the  state  space  by  a  -finite 
number  o-f  elements.  For  most  queueing  situations,  the 
steady— state  number  o-f  customers  in  the  system  -falls  within 
a  finite  range  with  a  probability  close  to  one.  For 
example,  it  might  be  a  reasonable  approx  i  mat.  i  on  to  postulate 
that  the  steady-state  number  of  customers  in  the  system 
never  exceeds  a  certain  number.  For  the  finite  case,  there 
is  always  one  dependent  equation  in  the  set  of  balance 
equations.  The  solution  technique  is  to  remove  one  of  the 
equations  and  then  add  an  equation  which  states  that  the 
probab i 1 i t i es  over  the  state  space  sum  to  one.  The  final 
step  is  to  solve  the  finite  system  of  equations  and  unknowns 
by  a  numerical  technique  such  as  Gauss-Seidel  iteration 
(4:123,158-165).  By  using  a  simple  formula  for  expected 
value,  it  is  then  possible  to  find  L,  the  expected  number  in 
the  system,  for  each  class.  Little’s  formula  may  then  be 
used  to  find  the  expected  waiting  time,  W,  for  each  customer 
class  ( ■  j* "  26 1  262 )  . 

The  actual  technique  developed  in  this  research  effort 
is  a  synthesis  of  the  three  queueing  techniques  discussed 
earlier.  The  LRUs  that  require  maintenance  compose  the 
first  class  of  customers;  the  ATE  stations  that  require 
maintenance  compose  the  second  class  of  customers.  The 
repair  time  at  intermediate  level  is  treated  as  Erlang  and 
the  method  of  phases  is  used.  The  transition  from 
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is  treated  as  a  network  of  queues. 


WQ^el.  Qescr i_gjt i_gn 

The  analytical  queueing  model  treats  B-1B  avionics 
maintenance  as  a  queueing  network;  this  network  is  shown  in 
Figure  5.  The  avionics  LRUs  are  assumed  to  arrive  to  the 
organi zati onal  maintenance  squadron  queue  according  to  a 
(homogeneous)  Poisson  process.  The  servers  -for  this  queue 
represent  pairs  o-f  OMS  maintenance  technicians.  Since  there 
are  38  OMS  maintenance  technicians  per  shift,  there  are  19 
pairs  o-f  technicians  and  thus  19  servers.  The  LRU  removal 
time  is  assumed  to  be  exponential  with  a  mean  o-f  3.5  hours. 
After  passing  through  the  OMS  queue,  the  LRUs  then  pass  on 
to  the  test  station  queues  in  one  of  two  ways.  One  way  is 
for  the  LRUs  to  go  into  awaiting  parts  (AWAP)  status,  which 
is  treated  as  an  infinite  server  queue  where  the  service 
time  represents  the  order  and  ship  time  for  the  SRU  spare. 
The  purpose  of  the  branch  is  to  make  sure  that  the  waiting 
time  for  SRIJ  spares  is  included  in  the  base  repair  cycle 
time.  The  other  way  is  for  the  LRUs  to  proceed  directly  to 
the  test  stations.  This  corresponds  to  the  case,  which 
happens  for  907.  of  the  repairs  and  for  all  RTOKs  and  NRTS 
actions,  where  it  is  not  necessary  to  order  a  SRU  spare 
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from  depot  level.  In  either  case,  LRUs  arrive  to  each  of 
the  test  stations  according  to  a  Poisson  process.  In 
addition,  ATE  failures  arrive  to  the  test  station  queues, 
also  according  to  a  Poisson  process.  The  LRU  failures  and 
the  ATE  failures  are  treated  as  two  classes  of  customers. 

The  ATE  failures  require  an  exponential  service  time  with  a 
mean  of  one  hour,  taut  the  LRU  failures  require  an  Erlang 
(with  k=2)  service  time  with  a  mean  which  is  computed  by  the 
model .  The  model  assumes  that  al 1  LRUs  have  the  same  mean 
repair  time;  this,  of  course,  is  only  a  crude  approximation. 
In  fact,  each  LRU  may  have  a  different  mean  repair  (service) 
time.  The  model  calculates  the  LRU  mean  repair  time  as  a 
(reliability)  weighted  average  of  all  mean  LRU  repair  times. 
Each  of  the  mean  LRU  repair  times,  in  turn,  is  a  weighted 
average  of  the  RTS,  RTOK  and  NRTS  mean  repair  times.  The 
detailed  task  flow  for  the  analytical  queueing  model  is 
shown  in  Figure  6. 


Solution  lechnigues 


The  organizational  maintenance  and  awaiting  parts 
queues  are  both  simple  exponential  (M/M/s)  queues.  Simple 
closed  form  expressions  for  the  state  probabilities  and  the 
waiting  time  for  the  M/M/s  queue  are  found  in  Ross 
(3:264-267);  the  model  has  a  FORTRAN  subroutine  which  uses 
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Get  data  -from  subroutine  BLOCK  DATA  through  FORTRAN 
COMMON  statement 

Add  LRU  -failure  rates  to  obtain  aggregate  Mean  Time 
Between  Demand  ( MTBD )  for  LRUs  assigned  to  each  station 
type 

Get  overall  aggregate  MTBD 
Solve  QMS  queue  as  M/M/ 19  queue 

Compute  mean  LRU  repair  time  by  station  type 

-  Weighted  average  by  RTS,  NRTS,  and  RTQK  percentages 

-  Weighted  average  by  LRU  reliability 

Solve  order  and  ship  queue  for  LRUs  awaiting  parts 

Solve  test  station  queues 

-  Solve  case  for  single  server 
Solve  case  for  two  parallel  servers 

Compute  total  network  waiting  time  for  LRUs  by  station 
type 


Figure  6.  Analytical  Queueing  Model 


these  expressions  to  calculate  the  li/M/s  state  probabilities 


and  the  queue  waiting  time.  It  is  also  true  that  the  output 
of  a  simple  exponential  (M/M/s)  queue  is  a  Poisson  process 
with  the  same  arrival  rate  as  the  original  input  Poisson 
process.  It  is  also  true  that  when  this  output  Poisson 
process  is  channeled  to  different  queues  with  fixed 
probabi 1 i ti es,  the  arrival  process  at  each  of  the  subsequent 
queues  is  a  Poisson  process  with  a  rate  equal  to  the 
original  overal 1  arrival  rate  multiplied  by  the  probability 
of  the  arrival  process  being  channeled  to  that  particular 
queue.  This  means  that  the  LRU  arrival  process  at  each  of 
the  test  stations  is  a  Poisson  process. 

The  solution  technique  for  the  test  station  queues  is 
somewhat  more  involved.  The  rate  diagram  for  the  single 
server  case  is  shown  in  Figure  7.  The  state  of  the  system 
is  now  a  vector.  The  I  TU  repair  process  is  divided  into  two 
phases.  Since  the  repair  time  has  an  Erlang  distribution 
(with  k=2)  ,  the  time  to  complete  each  phase  has  an 
exponential  distribution.  The  number  in  the  first  entry  of 
the  state  vector  represents  the  number  of  LRUs  yet  to  be 
served;  the  number  in  the  second  entry  represents  the  number 
of  LRUs  which  have  completed  the  first  phase  of  the  repair 
but  have  yet  to  complete  the  second  phase.  It  is  assumed 
that  work  on  an  LRU  which  has  completed  the  first  phase 
takes  priority  over  work  on  any  more  recent  LRU  arrivals. 
This  means,  for  the  one  server  case,  that  the  number 
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in  the  second  entry  at  the  state  vector  never  exceeds  one. 
The  third  entry  on  the  state  vector  is  0  when  there  are  no 
ATE  -failures  in  the  queue  (the  test  station  is  up)  and  1 
when  there  is  an  ATE  failure  in  the  queue  (the  station  is 
down).  ATE  test  station  repair  always  takes  priority  over 
any  LRU  maintenance.  In  addition,  an  ATE  test  station 
failure  preempts  any  on-going  LRU  maintenance.  For  such  an 
event,  because  of  the  memoryless  property  of  the  exponential 
distribution,  the  preemption  is  equivalent  to  having  the  LRU 
repair  phase  start  over. 

No  closed  form  expression  for  system  performance 
measures  may  be  found  for  the  queueing  system  in  Figure  7. 
This  is  due  to  the  preemptive  blocking  caused  by  ATE 
failures.  The  solution  technique  is  a  four  step  process. 
First,  the  state  space  is  approximated  by  a  finite  number  of 
terms.  In  this  case,  it  is  assumed  that  the  number  of  LRUs 
in  the  queue  never  exceeds  8.  Second,  the  rate  diagram  and 
the  associated  balance  equations  can  be  used  to  construct  a 
linear  system  of  equations  and  unknowns.  The  balance 
equations  are  based  on  the  "rate  in"  =  "rate  out"  principle. 
For  example,  the  balance  equation  for  the  vector  (3,1,1) 
would  be: 


In  general  ,  the  unknowns  are  the  probab 1 1 i t i es  of  each 


element  in  the  state  space.  The  coefficients  of  the 
unknowns  are  various  combinations  of  the  arrival  rates  and 
the  service  rates  of  the  LRU  failures  and  the  ATE  station 
failures  (denoted  by  and  l  * *  It  is  also 

necessary  to  remove  one  balance  equation  and  in  its  place 
add  an  equation  which  states  that  the  sum  of  all  the 
probabi 1 i t i es  is  one.  Third,  the  linear  system  of  equations 
and  unknowns  is  solved  numerically.  The  queueing  model 
developed  in  this  research  effort  uses  IMSL  routine  LEQT1F 
to  solve  the  system  of  equations  and  unknowns.  This  routine 
performs  standard  Gaussian  elimination  with  equi 1 i brat i on 
and  partial  pivoting  ( 19: LEQT 1 F- 1 , —3 ) .  Fourth,  the  state 
probabilities  are  then  used  to  find  ,  the  expected 
steady-state  number  of  LRUs  in  the  system.  This  is  a  simple 
calculation  based  on  the  idea  of  conditional  expectation. 
Little’s  Formula  can  then  be  used  to  calculate  the  total 
waiting  time  in  the  system  for  the  LRUs.  This  equation 
would  be: 

U 1L  -  Ll  / X±  ,7> 

This  computes  the  total  LRU  waiting  time  for  each  test 
stat i on  queue. 

The  model  then  computes  the  average  base  repair  cycle 
time  for  the  entire  LRU  maintenance  process  tor  each  station 
queue.  This  is  the  sum  of  (1)  the  total  time  in  the  system 


for  the  OMS  queue,  (2)  the  order  and  ship  time  multiplied  by 
the  -fraction  of  maintenance  actions  that  require  SRU  spares 
to  be  ordered  and  shipped,  and  (3)  the  total  time  in  system 
for  the  test  station  queue. 

The  case  for  two  servers  (two  test  stations)  is  handled 
very  similarly.  However,  the  rate  diagram  is  somewhat  more 
complicated.  A  portion  of  the  two  server  rate  diagram  (for 
states  (4,0,0)  through  (4,1,2))  is  shown  in  Figure  8.  Note 
that  ATE  maintenance  always  takes  priority  over  any  LRU 
maintenance.  Also,  LRU  repair  of  LRUs  in  the  second  phase 
always  takes  priority  over  LRUs  in  the  first  phase.  The 
solution  technique  for  the  two  server  case  follows  the  same 
steps  as  those  for  the  single  server  case.  The  FORTRAN 
computer  code  listing  for  the  analytical  model  is  included 


Model  Results 


The  analytical  queueing  model  was  run  tor  the  case  ot 


16  aircraft  per  base.  The 

resu its 

are  shown  in  1 ata 3 

Table  VII.  Analytical 

M  o d el  Resu 1 1 s 

Base  Repair 

Cycle  1 

1  me  (Hours) 

1 

Station 

Case  2  Station 

Digital 

2  3 .  hi 

18.9 

D  i  g  i  t  a  1  An  a  1  og  V  i  dec) 

26.  1 

1  9 .  3 

RF 

26 . 3 

20.  0 

R/EW 

26.  7 

20 . 0 

For  reasons  which  will  he  discussed  in  the  next  sub-section, 
the  queueing  model  estimates  tor  base  repair  cycle  time  are 
con  si  stent! v  lower  than  the  estimates  obtained  trom  the 
simulation  model.  Comparable  tiguree  tor  the  16  aircraft 
base  were  derived  from  the  simulation  model  and  are  shown  in 
Table  III.  fac::h  figure  is  the  average  of  10  replications. 


I  able  VIII.  Simulation  Model  Results 
H  a  s  e  R  e  p  a  i  r  ( 1  y  c::  1  e  I  i  m  e  ( H  o  u  r  s  ) 


1  Station  Case 

35.  1 

h  1 . 8 
69 . 2 


8 1  a  t 1  o  n 
25.  0 
24.9 
28.  7 


L  a  s  e 
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Analytical  versus  §lffiyl.§tion  Model. 

In  comparing  the  analytical  model  results  to  the 
simulation  model  results,  it  was  found  that  the  analytic 
estimates  of  base  repair  cycle  time  (total  time  in  syste 
were  consistently  lower  than  the  estimates  obtained  from 
simulation  model.  Moreover,  it  was  found  that  this  gap 
widened  as  the  average  station  utilization  increased.  I 
comparing  the  one  station  case  of  Table  VII  to  Table  VII 
it  can  be  seen  how  the  gap  gets  larger  as  the  average 
station  utilization  gets  larger.  The  major  reasons  for 
gap  follow: 

1 .  Far  most,  queueing  situations,  the  total  time  in  sys' 
depends  not  only  on  the  mean  of  the  service  time  but  alsi 
the  variability  of  the  service  time.  Increasing  this 
variability  causes  an  increase  in  the  total  time  in  systi 
There  is  good  reason  to  believe  that  the  analytical  mode 
underestimates  the  true  variability  associated  with  LRU  . 
ATE  repair.  First.,  for  LEU  repair,  the  analytical  model 
assumes  that  all  LRUs  have  the  same  mean  repair  time  and 
uses  weighted  averages  to  find  that.  mean.  In  fact.,  then 
are  different  means  for  RTS,  RTOK  and  NETS  actions,  and 
different  means  for  different  LRUs .  This  implies  that,  tl 
actua.  I  distribution  of  LEU  repair  time  is  not  a  simple 
Ft  land,  tint  i  at  tier  is  a  very  complicated  hyperexponent  1  a 


(with  a  larger  variability) .  Second,  a  similar  problem 
occurs  with  ATE  test  station  maintenance.  1  he  analytical 
model  allows  -tor  only  one  test  station  -failure  rate  and  one 
mean  test  station  repair  time.  In  fact.,  there  are  two 
different  test  station  service  times  depending  on  whether  or 
not.  a  SRU  spare  is  available  at  base  level.  This  means  that 
the  actual  station  service  time  also  has  a  hyperexponential 
distribution  (with  a  1 arger  variabi I i ty) . 

2.  The  simulation  model  has  a  variable  LRU  failure 
(arrival)  rate.  Specifically,  the  simulation  model  has  a 
different  arrival  rate  for  different  days  of  the  week. 
Moreover,  the  LRUs  arrive  in  two  spikes  during  each  day 
(once  from  the  morning  sorties  and  once  from  the  night 
sorties).  This  means  that  LRUs  are  more  likely  to  arrive 
with  other  LRUs,  causing  longer  waiting  times  in  queue.  On 
the  other  hand,  the  analytical  model  assumes  that  the 
avionics  LRUs  arrive  according  to  a  Poisson  process  which  is 
homogeneous  over  the  16  hour  maintenance  day.  This  has  the 
effect  of  leveling  out.  the  LRU  workload  and  causing  shorter- 
waiting  t i m e s  in  qu e u e . 

3.  The  analytical  model  may  have  a  downward  bias  in  total 
waiting  time  since  it  uses  a  finite  state  space  to 
approximate  the  actual  distribution  of  customers  (LRUs)  in 
the  system.  Although  the  probability  of  ever  having  more 


than  (say)  8  customers  in  the  system  is  quite  small,  the 
total  contribution  of  9  or  more  customers  to  the  true 
expected  number  o-f  customers  in  the  system  may  not  be 
insignificant.  After  all,  these  albeit  small  probabilities 
are  associated  with  1 arge  numbers  of  customers.  The  initial 
analytical  model  assumed  a  maximum  of  8  LRUs  in  the  system. 
The  model  for  the  single-server  case  was  also  run  with 
maxi  mums  of  12  and  24  LRUs  in  the  system.  This  caused  the 
value  of  base  repair  cycle  time  (total  time  in  system)  to 
increase  by  two  or  more  hours  in  each  case,  and  the  rate  of 
growth  did  not  appear  to  decline.  An  analysis  of  the 
simulation  model  output  indicated  that  this  problem  is  quite 
severe.  For  the  16  F'AA,  one  station  case,  the  maximum 
number  of  LRUs  in  each  queue  over  a  six  month  period  was 
found  to  be  4,  24,  58,  and  66  for  the  Digital,  Digital 
Analog  Video,  RF  and  R/EW  stations,  respectively.  The 
conclusion  drawn  from  this  is  that  the  finite-state  space 
approximation  causes  a  significant  downward  bias  in  the 
model  estimate  of  the  total  time  in  system. 


The  analytical  model,  overall,  requires  many  crude  and 
simplifying  assumptions;  much  more  so  than  does  the 
simulation  model.  Empirical  results  suggest  that  these 
differences  in  assumptions  cause  a  significant  difference 
between  the  output  of  the  two  models.  Since  the  assumptions 
of  the  simulation  model  are  much  more  realistic,  it  was 


assumed  that:  the  simulation  model  has  a  more  realistic 
output.  With  this  assumption,  the  simulation  model  was 
selected  over  the  analytical  queueing  model  as  the  best 
choice  of  the  B-1B  System  Program  Office. 

Possible  Further  Research  for  the  Analytical  Model 

The  previous  discussion  identified  three  serious 
problems  with  the  analytical  model  developed  during  this 
research  effort.  Further  research  could  be  accomplished  to 
see  if  these  problems  could  be  overcome.  The  first  problem 
was  that  the  present  analytical  model  underestimates  the 
true  variability  of  LRU  and  ATE  repair.  A  possible  solution 
is  to  model  these  repairs  by  using  hyperexponential 
distributions.  This  might:,  however,  cause  numerical 
problems,  since  the  present  two-server  algorithm  is  already 
at  72  equations  and  72  unknowns.  The  second  problem  was 
that  the  LRU  arrival  rate  is  not  constant,  which  the  present 
analytical  model  ignores.  One  possible  solution  might  be  to 
model  the  LRU  arrivals  as  a  hulk  arrival  process,  for  which 
limited  queueing  results  are  available  (18:199-219).  The 
third  problem  was  caused  by  the  use  of  a  finite  state 
approx i mat i on .  It  may  not  be  computationally  feasible  to 
simply  expand  the  size  of  the  finite  state  space.  However, 
one  practical  solution  might  be  to  assume  that  the 


probab i 1 i t i es  in  the  finite  space  do  not  sum  to  1.0  but 


rather  to  some  number  close  to  1.0  (say  0.951.  After 
obtaining  a  numerical  solution  to  the  probabilities  of  the 
finite  state  space,  it.  might  then  be  possible  to  devise  some 
kind  of  extrapolation  scheme  to  estimate  the  remaining 
probab i 1 i tes. 


VI  . 


Tradeoff  Studies 


D®t-.§CfQiDiQ(-3  Station  Quantity  Reg u  1  r e m e n t s 

A  statistically  significant  improvement  in  base  repair  cycle 
time  may  have  little  advantage  from  a  logistics  support  cost 
or  weapon  system  availability  point  of  view  as  discussed  in 
Section  IV.  The  purpose  of  this  section,  therefore,  is  to 
develop  criteria  and  techniques  for  determining  the  best- 
test  station  quantities  by  station  type  for  each  MOB. 

The  first  point  to  be  considered  is  that  the  test 
stations  must  be  able  to  accommodate  the  mean  avionics 
workload.  In  other  words,  at  a  minimum,  the  test,  station 
utilizations  must  always  be  smaller  than  1007..  If  this  were 
not  the  case,  the  queues  would  grow  indefinitely,  and  the 
LRU  base  repair  cycle  time  would  became  infinite.  The 
minimum  number  of  stations  which  achieve  station  utilization 
under  1007.  could  therefore  be  regarded  as  an  absolute  floor 
for  the  test  station  quantities.  In  some  cases,  however, 
this  approach  may  not  be  sufficient.  It  is  possible  that  in 
certain  situations  the  base  repair  cycle  time,  although 
finite,  may  nevertheless  be  "excessive"  is  some  sense.  If  a 
rule  or  technique  could  be  developed  which  could  indicate 
when  a  predicted  base  repair  cycle  time  was  "excessive",  it 
would  then  be  appropriate  to  select,  the  minimum  number  of 
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test;  stations  such  that,  no  base  repair  cycle  time  was 
" eMcessi ve " . 

One  approach  would  be  to  compare  the  predicted  base 
repair  cycle  time  to  some  established  standard.  Such  a 
standard  should  not  be  an  arbitrary  number,  but  should  be 
selected  to  achieve  a  necessary  level  of  support.  Such 
support  is  in  the  form  of  the  LRU  base  repair  pipeline.  For 
example,  suppose  that  at  a  given  bast;,  two  avionics 
processors  fail  per  day.  In  addition,  suppose  that  the 
planned  or  desired  base  repair  cycle  time  is  five  days. 

This  means  that  the  base  would  require  an  LRU  pipeline  ot 
ten  spare  processors  at  base  level.  Of  course,  the  actual 
spares  level  computation  would  be  more  complex  than  this. 
First,  it  would  have  to  account  for  the  small  percentage  of 
time  that  the  avionics  LRUs  are  NRTSed  (sent  to  depot)  for 
repair.  Second,  it  would  have  to  make  some  distributional 
assumptions  about,  the  variability  of  the  LRU  arrivals  per 
day  and  about  the  variability  of  the  base  repair  cycle  time. 
The  spares  level  computation  would  then  add  a  safety  stock: 
level  to  the  average  (expected)  pipeline  quantity.  However, 
the  principle  remains  the  same.  The  approach  is  to 
postulate  that  LRU  spares  will  be  procured  based  on  the 
assumption  of  a  planned  base  repair  cycle  time.  This 
planned  base  repair  cycle  time  could  then  become  the? 
standard  to  judge  whether  the  test  station  quantities  were 


sufficiently  large. 


Early  in  a  program  when  little  or  no 


operational  data  is  available,  LRU  spares  are  typically 


procured  based  on  a  base  repair  cycle  time  of  -four  calendar 


days  (20).  It  would  then  be?  a  simple  matter  to  compare  the 


simulation  model  output  to  the  standard  at  -four  days  to 


judge  i f  a  given  base  repair  cv_le  time  were  "excessive". 


The  model  output  is  based  on  -five  working  days  per  week  and 


excludes  weekends.  Thus,  the  four  calendar  day  standard 


must  be  converted  to  2.857  working  days  (by  multiplying  by 


five-sevenths)  or  68.57  hours. 


Given  that  tour  calendar  days  is  to  be  the  standard, 


there  is  still  some  room  tor  judgement  as  to  how  to  compare 


the  simulation  output  to  the  standard.  One  way  would  be  to 


perform  a  one-sided  statistical  test  of  hypothesis. 


Suppose,  for  example,  that  the  simulation  model  were  very 


close  to  the  standard  of  four  days.  Then,  because  of  the 


statistical  noise  present  in  the  simulation  model  output,  it 


would  be  impossible  to  determine  if  the  actual  (expected) 


base  repair  cycle  time  were  in  fact  below  the  standard.  For 


a  one-sided  statistical  test,  suppose  that  the  null 


hypothesis  to  be  tested  is  that  the  base  repair  cycle  time 


is  less  than  or  equal  to  four  days,  and  that  the  alternative 


hypothesis  is  that  base  repair  cycle  time  is  greater  than 


four  days.  This  approach  implicitly  assumes  that  it  is  much 


worse  to  buy  an  unnecessary  test  station  (make  a  type  1 


error)  then  it  is  to  buy  an  insufficient  quantity  of  test 


stations  (make  a  type  II  error).  Similarly,  the  role  of  the 
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null  and  of  the  alternative  hypotheses  can  be  reversed,  and 


the  judgement  at  the  relative  importance  o-f  the  errors  would 
therefore  change  as  well.  For  cases  where  the  station 
utilization  is  under  1007.,  however,  it  is  not  at  all  obvious 
as  to  which  type  ot  error  is  the  more  serious.  If  the  two 
types  of  errors  were  to  be  weighed  equally,  then  the 
one-sided  test  of  hypothesis  would  not.  be  the  correct 
approach.  Assuming  that  the  simulation  model  output  is 
normally  (or  at  lease  symmetrically)  distributed,  then  the 
way  to  weigh  the  two  types  of  errors  equally  is  to  simply 
compare  the  siample  mean  of  the  simulation  output  to  the 
standard.  This  was  the  approach  taken  for  the  remainder  of 
this  research  effort. 


Baseline  Case  and  Sensitivity  Analysis 

The  simulation  model  was  run  for  the  baseline  case 
(described  in  Section  III)  for  each  of  the  MOEis  and  for  each 
of  the  test  station  types.  The  mean  of  10  replications  was 
compared  to  the  four  calendar  day  base  repair  cycle  time 
standard;  minimum  test  station  quantities  were  selected  so 
that  the  standard  could  be  achieved.  These  results  are 
shown  in  Table  IX.  The  baseline  case  shows  an  operational 
requirement  for  24  test  stations  in  total. 
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Table  IX. 


Baseline  Test  Station  Quantities 


Digital 


MOB 

Digital 

An  Video 

RF 

r/ew 

#1 

-  26 

PAA 

1 

1 

#2 

-  16 

PAA 

1 

1 

rj 

#3 

-  32 

PAA 

1 

1 

2 

*“> 

#4 

-  16 

PAA 

1 

1 

n 

n 

TOTAL  4  4  8  8 


The  above  baseline  computation  is  based  on  preliminary 
contractor  data  for  reliability  and  maintainability 
estimates.  A  second  set  of  reliability  figures  was  obtained 
for  HQ  SAC  personnel  (21)  which  reflects  current  experience 
on  like  items  in  current  inventory  aircraft.  This 
sensitivity  is  labeled  HQ  SAC  MTBD  and  is  shown  in  Table  X. 
The  first  baseline  shown  above  also  used  the  GDE  estimates 
for  LRU  performance  test  times.  Another  sensitivity 
analysis  was  run  using  the  RI  performance  test  times  (as 
discussed  in  Section  III).  The  results  for  this  case  are 
labeled  RI  performance  test  times  and  are  shown  in  Table  XI. 
The  HQ  SAC  MTBD  case  shows  an  operational  requirement  t  •  33 
test  stations  in  total,  and  the  RI  performance  test  time 
case  shows  an  operational  requirement  for  29  test  stations 


I  able  X.  HQ  SAC  MTBD  Estimates 

D  i  g  i  t  a  1 


MOB 

Digital 

An  Video 

RF 

R/EW 

#1 

-  26 

PA  A 

1 

V. 

|  -( 

#2 

-  1 6 

PAA 

1 

3 

V 

3 

#3 

—  3  3 

PAA 

3 

3 

3 

#4 

-  16 

PAA 

:! 

2 

2 

2 

TOTAL.. 
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1 0 

1  0 

T  ab  1  e 

X  I  .  PI 

Perfo  r  m  a  n  c  e  T  e  s  t 

1  i  mes 

D  i  g  i  t  a  1 

MOB 

yiaitai 

An  Video 

HP 

R/EW 

#1 

-  26 

PAA 

i 

2 

2 

#2 

-  16 

PAA 

i 

1 

3 

3 

#3 

—  33 

PAA 

3 

3 

3 

3 

#4 

—  16 

PAA 

i 

i 

2 

2 

TOTAL. 

5 

6 

9 

9 

Other  sensitivity  analysis  was  pert'.rmed  -for  the  Is  Phh 
case.  The  baseline  test  station  quani, t; es  for  this  case 
are  one  digital  station,  one  digital  analog  video  station, 
two  radio  frequency  stations,  and  two  radar  EW  stations. 

The  factors  varied  were  ( 1 >  changing  the  SRU  spare 
availability  from  90%  to  80%,  (2)  reducing  the  Alt 

reliability  by  50%,  (3)  improving  the  Oil-  l  eliabi  lit  y  by 

100%,  (4)  increasing  the  Alt:  repair'  times  is,  1  o<J% ,  and  0:0 

reducing  l  tie  Alt'  repair  times  by  50%.  Only  the  first  two 
a  significant  impact  on  I  ..HU  base  repair  cycle 


factors  had 


LRIJ  Bdrif?  h'epdi  r  Uyc  i.  e  }  i  roc?  in  Hour 
16  f-'AA  per  Rase 


Ben si t i vi t v 


Basel i ne 

SRU  Availability 
at  807. 

Reduce  ATE 
Eel i ab i 1 1 1  y 

Imp) rove  ATE 
Rel i  abi 1 i ty 

Increase  ATE 
Repair  Times 

Reduce  ATE 
Repai r  Times 


I)  i  q  i  t  a  1 

Digital  An  Video  Et 

Qty  •”  1  Qty  =  1  Qty  = 


46.8 


R  /  E  W 


51.7 


42.8 


92.  1 


5 1 . 2 


46.  9 


29.  0 


28 . 8 


47.  1 


26.4 
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Logistics  Support  Cost  Tradeoff 


The  first  approach  in  determining  test  station 
quantities  was  to  buy  a  sufficient  quantity  of  test  station 
such  that  the  predicted  base  repair  cycle  time  was  always 
under  the  standard  of  four  calendar  days.  An  alternative 
approach  is  to  compare  the  costs  of  LRU  spares  (pipeline  an 
safety  stock  spares)  to  the  cost  of  test  stations. 
Specifically,  an  add i 1  tonal  test  station  would  be  procured 
as  long  as  the  savings  in  l  E'lJ  spares  (due  to  the  reduced 
base  rep  a  t  r  rye  l  o  1  i  me)  were  greater  It  tan  '.he  <.  ( <  l  of  t  he 
atitlt  I  i  ana  I  station.  1  r  <  order  1  or  this  approach  to  be  vaJ  id 


tJ  7 


it.  is  necessary  to  assume  that  the  actual  procur  emom: 
spares  will  be  based  on  the  actual  base  repair  cvcle  t i me , 
and  not  some  standard  factor . 

the  cost  of  L.RL.I  spares  was  calculated  by  use  of  t  he  Air 
Force  Acquisition  Logistics  Division  •' AFALD)  Loiji  sties 
Supporl'.  Cost  ( L..  SC )  model  ,  I  he  Mean  lime  Between  Demand 
( Ml  BD )  tor  each  LRU,  the  unit',  cost  of  the  LRU,  and  the  LRU 
base  rep)  air  cycle  time  were  the  critical  inputs,  I  he  M  I  BD 
e  s  t  i  rn  a  t  e  s  w  e  r  e  t  h  e  c  o  n  t  r  a  c  t  o  r  e  s  t  i  m  a  fc  e  s  t  h  a  t  w  e  i  e  d  e  s  c  r  1  b  e  d 
in  Section  ill.  The  unit  costs  of  each  LK'U ,  measured  in 
FY  81  dollars,  were  also  obtained  from  contractor  estimates. 
The  base  repair  cycle  time  was  obtained  from  the  output  of 
the  simulation  mode)  for  each  of  the  ID  situations.  I  tie  Lot 
model  ,  which  calculates  the  I..RU  spares  (pipeline  and  safety 
stoci  combined)  quantity,  was  then  used  to  compute  the  total 
LRU  spares  dollar  investment  required  for  each  station  tvpe 
and  station  quantity  at  each  base.  beta! I s  on  the  model  cf n 
be  found  in  the  I  oqisti.es  Support  Lost  Model  User  '  s  l-landbool 
(20),  Other  LRU  logistics  support  costs,  such  as 
maintenance  manhour  costs  or  inventor  v  manaqemen  l  rusts,  cud 
not  vary  wi  lit  >  r  Li a  n  q  e  s  in  I  os  1  si  a  t  i  cm  i  quant  j  1  i  es  and  ripi  r  ■ 
not  considered  in  I  I  us  imi  .  i  ;  isini  ■  ■  l  he-.  u  e  const  ■ants  1  . 

A  1  E  test  s  t  at  i  m  >■;  wor  e  a  s--5!  iim—  d  i  i  «  r  <  c.  !  1  ..  Li  m  II  oi  i  dol  1  ar  ; 

(  in  BY  Li  I  dii !  |  ii  ■  i  {til  i  ■  i  h  it  I  he  lour  at.  a  l  ion  t  vises,  in 

addi  1  i '  ii  i ,  i  I.  i/i.i  i  --  1 1 1 1  it  •  1 1  l  h  d  I  In  "  1 1 1  1 1  ‘  i  i  rial  t  osl  tin  All  hRU 

~.  r  i ..  I  r  c  i  i  i  al  i  - 1 !  i.-t  i  I  h  i  i  •!  olid  cm  I  tut  1. 1  1  ;d  ion  was  ret  u. 


indicated  by  the  first  approach. 


Thus,  the  first  approach 


rpf  its  in  an  overall  quantity  of  24  stations  for  the  B-1B 
fleet,  while  the  second  approach  results  in  an  overall 
quantity  of  26  stations  for  the  B-1B  fleet.  Either 
approach,  however,  is  based  on  very  preliminary  data  and 
must  be  regarded  as  tentative.  Since  the  results  are  nearly 
identical,  and  since  the  first  approach  involves  fewer 
assumptions  than  does  the  second,  the  first  approach  was 
selected  as  the  solution  technique  for  the  remainder  of  this 
research  effort.. 


Table  XIII.  Baseline  Case;  16  Aircraft  Per  Base 
Logistics  Support  Cost  Tradeoff 
FY  81  $  in  Millions 


Digital  Test  Station 
QTY  =  1  QTY  =  2  QTY  =  3 


LRU  Spares 

28.  2 

27 . 0 

26.6 

Test 

Stat i ons 

2.8 

5.6 

8.4 

Total  Cost 

3 1 . 0 

32.6 

35 . 0 

Radi  o 
QTY  =  1 

Freq  Station 

QTY  =  2  QTY  =  : 

LRU  Spares 

35 . 3 

30.  3 

30 . 3 

Test 

Stat i ons 

2.8 

5.  6 

8.4 

Total  Cost 

38.  1 

35.9 

38.  7 

85 


Digital  Analog 
Video  Test  Station 
QTY  =  1  QTY  =  2  QTY  = 


18.0  15.9  15. 


2.8  5.6  8.  ■ 

20.8  21.5  24 .  : 


Radar  EW  Station 
QTY  =  1  QTY  =  2  QTY  =  : 

22.6  18.4  18.: 


2.8  5.6  8.  ■ 

25.4  24.0  26. < 


Table  XIV.  Baseline  Case;  26  Aircraft  Per  Base 
Logistics  Support  Cost  Tradeoff 
FV  81  $  in  Millions 

Digital  Analog 

Digital  Test  Station  Video  Test  Station 

QTY  =  1  QTY  =  2  QTY  =  3  QTY  =  1  QTY  =  2  QTY  = 


LRU  Spares 

38.5 

35.  9 

35.9 

20.  1 

17.6 

17.6 

Test 

Stati ons 

2.8 

5.6 

8.4 

2.8 

5.6 

8.4 

Total  Cost 

41.3 

41.5 

44.3 

22.9 

23.2 

26 . 0 

Radio 
QTY  =  1 

Freq  Station 

QTY  =  2  QTY  =  3 

Radar 
QTY  =  1 

EW  Station 
QTY  =  2  QTY 

=  3 

LRU  Spares 

X 

45.  5 

43.  7 

X 

26.9 

'->A 

i.O.  /L 

Test 

Stati ons 

X 

5.6 

8.4 

X 

5.6 

8.4 

Total  Cost 

X 

51. 1 

52.  1 

X 

32 .  5 

34 . 6 

BTOfSciwSt? 


Wl 


Table  XV.  Baseline  Case;  o2  Aircraft  Per  Base 
Logistics  Support  Cost  Tradeo-f-f 
FY  81  $  in  Millions 

Digital  Analog 


Digital  Test 
QTY  =  1  QTY  =  2 

Stat i on 
QTY  =  3 

Vi  deo 
QTY  =  1 

Test  Station 

QTY  =  2  QTY  =  3 

LRU  Spares 

46.  1 

42.  4 

42.3 

r^r-y  -J 

i.i.1  X 

18.7 

18.  2 

Test 

Stations 

2.8 

5.6 

8.4 

2.8 

5.6 

8.4 

Total  Cost 

48.9 

48.0 

50.  7 

24.9 

24.3 

26.  6 

Radio 
QTY  =  1 

Freq  Station 

QTY  =  2  QTY  =  3 

Radar 
QTY  =  1 

EW  Station 

QTY  =  2  QTY  =  3 

LRU  Spares 

X 

53.4 

51 . 3 

X 

31.8 

31.8 

Test 

Stations 

X 

5.6 

8.4 

X 

5.6 

8.4 

Total  Cost 

X 

59.0 

59.7 

X 

37.4 

40.  2 

Q®Oi.Q¥®!=Qt  for  §  Major  Conventional  War 

All  analysis  discussed  so  -far  has  been  restricted  tc 
peacetime  maintenance  and  support.  The  peacetime  scenario 
calls  -for  aircraft  on  daily  alert  and  for  aircrew  training 
missions.  In  times  of  crisis,  however,  some  B-lBs  might  be 
dispersed  to  a  satellite  base.  Only  organisational  level 
(and  not  intermediate  level)  maintenance  would  be  performed 
at  these  bases,  thus  no  additional  test  stations  would  be 
required  to  support  such  a  dispersal  (15:26-27).  There  is 
also  no  plan  to  use  the  B-1B  ATE  to  support  operations 
during  a  sustained  nuclear  war.  Again,  only  organi zati onal 
maintenance  would  be  performed  to  support  such  operations. 
After  all,  the  B-1B  ATE  is  not  hardened  for  protection 
against  electromagnetic  pulse  (EMP)  events.  Another 
possibile  wartime  mission  would  be  a  major  conventional  war. 
Although  there  is  no  formally  documented  requirement  for 
such  a  mission,  the  potential  conventional  role  for  the  B-1B 
is  largely  undefined.  The  analysis  that  follows  is  entirely 
hypothetical  and  does  not  constitute  any  formal  SAC  plan  for 
actual  usage  of  the  B-1B. 

The  B-1B  fleet  consists  of  5  squadrons  of  16  PAA  each 
and  1  squadron  of  10  PAA  for  combat  crew  training.  In  this 
analysis,  it  was  assumed  that  a  flight  of  6  PAA  would  be 
taken  from  each  of  the  first  5  squadrons  and  that  10  PAA 
would  remain  behind  for  purposes  of  strategic  alert.  It  was 


also  assumed  that  the  10  PAA  squadron  would  also  remain  at 
its  normal  location.  Thus,  a  total  o-f  30  aircraft  would  be 
deployed  while  60  aircraft  would  remain  at  the  MOBs.  It  was 
also  assumed  that  the  30  deployed  aircraft  would  be  sent  to 
two  sites  overseas;  one  site  would  receive  18  aircraft,  and 
the  other  would  receive  12  aircraft.  The  final  allocation 
of  aircraft  is  shown  in  Table  XVI. 

Table  XVI.  Deployment  of  B-1B  Fleet  to 
Support  Major  Conventional  War 


Before 

After 

Deployment 

Deployment 

Location 

PAA  Aircraft 

PAA  Aircraft 

Mission 

MOB  #1 

26 

20 

Strategic  alert, 
training 

MOB  #2 

16 

10 

Strategic  alert 

MOB  #3 

32 

20 

Strategic  alert 

MOB  #4 

16 

10 

Strategic  alert 

Deployment 

Site  #1 

0 

18 

Conventional  war 

Deployment 

Site  #2 

0 

12 

Conventional  war 

At  the  deployment  sites,  it  was  assumed  that  the 
aircraft  would  fly  three  times  as  many  flying  hours  as  in 
peacetime,  and  that  the  aircraft  would  fly  the  same  amount 
every  day  for  seven  days  per  week.  It  was  also  assumed  that 
the  maintenance  shifts  would  be  expanded  from  two  8-hour 
shifts  per  day  to  two  12-hour  shifts  per  day,  and  that  the 
maintenance  shifts  would  be  expanded  from  five  days  per  week 
to  seven  days  per  week.  It  was  also  assumed  that  the 


reliability  and  maintainability  characteristics  of  the  LRUs 


and  ATE  test  stations  would  be  the  same  in  wartime  as  in 
peacetime.  The  model  was  adjusted  to  account  for  the 
expanded  maintenance  shifts  and  used  to  determine  the 
minimum  test  station  quantities  required  to  support  a  four 
calendar  day  base  repair  cycle  time.  The  total  requirement 
to  support  the  hypothetical  deployment  was  found  to  be  31 
test  stations.  These  results  are  shown  in  Table  XVII. 

Table  XVII.  Test  Station  Quantities  to 
Support  Hypothetical  Deployment 


Site 

Di.gi.tad 

Digital 

An  Vi^deo 

BE 

r/ew 

MOB  #1 

1 

1 

2 

4. 

MOB  #2 

1 

1 

1 

1 

MOB  #3 

1 

1 

2 

'“l 

4L 

MOB  #4 

1 

1 

1 

1 

Deployment 

Site  #1 

1 

1 

2 

'~1 

dm 

Deployment 

Site  #2 

1 

1 

1 

n 

dm 

Total 

6 

6 

9 

10 

90 


VII.  Summary  and  Recommendations 


Overview 

The  purpose  of  this  research  effort  was  to  develop  a 
technique  to  determine  B— IB  automatic  test  equipment  (ATE) 
station  quantities  required  to  support  the  B-1B  fleet  at 
base  level.  The  ATE  stations  are  essential  for  the  support 
of  over  100  B— IB  line  replaceable  units  (LRUs).  The 
decision  as  to  how  many  stations  to  procure  not  only  affects 
whether  the  B-1B  avionics  LRUs  can  be  supported,  but  also 
involves  expenditures  of  several  millions  of  dollars.  The 
remainder  of  this  section  provides  a  summary  of  the  major 
results  of  this  research  effort  and  presents  the  B-1B  System 
Program  Office  (SPO)  with  a  package  of  options  for  test 
station  acquisition  strategy.  The  summary  describes  the 
simulation  model  of  B-1B  avionics  maintenance  that  was 
developed,  describes  the  techniques  developed  to  determine 
station  quantity  requirements,  and  compares  the  simulation 
model  to  a  second  model,  based  on  analytical  queueing 
theory,  also  developed  as  part  of  the  research  effort. 

Summary 

The  most  significant  product  of  this  research  effort  is 


a  detailed  simulation  model  developed  in  the  Q-GERT 
simulation  language.  Had  this  model  not  been  developed,  the 
B-1B  SPO  would  have  had  to  rely  on  the  General  Dynamics 
station  loading  model  described  in  Section  I.  There  are 
several  deficiencies  in  the  contractor  model  which  have  been 
remedied  in  the  Q-GERT  model.  First,  the  contractor 
estimates  of  mean  LRU  maintenance  times  failed  to  account 
for  manual  troubleshooting  or  additional  testing  which 
sometimes  occur  in  real-world  maintenance.  The  contractor 
estimates  have  been  adjusted,  as  explained  in  Section  III, 
to  account  for  both  effects.  Second,  the  contractor  model 
assumes  that,  the  avionics  LRU  arrival  rate  is  constant  over 
the  maintenance  day.  However,  current  SAC  planning  suggests 
that  the  B-1B  flying  hours  are  not  spread  evenly  over 
different  times  of  the  day  or  different  days  of  the  week  as 
explained  in  Section  III.  The  Q-GERT  simulation  model  was 
developed  to  include  this  distribution.  Third,  the 
contractor  model  fails  to  include  the  technician  travel  time 
to  obtain  a  shop  replaceable  unit  (SRU)  spare  necessary  to 
complete  a  LRU  repair.  Estimates  of  these  travel  times  were 
obtained  from  HQ  SAC  personnel  and  incorporated  into  the 
Q-GERT  model  as  described  in  Section  IV.  Fourth,  the 
contractor  model  implicitly  assumes  that  SRU  spares  are 
available  on  base  1007.  of  the  time.  The  Q-GERT  model,  on 
the  other  hand,  is  more  realistic;  it  explicitly  assumes  SRU 
spares  are  available  only  907.  of  the  time.  This  is 


important  because  the  SRU  spare  stockouts  have  a  significant 
effect  on  station  avai  1  abi  1  i  ty .  The  main  point  overall  is 
that  the  contractor  model  is  consistently  optimistic  in 
estimating  the  workload  of  avionics  LRUs  on  ATE  test 
stations.  The  Q-GERT  model  was  developed  in  an  iterative 
•fashion  with  a  great  deal  of  feedback  and  constructive 
criticism  from  personnel  from  HQ  SAC  and  the  B-1B  SF'O.  This 
helped  ensure  model  realism. 

The  second  major  product  of  this  research  effort  was 
the  development  of  techniques  to  determine  the  best  choice 
of  test  station  quantities.  To  accomplish  this,  the 
simulation  model  was  expanded  to  include  organi zati onal 
(flight-line)  maintenance  and  various  admi ni strat i ve  delays 
so  that  the  model  output  is  the  (complete)  LRU  base  repair 
cycle  time.  This  is  the  time  from  LRU  failure  to  repair  and 
return  to  base  supply.  Two  approaches  were  then  developed 
in  determining  the  best  choice  of  station  quantities.  The 
first  approach  was  to  buy  sufficient  quantities  such  that 
the  LRU  base  repair  cycle  time  was  shorter  than  some 
established  standard.  This  standard  would  be  the  planned 
base  repair  cycle  time  used  to  determine  the  avionics  LRU 
spare  (pipeline  and  safety  stock)  quantities.  The  second 
approach  was  to  perform  a  cost-benefit,  analysis  on 
procurement  of  additional  test  stations  versus  procurement 
of  additional  LRU  spares.  This  approach  compared  the  costs 
of  additional  test  stations  (and  the  benefits  of  shorter 


base  repair  cycle  times)  to  the  costs  of  additional  LRU 
spares  (and  the  benefits  of  -fewer  test,  stations).  The  two 
approaches  yield  nearly  identical  results. 

A  second  model  o-f  B-1B  avionics  maintenance  was  also 
developed  during  this  research  e-f-fort.  The  analytical  model 
is  based  on  the  theory  o-f  queueing  networks  and  other 
queueing  techniques.  The  analytical  model  still  required 
many  more  crude  and  simplifying  assumptions  than  did  the 
simulation  model,  and  the  results  obtained  from  the 
analytical  model  were  not  very  sati sf actory.  The  reasons 
for  the  poor  results  (and  recommendations  for  future 
research  which  might  correct  these  problems)  are  explained 
in  Section  V.  The  Q-GERT  simulation  model  was  selected  as 
the  best  choice  for  future  use  by  the  B-1B  SPO. 

B§£9(Dffl§Qded  Management  Opt  .ions  for  the  B- 1.B  SPO 

The  first  choice  to  be  made  concerns  the  selection  of 
reliability  and  maintainability  inputs.  Relying  on  the 
associate  contractor  estimates  for  LRU  reliability  and  the 
ATE  vendor’s  estimates  for  LRU  performance  test  times  (as 
defined  in  Section  III)  is  financially  prudent  since  these 
inputs  lead  to  the  minimum  required  quantities.  However, 
there  is  some  risk  for  B-1B  avionics  supportabi 1 i ty  since 
these  estimates  might  be  excessively  optimistic.  If  this 


were  to  happen,  the  station  quantities  would  not  be 
sufficient  to  support,  the  avionics  workload.  One  course  o-f 
action  would  be  to  accept  the  contractor  estimates  -for  the 
time  being,  and  therefore  make  a  tentative  decision  to  buy 
only  the  minimum  number  of  stations.  With  this  course  of 
action,  the  B-18  SPO  should  still  retain  some  flexibility  to 
obtain  additional  test  stations  should  the  need  arise.  This 
flexibility  could  be  preserved  by  a  provision  for  separately 
priced  contractual  options  for  additional  test  stations. 
Another  way  to  preserve  flexibility  would  be  to  procure 
stations  collectively  for  i  ritermedi  ate  level  maintenance 
(which  was  addressed  in  this  research  effort),  depot  level 
maintenance,  contractor  use,  and  maintenance  training  (which 
were  not  addressed  in  this  research  effort)  with  a  few 
additional  stations  designated  for  management  reserve.  A 
second  course  of  action  would  be  to  procure  test  stations 
based  on  the  more  pessimistic  reliability  and 
mai ntai nabi 1 i ty  inputs  explained  in  Section  III.  This 
course  of  action  involves  less  risk  for  B-1B  avionics 
supportabi  1  i  ty ,  but.  also  leads  to  greater  quantities  of 
stations  and  therefore  larger  expenditures.  Some  of  these 
expenditures  might  turn  out  to  have  been  unjustified  if  the 
baseline  contractor  estimates  turn  out  to  be  reasonable. 

The  first  course  of  action,  therefore,  is  probably  the 
wisest  as  long  as  the  B-1B  SPG  is  able  to  preserve 
flexibility  to  procure  additional  stations  as  required. 


Assuming  that  the  B--1B  h>FU  indeed  retains  some 
flexibility  for  a  range  of  station  quantities,  the  B--1B  BF'U 
should  periodically  review  and  refine  its  estimates  of  test, 
station  requirements.  It  is  likely  that,  over  time,  better 
quality  estimates  of  the  reliability  and  maintainability 
inputs  will  become  available.  Presumably,  the  associate 
contractors  estimates  should  become  more  refined  as  the  B--1E- 
avionics  go  through  various  reliability  qualification  tests 
and  maintainability  demonstrations.  A  review  of  test, 
station  requirements  could  take  place  periodically  (say 
every  six  months)  by  running  the  Q—6ERT  simulation  model 
with  the  most  recent  information  available.  In  addition,  it 
might  be  desirable  to  task  staff  organizations  (such  as 
ASD/EN  or  AFALC/PT)  for  assistance  in  determining  that  the 
contractor  estimates  are  reasonable,  or  in  obtaining  better 
estimates  if  needed. 

The  second  choice  to  be  made  is  the  methodology  to 
determine  test  station  quantities.  As  explained  in  Section 
VI,  the  first,  method  is  to  buy  sufficient  test  stations  to 
support  a  four  calender  day  base  repair  cycle  time;  the 
alternative  second  method  is  to  buy  sufficient  test  stations 
to  achieve  the  lowest  possible  logistics  support  cost 
(including  the  cost  of  ATE  test  stations  and  avionics  LRU 
spares) .  There  are  many  issues  to  be  considered  before  such 
a  choice  could  be  made.  The  first  issue  is  the  potential 
impact  of  funding  constraints.  For  example,  if  Ei-1B  spares 


funding  were  to  be  severely  constrained  in  the  future,  it 
would  tend  to  favor  the  alternative  second  method  (lowest 
logistics  support  costs).  This  second  method  leads  to 
larger  station  quantities,  shorter  repair  cycle  times,  and 
therefore  reduced  LRIJ  spares  requirements.  Qn  the  other 
hand,  if  B-1B  peculiar  support  equipment  (which  includes 
ATE)  funding  were  to  be  severely  constrained  in  the  future, 
it  would  tend  to  favor  the  first  method  (support  a  four  day 
base  repair  cycle  time)  since  it  leads  to  lower  test  station 
quantities.  The  second  issue  to  be  considered  is 
coordination  with  other  organi z at i ons .  The  first  method 
(support  a  four  day  base  repair  cycle  time)  could  be 
accomplished  unilaterally  by  the  B-1B  SF’O.  The  alternative 
second  method  would  require  an  integrated  ATE/LRU  spares 
acquisition  strategy.  Both  ATE  procurement  and  avionics  LRU 
spares  procurement  would  have  to  be  based  on  the  same 
logistics  factors  (LRU  reliability  and  base  repair  cycle 
time)  to  actually  achieve  the  lowest  logistics  support  cost. 
Accompl i shi ng  such  an  integrated  strategy  would  require 
close  coordination  between  the  SPO,  the  B-1B  System  Manager, 
and  many  equipment  specialists  and  item  managers  located  at 
many  ALCs. 

The  third  choice  to  be  made  concerns  buying  test, 
station  quantities  to  support  peacetime  or  wartime.  One 
approach  to  this  issue  would  be  to  only  procure,  at  least 
•for  the  time  being,  test  stations  to  support  peacetime  only. 


However,  analysis  indicates  that  peacetime  quantities,  even 
with  expanded  maintenance  shifts,  will  not  be  sufficient  to 
support  a  major  conventional  deployment.  Should  the 
requirement  for  a  conventional  deployment  ever  emerge, 
additional  test  stations  would  have  to  be  borrowed  (from 
depot  or  other  operational  sources)  or  procured  at  that 
time.  An  alternative  approach  would  be  to  develop  the  most 
likely  conventional  scenario  and  procure  sufficient  test 
stations  to  support  such  a  scenario.  It  would  be  difficult 
to  justify  adapting  this  alternative  approach  unless  the 
requirement  for  a  conventional  deployment  were  formally 
documented  in  an  official  planning  document  such  as  the  B-1B 
Program  lianagemet  Directive  (PMD)  . 

Even  if  the  O-GERT  simulation  model  were  perfectly 
accurate,  there  is  simply  too  much  uncertainty  regarding  the 
reliability  and  maintainability  inputs  to  determine  the 
precise  base  level  test  station  requirements.  Rather,  at 
this  time,  it  is  only  possible  to  determine  a  reasonable 
range  of  quantities.  The  actual  product  of  this  research 
effort  is  therefore  not  the  final  estimate  of  test  station 
quantities,  but  rather  a  simulation  model  with  which  other 
techniques  can  be  used  by  the  B-1B  SPO  over  the  next  several 
years  to  continually  update  and  refine  its  estimate  of  test 
station  quantities.  In  addition,  the  B-1B  SPD  can  now 
select  from  a  variety  of  approaches  (which  were  described  in 
Sections  IV  and  VI)  to  determine  the  best  assessment  of 


station  quantities.  The  selection  of  the  approach  to  this 
problem  is  a  matter  of  management  judgement,  which  of  course 
should  be  reserved  for  the  B--1B  SF‘Q.  A  summary  of  the 
results  for  the  various  approaches  is  shown  again  for  review 
in  Table  XVIII. 


Table  XVIII.  Summary  of  Test  Station  Quantities 
to  Support  B-1B  Operational  Fleet  at  Base  Level 


Approach 

Digital 

Di gi tal 
Analog 
Video 

rf 

r/ew 

Total 

Basel i ne 

4 

4 

8 

8 

24 

Lowest  Logistics 
Support  Cost 

Id 

5 

8 

8 

26 

HQ  SAC  MTBD 

Est i mates 

5 

8 

10 

10 

33 

R I  Performance 

Test  Time  Estimates 

5 

6 

9 

9 

29 

Conventi onal 
Deployment 

6 

6 

9 

10 

31 

Appendix  A: 


Simulation  Model  User’s  Guide 


of  iDEyts 

The  -following  inputs  are  stored  in  FORTRAN  data 
statements  in  FORTRAN  Subroutine  UI.  The  ATE  test  station 
inputs  are  one-dimensional  arrays  of  size  four;  one  entry 
corresponds  to  each  station  type.  The  inputs  must  be  stored 
in  the  correct  order  (Digital,  Digital  Analog  Video,  RF,  and 
REW) .  The  ATE  test  station  inputs  are: 

TSMTBF  -  Test  Station  Mean  Time  Between  Failures 
ATEQTY  -  Quantity  of  stations  by  station  type 
ATERF'l  -  Average  time  for  setup  and  fault  isolation 
ATERF'2  -  Average  time  for  station  repair  and  teardown 
The  LRU  inputs  are  one-dimensional  arrays  of  (up  to) 
size  200.  Each  entry  in  each  array  corresponds  to  one  LRU. 
The  inputs  may  be  stored  in  any  order,  but  the  order  must  be 
the  same  for  all  arrays.  The  LRU  inputs  are: 

XMTBD  -  Mean  Time  Between  Demand  (in  flying  hours) 

QPA  antity  per  aircraft 

STNRE  .  Station  Requirement: 

)  =  Digital 

j  =  Digital  Analog  Video 
0  =  Radio  Frequency 
—  4.0  =  Radar  Electronic  Warfare 


AFSC 


—  Intermediate  Level  Technician  Requirement 

—  1.0  -  Communication  Navigation 

—  2.0  =  Auto  Fit  Controls  Instruments 

—  3.0  =  CJ-f tensive  Avionics 

—  4.0  =  Defensive  Avionics 

SETUP  -  Time  to  set-up  LRU  on  station 
TSSITA  -  Time  to  run  ITA  confidence  check 

WARMUP  -  Time  to  warm  up  LRU 

F'ERFT  -  Time  to  run  one  complete  performance  test 

FRRTS  -  Fraction  of  maintenance  actions  that  are 

Repairable  This  Station 

FRRTOK  -  Fraction  of  maintenance  actions  that  are 
Retest  Okay 

FRNRTS  -  Fraction  of  maintenance  actions  that  are 
Not  Repairable  This  Station 

TEARDN  -  Time  to  tear-down  LRU  from  station 

REPAIR  -  Time  to  repair  LRU  (remove  and  replace 
failed  SRIJ)  . 

Other  inputs  stored  in  FORTRAN  data  statements  in  UI 


ACBASE  -  Number  of  aircraft  per  base 

SRUAVL  -  Fraction  of  the  time  that  a  required  SRU 
spare  will  be  available  on  base 

FHACMO  -  Flight  hours  per  aircraft  per  month 

FLYPER  -  Time  length  of  interval  over  which  sorties 
are  launched 

FACTDY  -  Fraction  of  weekly  flying  hours  by  day  of  the 
week;  an  array  of  size  five 

FACTTM  -  Fraction  of  flying  hours  by  time  of  day;  an 
array  of  size  two  (for  morning  sorties  and 
night  sorties) 

NLRU  --  Number  of  LRUs;  must  be  less  than  or  equal 


Conversion  of  Flying  Hours  to  Real.  lime 


The  LRU  MTBD  values  must  be  input  in  -flying  hours. 
FORTRAN  Subroutine  UI  takes  the  reciprocal  of  each  MTBD  to 
find  the  LRU  arrival  rate.  These  LRU  arrival  rates  are  then 
added  to  determine  the  aggregate  arrival  rate.  The 
reciprocal  of  this  aggregate  arrival  rate  is  the  aggregate 
MTBD  in  flying  hours;  this  is  FORTRAN  variable  DEMFH  in  the 
model.  The  aggregate  MTBD  is  then  converted  to  real  time  by 
the  following  equation: 


DEMSH  =  _ DEMFH  *  < FLYF'ER  ♦  4.533) _ 

FACTDY ( I AT  1 )  *  FACTTM ( IFN)  *  ACBASE  *  FHACMC 


4.333  is  the  number  of  weeks  per  month.  I ATI  is  an  integer 


from  1  through  5  to  denote  the  day  of  the  week.  IFN  i s  an 
integer  (either  1  or  2)  to  denote  the  time  of  day  (morning 


or  night).  The  time  interval  (FLYF'ER)  for  LRU  arrivals  is 
assumed  to  be  4.0  hours  as  explained  in  Section  IV. 


How  to  Run  the  Model 


The  simulation  model  is  run  in  batch  mode  on  the  ASD 
CDC  Cyber.  For  user  convenience,  the  model  will  be 
furnished  in  three  versions  for  the  16  F'AA,  26  F'AA,  and  32 
F'AA  cases.  These  models  are  named  M0DEL16,  MQDEL.26,  and 
MODEL-32.  Output  must  be  obtained  in  the  computer  room  in 
AFIT/EN,  Bldg  640  or  wherever  the  output  is  designated.  The 


batch  command  for  AFIT/EN  would  be:  BA, MODEL  16, AF . 


How  to  Change  Station  Quantities 

Unfortunately,  changing  station  quantities  requires 
changes  to  both  FORTRAN  and  Q-GERT  code.  The  quantities 
must  be  input  in  the  correct  order  (Digital,  Digital  Analog 
Video,  RF  and  R/EW) .  The  FORTRAN  code  which  must  be  changed 
is: 

DATA  ATEQTY /1.0, 1.0, 1.0, 1.0/ 

Each  entry  represents  the  number  of  stations  at  the  base  by 
station  type. 

The  Q-GERT  lines  which  must  be  changed  are: 

RES, 2/DIGITAL, 1,66,28* 

RES, 3/DIGAV, 1,67,29* 

RES , 4/RF , 1,68,30* 

RES, 5/RAD AREW, 1,69,31* 

ALT, 11 2, D, 2,— 1, 66 , 28* 

ALT, 1 13, D, 3, -1,67, 29* 

ALT, 1 14, D, 4, -1,68, 30* 


ALT, 

1 

15, 

D, 

5, 

— 

1 

,69 

,31 

ALT, 

1 

30, 

D, 

o 

•*-  f 

1 

•I 

66 , 

28* 

ALT, 

1 

31, 

D, 

1 

67, 

29* 

ALT, 

1 

32, 

D, 

4, 

1 

1 

68, 

30* 

ALT, 

1 

33, 

D, 

5? 

1 

3 

69, 

31* 

There  are  thus  twelve  lines  of  Q-GERT  code  which  must  be 
changed  to  change  input  test  station  quantities.  For  the 
RESOURCE  lines  (which  start  RES,),  the  quantity  of  stations 
is  the  third  entry  (1  in  this  example) .  For  the  ALTER  node 
lines  (which  start  ALT,),  the  quantity  of  stations  is  the 


fifth  entry.  For  the  ALTER  nodes  112  through  115,  the  test 
station  quantities  are  input  as  negative  integers  to  denote 
the  end  of  a  maintenance  shift.  For  the  ALTER  nodes  130 
through  133,  the  test  station  quantities  are  input  as 
positive  integers  to  denote  the  start  of  a  maintenance 
shift. 

The  remainder  of  this  appendix  consists  of  the  network 
D— GERT  diagram  and  the  model  listing. 
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LIST, ALL 


1 00=LMR , CM 1 77000 , T400 , I 0200 .  T830662 , ROARK , 4533 . 

1 10=ATTACH, IMSL, IMSL, ID=LIBRARY, SN=ASD. 

120=LIBRARY, IMSL. 

1 30=ATT ACH , PROOF  I L , QGERTPROC , I D= AF I T . 

140=FTN5, ANSI=0. 

150=BEGIN, QGERT, PROOF IL, M=LGO, MODE=X . 

160=*E0R 

170= 

1 80=  FUNCT ION  UF ( I FN > 

190= 

200=  COMMON / QV AR / NDE, NFTBU  <500) , NREL (500) , NRELP  <500) , 

210=  1  NREL2  <  500  >  ,  NRUN ,  NRUNS ,  NTC  <  500 )  ,  F'ARAM  <100,4),  TBEG ,  TNOW 

220= 

230=  COMMON/UCOM1  /DEMFH,  FLYF'ER,  ACBASE,  FHACMO, 

240=  1  FACTDY  <5) , FACTTM  <2) 

250= 

260=  COMMON /UC0M2/CUMPRB <200) , STNREQ <200) ,AFSC<200) ,FRRTS<200) , 

270=  1  FRRTOK <  200 ) , FRNRTS <  200 ) , SETUP  <  200  > , TSS I TA  <  200 ) , 

280=  2  WARMUP  <200) , PERFT  <200) , REPAIR  <200) , TEARDN  <200) 

290= 

300=  COMMON/UCOM3/DSEED, SRUAVL 

310= 

320=  DOUBLE  PRECISION  DSEED 

330=  DSEED  =  DSEED  +  < 100000. *NRUN> 

340= 

350=C  ATTRIBUTE  1  IS  DAY  OF  THE  WEEK 

360=  ATI  =  GATRB  < 1 ) 

370=  I ATI  =  ATI 

380= 

390=C  CONVERT  RELIABILITY  IN  FLYING  HOURS  TO 

400=C  RELIABILITY  IN  ACTUAL  HOURS 

410=  DEMSH  =  <  DEMFH*FL YPER*4 . 333 ) / <  FACTDY  < I AT 1 ) * 

420=  1  FACTTM  < IFN) *ACBASE*FHACMO) 

430= 

440=C  UN IF  IS  UNIFORM <0,1) 

450=C  UF  IS  SIMULATED  EXPONENTIAL  ARRIVAL  TIME 

460=  UN IF  =  GGUBFS  < DSEED ) 

470=  UF  =  - 1 . 0*DEMSH* ALOG  <  UN I F  > 

480=  IF  <UF  .GT.  16.0)  UF  =  16.001 

490= 

500=C  DETERMINE  WHICH  LRU  HAS  FAILED  BY  SECOND  UNIF<0,1) 

510=  UN IF  =  GGUBFS < DSEED) 

520=  I  =  1 

530=  IF  <UNIF  .LT.  CUMPRB  <  I )  )  GO  TO  HO 

540=  100  1=1+1 

550=  J  =  I  -  1 

560=  I F  <  <  CUMPRB  <  J ) . LE . UN IF)  .AND.  <UNIF. LT. CUMPRB < I )) )  GO  TO  110 

570=  GO  TO  100 

580=  110  CONTINUE 

590= 

600=C  ATTRIBUTE  1  —  STATION  REQUIREMENT 
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I 


1 


s 


610* 
620= 
630= 
640=C 
650= 
660= 
670= 
680=C 
690= 
700= 
710= 
720=C 
730=C 
740= 
750= 
760= 
770= 
780= 
790= 
800= 
810= 
820=C 
830= 
840=C 
850= 
860= 
870 =C 
880= 
890= 


A1  =  STNREQ ( I ) 
CALL  PATRB ( A 1 , 1 ) 


ATTRIBUTE  2  —  T 
A2  =  AFSC(I) 

CALL  PATRB ( A2, 2) 


-  TECHNICIAN  (AFSC)  REQUIREMENT 


THREE  FOR  LRU  REPAIR 


(RTS, RTOK,  OR  NRTS) 


ATTRIBUTE  4  —  7 

A4  =  3.0 

CALL  PATRB  <A4, 4) 


DETERMINE  TYPE  OF  LRU  REPAIR 
BY  UN IF (O,  1) 

UN IF  =  GGUBFS ( DSEED) 

XLOW  =  FRRTS(I) 

XMID  =  FRRTS(I)  +  FRRTOK < I ) 
XHIGH  =  XMID  +  FRNRTS ( I ) 

IF  (UN  IF  .LT.  XLOW)  GO  TO  120 
IF (UN IF  .LT.  XMID)  GO  TO  130 
IF  <UNIF  .LE.  XHIGH)  GO  TO  140 

IF  RTS  (REPAIR  THIS  STATION) 
120  A5  =  1.0 

ATTRIBUTE  5  —  1.0  FOR  RTS 
CALL  PATRB  ( A5, 5) 


SUM  IS  MEAN  REPAIR  TIME  FOR  ITH  LRU 

SUM1  =  SETUP! I)  +  TSSITA(I)  ■+•  WARMUP ( I )  +  0.5*PERFT(I> 
SUM2  =  (0.75)*REPAIR(I)  +  (O. 25) *2. 0*REPAIR ( I ) 


900= 

1  + 

PERFT(I)  +  TEARDN ( I ) 

910= 

SUM  =  SUM1  +  SUM2 

920= 

930=C 

SIMULATE 

ERLANG  WITH  K  =  2 

940= 

ERLANG  = 

0.0 

950= 

DO  125  K 

-  1,2 

960= 

UN  IF  = 

GGUBFS (DSEED) 

970= 

DRWEXP 

=  —  1 . 0* (SUM/ 2. 0) *ALOG (UN IF) 

980= 

ERLANG 

=  ERLANG  +  DRWEXP 

990= 

125  CONTINUE 

1000= 

1010=C 

1020= 

1030= 

1040= 

1 050=C 
1060= 
1070= 
1080= 
1090=C 
1100= 
1110= 
1120= 
1130= 
1140= 


ATTRIBUTE  6  —  TIME  FOR  FIRST  PART  OF  REPAIR 
A6  =  (SUMl/SUM) *ERLANG 
CALL  PATRB (A6, 6) 

ATTRIBUTE  7  —  TIME  FOR  SECOND  PART  OF  REPAIR 
A7  =  (SUM2/SUM) *ERLANG 
CALL  PATRB (A7, 7) 

DETERMINE  IF  SRU  SPARE  WILL  BE  AVAILABLE 
UN IF  =  GGUBFS (DSEED) 

I F ( UN I F  .LT.  SRUAVL)  AS  =  1.0 
IF (UN IF  .GE.  SRUAVL)  A8  =  0.0 
CALL  PATRB (A8, 8) 
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RETURN 


1150= 

1160= 

1 170=C  IF  RETEST  OKAY  (RTOK) 

1180=  130  A5  =  2.0 

1190=  CALL  PATRB  <  A5,  5) 

1200= 

1210=C  SUM  IS  MEAN  RTOK  TIME  FOR  ITH  LRU 

1220=  SUM1  =  SETUP ( I )  +  TSSITA(I)  +  WARMUP < I )  +  PERFT ( I ) 

1230=  SUM2  =  0.3*PERFT(I)  +  TEARDN  K I ) 

1240=  SUM  =  SUM1  +  SUM2 

1250= 

1260=C  SIMULATE  ERLANG  WITH  K  =  5 

1270=  ERLANG  =  0.0 

1280=  DO  135  K  =  1,5 

1290=  UN IF  =  G6UBFS  <  DSEED) 

1300=  DRWEXP  =  — 1 . O* (SUM/ 5. O) *ALOG (UNIF) 

1310=  ERLANG  =  ERLANG  +  DRWEXP 

1320=  135  CONTINUE 

1330= 

1340=C  ATTRIBUTE  6  —  TIME  FOR  FIRST  PART  OF  REPAIR 

1350=  A6  =  ( SUM 1 /SUM) *ERLANG 

1 360=  CALL  PATRB <  A6 , 6 ) 

1370= 

1380=C  ATTRIBUTE  7  —  TIME  FOR  SECOND  PART  OF  REPAIR 

1390=  A 7  =  (SUM2/SUM) *ERLANG 

1400=  CALL  PATRB  <  A7 ,  7 ) 

1410= 

1420=C  ATTRIBUTE  8  —  SRU  SPARE  NOT  NEEDED  FOR  RTOK 

1430=  A8  =  2.0 

1440=  CALL  PATRB  ( A8, 8) 

1450= 

1460=  RETURN 

1470= 

1480=C  IF  NOT  REPAIRABLE  THIS  STATION  (NRTS) 

1490=  140  A5  =  3.0 

1500=  CALL  PATRB (A5, 5) 

1510= 

1520=C  SUM  IS  MEAN  NRTS  TIME  FOR  ITH  LRU 

1530=  SUM1  =  SETUP ( I )  +  TSSITA(I)  +  WARMUP ( I )  +  0.5*PERFT(I) 

1540=  SUM2  =  TEARDN ( I ) 

1550=  SUM  =  SUM1  +  SUM2 

1 560= 

1570=C  SIMULATE  ERLANG  WITH  K  =  5 

1580=  ERLANG  =0.0 

1590=  DO  145  K  =  1,5 

1600=  UNIF  =  GGUBFS (DSEED) 

1610=  DRWEXP  =  — 1 . 0* (SUM/ 5. O) *ALOG (UNIF) 

1620=  ERLANG  =  ERLANG  +  DRWEXP 

1630=  145  CONTINUE 

1640= 

1650=C  ATTRIBUTE  6  —  TIME  FOR  FIRST  PART  OF  NRTS 

1660=  A6  =  (SUM1 /SUM) *ERLANG 

1670=  CALL  PATRB (A6, 6) 


TOWVTOT 1* .  i'.  r  *  * A’:j  V  **  v 


k:v 


IV 


1690=C 

1700= 

1710= 

1720= 

1730=C 

1740= 

1 750= 
1760= 
1770= 
1780= 
1790= 

1 800= 
1810= 
1820= 
1830=C 
1840= 
1850= 
1860= 
1870= 
1880= 
1890= 
1900= 
1910= 

1 920=C 
1 930=C 
1940= 

1 950= 
1960= 
1970= 
1980=C 
1990= 
2000= 
2010= 
2020=C 
2030=C 
2040=C 
2050=C 
2060=C 
2070= 
2080= 
2090= 
2100= 
21 10= 
2120= 
2130= 
2140= 
2150= 

2 1 60=C 
2170= 
2180= 
2190= 
2200=C 
2210= 
2220= 


ATTRIBUTE  7  —  TIME  FOR  SECOND  PART  OF  NRTS 
A7  =  (SUM2/SUM) *ERLANG 
CALL  PATRB ( A7 , 7 ) 

ATTRIBUTE  8  —  SRU  SPARE  NOT  NEEDED  FOR  NRTS 

A8  =  2.0 

CALL  PATRB (A8, 8) 

RETURN 

END 


SUBROUTINE  US ( ISN, DT IM> 

ISN  IS  DUMMY  ARGUMENT:  DT IM  IS  TIME  BETWEEN  ATE  FAILURES 

C0MM0N/QVAR/NDE,,NFTBU(500)  ,  NREL  (500)  ,  NRELP(500)  , 

1  NREL2 <  500 )  , NRUN , NRUNS , NTC ( 500 ) , P ARAM ( 1 00 , 4 ) , TBEG , TNOW 

COMMON /UC0M3/ DSEED, SRUAVL 

C0MM0N/UC0M4/ATEREL,  ATEPRB  (4  )  ,  ATERP 1  (4)  ,ATERP2(4> 

ATEREL  IS  AGGREGATE  STATION  RELIABILITY 

ATEPRB  ( I )  IS  CUM  7.  RELI ABLILTY  FOR  ITH  STATION  TYPE 

DOUBLE  PRECISION  DSEED 
DSEED  =  DSEED  +  ( 100000. *NRUN> 

SIMULATE  TIME  BETWEEN  ATE  FAILURES 
UN IF  =  GGUBFS (DSEED) 

DTIM  =  - 1.0* ATEREL* ALOG (UN IF) 

DETERMINE  WHICH  STATION  TYPE  FAILED 

1  =  DIGITAL 

2  =  DIGITAL  ANALOG  VIDEO 

3  =  RADIO  FREQ 

4  =  RADAR  EW 

UN IF  =  GGUBFS (DSEED) 

I  =  1 

IF (UN IF  .LT.  ATEPRB ( 1 ) )  GO  TO  220 
210  1=1+1 
J  =  I  -  1 

IF( (UNIF.GE. ATEPPB(J) )  .AND.  (UNIF. LT. ATEPRB ( I )) )  GO  TO  220 
GO  TO  210 
220  CONTINUE 

ATTRIBUTE  1  —  TYPE  OF  STATION 
A 1  =  I 

CALL  PATRB (A 1,1) 

ATTRIBUTE  2  —  TECHNICIAN  REQUIREMENT  ( PMEL ) 

A2  =  5.0 

CALL  PATRB (A2, 2) 
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2230= 

2240=C  ATTRIBUTE  4  —  1.0  FOR  ATE  MAINTENANCE 

2250=  A4  =  1.0 

2260=  CALL  PATRB(A4,4) 

2270= 

2280=C  ATTRIBUTE  5  —  1.0  FOR  ATE;  NO  RT OKS  OR  NRTS  ASSUMED 

2290=  A5  =  1.0 

2300=  CALL  F'ATRB  <  A5, 5) 

2310= 

2320=  SUM1  =  ATERP1 ( I ) 

2330=  SUM2  =  ATERP2  (  I ) 

2340=C  SUM  IS  MEAN  REPAIR  TIME  FOR  ITH  STATION 

2350=  SUM  =  SUM 1  +  SUM2 


2360=C 

SIMULATE 

ERLANG  WITH  K  =  2 

2370= 

ERLANG  = 

0.0 

2380= 

DO  225  K 

=  1,2 

2390= 

UN  IF  = 

GGUBFS (DSEED) 

2400= 

DRWEXP 

=  -1.0* (SUM/2. 0)*AL0G(UNIF) 

2410= 

ERLANG 

=  ERLANG  +  DRWEXP 

2420= 

225  CONTINUE 

2430= 

2440=C 

ATTRIBUTE 

6  —  TIME  FOR  FIRST  PART  OF 

2450= 

A6  =  ( SUM 1 /SUM ) *ERLANG 

2460= 

CALL  PATRB (A6, 6) 

2470= 

2480=C  ATTRIBUTE  7  —  TIME  FOR  SECOND  PART  OF  REPAIR 

2490=  A7  =  (SUM2/SUM) *ERLAN6 

2500=  CALL  PATRB <  A  7 , 7 ) 

2510= 

2520=C  DETERMINE  SRU  SPARE  PART  AVAILABILITY 

2530=  UN IF  =  GGUBFS  < DSEED ) 

2540=  IF (UN IF  . LT.  SRUAVL)  A8  =  1.0 

2550=  IF (UN IF  .GE.  SRUAVL)  A8  =  0.0 

2560=  CALL  PATRB (A8, 8) 

2570= 

2580=  RETURN 

2590=  END 

2600= 

2610= 

2620=  SUBROUTINE  UI 

2630= 

2640=  C0MM0N/UC0M1 /DEMFH, FLYPER, ACBASE, FHACMO, FACTDY (5) , FACTTM (2) 

2650= 

2660=  CQMMQN/UC0M2/CUMPRB (200) , STNREQ (200) , AFSC(200) ,FRRTS<200) , 

2670=  1  FRRTOK  <  200) , FRNRTS ( 200) , SETUP (200) , TSSITA (200) , 

2680=  2  WARMUP ( 200 ) , PERFT <  200 ) , REPA I R ( 200 ) , TEARDN ( 200 ) 

2690= 

2700=  COMMON/UCOM3/DSEED, SRUAVL 

2710= 

2720=  COMMON /UC0M4/ATEREL , ATEPRB (4 ) , ATERP 1 (4) ,ATERP2<4) 

2730= 

2740=  REAL  TSMTBF (4) , ATEQTY (4) , XMTBD (200) , QFA (200) 

2750= 

2760=  DOUBLE  PRECISION  DSEED 


2780= 

2790= 

2800= 

2810= 

2820= 

2830= 

2840= 

2850= 

2860= 

2870= 

2880= 


2900= 

2910= 


2930= 

2940= 

2950= 

2960= 

2970= 
2980= 
2990= 
3000= 
3010= 
3020= 
3030= 
3040= 
3050= 
3060= 
3070= 
3080= 
3090= 
3 1 00= 
3110= 
3120= 
3 1 30= 
3140= 
3150= 
3160= 
3170= 
3180= 
3190= 
3200= 
3210= 
3220= 
3230= 
3240= 


DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 


DATA 


DATA 

DATA 


DATA 


TSMTBF/256. , 145. , 167. , 130. / 

ATEQTY/1.0, 1.0, 1.0, 1.0/ 

ATERP 1/4*0. 5/ 

ATERP2/4*0. 5/ 

ACBASE/ 16. 0/ 

SRUAVL/0. 9/ 

FHACMO/29. 2/ 

FLYPER/4 . 0/ 

FACTDY/0 . 14,2*0.29,0. 19,0.09/ 

FACTTM/2*0. 50/ 

XMTBD/ 16458. ,1199. ,707. , 1488. ,2024. , 1576. ,952. , 100000 

1923. .  11731. . 151. .560. .590. .780. .2260. .5690. .330. . 101 

620. . 740. . 1800. .460. .232. .272. .3803. .4218. .424. . 1180. 

1923. .  2085. . 446. . 2175. . 1931 . . 3653. . 3538. . 372. .2019. . 3 

13600. .  100000. . 1530. . 41 . . 14223. . 8509. . 1391 . . 2085. . 861 

9667. . 5856. .967. . 5186. . 12171 . .5877. .6385. .3176. .5037. 

2203. .  36000. . 3030. . 198. . 289. . 478. . 282. . 418. . 1391 . . 437 

245. .  1597. . 1597. . 221 . . 7838. . 847. . 1380. . 3030. . 3030. , 

159. .  225. . 208. . 486. . 191 . . 269. . 225. . 156. . 109. . 344. . 185 

157. . 223. . 136. .111. .265. . 198. . 136. . 159. .58. .449. .449. 

422. . 422. .422. .216. .216. . 175. . 175. .50. .2906. .82. , 

706. . 4480. / 

QF'A/31 .  ,  3*2.  ,  1 .  ,  3*2.  ,  5.  ,  1 .  ,  3.  ,  8.  ,  1 .  ,  3.  ,  2.  ,  1 .  ,  1 .  ,  3.  , 
8*1 .  ,  2*1 .  ,  2. , 1 . , 2*2. , 1 . , 3. , 2*1 . , 2. , 4*1 . , 2. , 1 . , 3. . 

2*1. ,2. ,5*1. ,2. ,5*1. ,2. ,7*1. ,4*2. , 10. , 1 !, , 2. , 1 . , 3. , 

1 . .  2. . 3. . 3*3. , 7*1 . , 2*3. , 1 . , 2. , 2*1 . , 3. , 1 . , 2. , 1 . , 

1  1  doi  o  o  i  1/ 

A  ■  <y  i  ■  y  A.  a  j  X  a  j  A-  a  j|  X.  ■  f  Xa  J  X.  «  j|  X.  •  y  X  a  y  Xa  / 

STNREQ/4 1 *2 . 0, 34*1 . 0, 18*4. 0, 16*3. 0/ 

AFSC/3. 0, 6*2. 0, 3*3. 0,4. 0, 9*3. 0,4.0, 1 . 0, 4*3. 0, 7*2. 0, 
3*3. 0, 2. 0, 4. 0, 4*3.0, 13*2.0, 3*3. 0, 6*4 . 0, 5*2. 0, 2*3. 0, 
2.0, 19*4.0,2*3.0, 14*4.0,2*3.0/ 


ETUP/2*. 

y 

2*. 33, .28 

•1 

.  18 

J 

a 

3,  . 

18, 

A. 

*.2 

8 

2*. 

33, 

Ox 

A.  «  .  . 

33, 

.  18,  . 

33  ? 

2*. 18, .28 

!» 

a  •-»  -3 

J 

. 

8,  . 

33, 

a 

28, 

a 

1 

8,  . 

28, 

a  •33 

*.2 

8,3*. 

33, 

.28, .23, . 

3,  - 

1 

8 

.  33 

,  2* 

a 

18, 

n 

A. 

* 

a 

,  3* 

.  23 

18, 

.23,  . 

18, 

n  u.  n  u 

A.  “  a  A-  •_»  y  X-  “  a 

1 

8,2 

* 

a 

i 

8,  . 

23, 

a 

28, 

a 

8,  . 

18, 

.  18 

33, 

.33,  . 

28, 

a  33  a  28  y  a 

/-> 

A 

3,  . 

j 

.  28 

,  -  1 

8 

iy  a  Al 

.  j» 

* 

a  33 

•1  a 

3? 

23, 

5*.  18 

,  .2 

.  -3  y  •  A  •.-*  y  a 

8,  . 

1 

8 

!» 

5*. 

28, 

-3 

*  a  3 

•3 

2*. 

28, 

•  X  •—"*  y  vJ  *  a  X  D  y  a  X-  J|  uJ  *  «  A.  <|  a  A  O  i|  a  X  D  y  wJ  *  a  A  C  j  y  *  a  y  xl.  *  a  A.  O  y 

6  7*. 18, 2*. 23, .23,2*. 18, .28, . 18/ 

DATA  TSSITA/ 109*0. 17/ 

DATA  WARMUP/ 109*0.0/ 

DATA  F'ERFT /  .  20,  .  59 ,  1 . 44  ,  .  37  ,  1 . 08,  .29,  1 . 31 ,  .  16,  .  75,  .61  ,  .  05, 
1  . 83 , . 55 , . 55 , . 75 , . 74 , . 95 , . 95 , . 89 , 1 . 25 , . 5 1 , . 65 , 1 . 63 , 


83, 

a  55  y 

a  55  y 

.  75, 

.74, 

.  95 

57, 

a  63, 

a  55  y 

.  64, 

.  68, 

.  55 

77, 

07 

a  X—  sJ1  <y 

.  13, 

.27, 

.27, 

.  71 

29, 

nc r 

a  , 

.  32, 

n  t 
•  ^  ->  j 

.  36, 

.  36 

24, 

.  19, 

.28, 

.24, 

.  09, 

.  1  3 
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3250= 

3260= 

3270= 

3280= 

3290= 

3300= 

3310= 


3350= 
3360= 
3370= 
3380= 
3390= 
3400= 
34 1 0= 
3420= 
3430= 
3440= 
3450= 
3460= 
3470= 
3480= 
3490= 
3500= 
3510= 
3520= 
3530= 
3540= 
3550= 
3560= 
3570= 
3580= 
3590= 
3600= 
3610= 
3620= 
3630= 
3640= 
3650= 
3660= 
3670= 
3680= 
3690= 
3700= 
3710= 
3720= 
3730= 
3740= 
3750= 
3760= 
3770= 
3780= 


6  . 99 ,.89 ,1.16,1. 16, .73, .81, .21,1. 04 , .  2 1 , .  08 ,  1 . 

7  .38,1.29,1.11,1.16,1.59,. 20 , . 59 , . 99 , 1 . 33 , . 53 , 

8  2.19,1. 37, . 31 , . 31 , . 60, . 60, . 60, . 63, . 63, . 79, . 79 

9  .27, . 17, .29,2.65, .41/ 

DATA  FRRTS/.8, .45, .45, .5, .65, .45, .5, .85, .50, .55, .5 

1  .75, .8, .65, .8, .5, .85, .5, .45, .45, .6, .85, .8, .75 

2  .65, .7, .55, .7, .55, .7, .7, .75, .85, .45, .6, .5, .8, 

3  . 85, . 75, . 75, . 75, .75, .5, .6, .8, .8, .8, .8, . 65, 3*. 

4  . 75 , . 85 , . 8 , . 75 , . 55 , . 5 , . 6 , . 55 , .8, .7, .5, .6, .6, . 

5  .  45,  .  45,  .  55,  .  60,  .  75,  .  75,  .  8,  .  8,  .  8,  .  7,  .  8,  .  7,  .  7, 

6  . 65, . 55, . 55, . 5, . 55, .7, .7, . 65, .6, .8, .8, .6, . 65, 

7  . 80, .8, .8, .8, .8, .8, .7, .7, .7, 2*. 8, 2*. 85/ 

DATA  FRRTOK/ . 1 , . 45, . 45, . 40, . 25, . 45, . 4, . 05, . 4, . 35, . 

1  . 1 5 , . 1 , . 25 , . 1 , . 4 , . 05 , . 4 , . 45 , . 45 , . 3 , . 05 , . 1 , . 1 5 , . 1 , 

2  .  25,  .  2,  .  35,  .  2,  .  35,  .  2,  .  2,  .  15,  .  05,  .  45,  .  3 ,  .  4 ,  .  1 ,  .  1  ,  .  05 

3  .15,.15,.15,.15,.4,.3,.l,.l,.l,.l,. 25 , 3* . 3 , . 15, . 05 , 

4  . 1 , . 1 5 , . 35 , . 4 , . 3 , . 35 , . 1 , . 2 , . 4 , . 3 , . 3 , . 1 , . 45 , . 4 

c r  t  j  c  -jcr  4  4  h  o  1  '")  cr  -7  izr  ~scr 

W  I  -2*  m  X  U  4  •  A  vJ  jL.  f  ■  J.  4  ■  X—  f  f  ■  2  wJ  f  at  •  V_j  q  m  vJ  a|  a 

/  i-y  t-y  ^CT  ~.r  1  4  ^  -I  H  -4  -J  4 

C)  ■  x.  |  ■  xL  |  ■  x.  2  q  »  'J1  q  ■  X  <1  ■  1  ,  ■  2*  <1  «  X-  > — I 

7  . 2, . 2, . 1 , . 1 , . 05 , . 05/ 

DATA  FRNRTS/ 109*0. 10/ 

DATA  TEARDN/2*. 20, 2*. 27, .23, . 17, .27, . 17, .23, .23, .2 
1  .23, .23, .27, .20, .27, . 17, . 17, .23, .27, .23, .27, . 

07  07  h»  *r  07  r>r>  07 

*•  m  X-  /  <1  m  X—  /  q  ■  X  q  ■  X-  q  >  X<J  j|  ■  X  /  q  a  X  /  q  a  *.  /  q  a  ||  a  XV  j  a  X  /  q  ■ 

3  .  1 7 ,  .  1 7  ?  .  27 5 . 27 ?  .  20 ,  .  20 .  20 ,  .17,  .  20 ,  .  1 7 ,  .  20 ,  - 

41  7  1  ~7  r:»  t  r>T  17  17  07  r>7  o'T  o~7 

CT  ^  r-y~r  *  y  y  f-y-y  y-j  4  -y  1  ~7  4  ~7  “7  T  '7 

m  X  4  it  X  ■  4|  a  X  /  q  a  X  a|  ■  X  /  q  a  X  /fa  X  faJ./faJl  /  f  a  J.  /  f  a  f  ■  •— *  fa*' 

6  . 23, . 23, . 23, 4*. 23, 3*. 27, . 23, .23,7*. 17,2*. 20, . 

7  . 17, . 17, .23, . 17/ 

DATA  REF  A I R/ 2* . 37 , 2* . 50 ,  . 43 ,  . 30 , . 50 ,  . 30 ,.43,2*. 50 , 

1  . 50 , . 37 , .50,2*. 30 , . 42 , . 50 , . 43 , . 50 , . 43 ,  .  23 ,  .  43 

2  3* . 43 , 3* . 50 , . 43 , . 37 , . 50 , . 30 , . 50 , 2* . 30 , . 43 ,  . 50 

3  «  30,  •  37,  •  30,  .  37,  •  37,  .  30 ,  .  30,  2*.  30,  .  37,  .  43,  •  43 

4  . 30 , 2* . 50 , . 43, . 50 , . 43 , 2* . 37 , . 43 , . 30 ,.37,2*. 50 

5  5* «  3 0 , 6* • 37 , . 43 ,  • 30 , 5* • 43 , 3* . 50 , 2* . 43 , 7* . 30 , 3 

6  2*. 30, . 43, . 30/ 

NLRU  =109 
DSEED  =  123457. DO 

SUM  =  0 . 0 
DO  300  I  =  1 , NLRU 

SUM  =  SUM  +  <  QPA ( I ) /XMTBD ( I )  ) 

300  CONTINUE 

DEMFH  =  1.0/ SUM 

CUMF’RB  ( 1  )  =  QPA  (  1  )  *DEMFH/ XMTBD  <  1  ) 

DO  310  I  =  2, NLRU 
J  =  I  -  1 

FROB  =  QF'A(I)*DEMEH/XMTBD(I) 

CUMF'RB  (  I  >  =  CUMF'RB  ( J  )  +  F'ROB 
310  CONTINUE 

NTYF'E  =  4 
FAILRT  =  0.0 


r> 

p> 


3790= 
3800= 

38 1 0= 
3820= 
3830= 
3840= 
3850= 
3860= 
3870= 
3880= 
3890= 
3900= 
3910= 
3920= 
3930=*E0R 


o 


330 


DO  320  1  =  1,  IMTVPE 

FAILRT  =  FAILRT  +  < ATEQTY < I ) / TSMTBF ( I ) ) 
CONTINUE 

ATEREL  =  l.O/FAILRT 

ATEPRB ( 1 )  =  ATEQTY ( 1 ) *ATEREL/ TSMTBF ( 1 ) 

DO  330  I  =  2, NTYPE 
J  =  I  -  1 

PROS  =  ATEQTY ( I ) *ATEREL/ TSMTBF  (  I ) 

ATEPRB ( I )  =  ATEPRB ( J )  +  PROB 
CONTINUE 
RETURN 
END 


3940=GEN, ROARK, MODELS, 10, 
3950=S0U, 1,0, 1* 

3960=ACT, 1,1,00,120. 0* 
3970=ACT , 1 , 2 , CO ,  O .  O* 
3980=ACT ,1,3,00,24. O* 
3990=ACT ,1,4, CO, 48. 0* 
4000=ACT , 1 , 5 , CO , 72 . O* 

40 1 0= ACT , 1,6,00,96.  0* 
4020=REG , 2 , 1,1* 

4030= VAS, 2,1,00,1.0* 
4040=REG, 3, 1 , 1* 

4050= V AS , 3 , 1 , CO , 2 . O* 
4O60=REG, 4, 1,1* 

4O70=VAS , 4 , 1 , CO , 3 . O* 
4080=REG, 5, 1,1* 

4090= VAS , 5 , 1 , CO , 4 . 0* 

4 1 00=REG , 6 , 1 , 1 * 

41 10= VAS, 6, 1, CO, 5.0* 

4 1 20= ACT ,2,7* 


4130=ACT, 3, 7* 

4 1 40=ACT ,4,7* 

4 1 50=ACT ,5,7* 

4 160= ACT, 6, 7* 

4 1 70=REG , 7 , 1, 1* 

4 180= ACT, 7,8, CO, 11.0* 

4 1 90=REG , 8 , 1 , 1 , D , M* 
4200=ACT ,8,9, UF ,  1,1* 
4210=REG, 9, 1, 1,D,M* 
4220=ACT , 9 , 9 , UF , 1 , 2* 
4230= ACT, 9, 14* 
4240=REG, 10,1,1* 
4250=M0D, 20,9, 10* 
4260=M0D, 19,10,9* 
4270=ACT , 7 , 1 1 , CO , 25 . 0* 
4280=REG ,  1  1  , 1 ,  1 , D, M* 
4290=ACT , 1 1 , 1 2 , UF , 2 , 3* 
4300=REG, 12, 1 , 1 , D, M* 

43 1 0=ACT ,12,12, UF ,2,4* 
4320=ACT , 12,14* 


1983,8, , , 


6240. O, 5, E, O, 8* 

START  DAILY  ARRIVAL  PROCESS 

RE-START  EVERY  FIVE  DAYS 

MONDAY 

TUESDAY 

WEDNESDAY 

THURSDAY 

FRIDAY 

MONDAY 


TUESDAY 


WEDNESDAY 


THURSDAY 


FRIDAY 


MERGE  DAILY  ARRIVAL  STARTS 
ELEVEN  AM 

START  AM  POISSON  PROCESS 

CONTINUE  AM  POISSON  PROCESS 
ROUTE  TO  OMS  QUEUE 
DEAD  END  TO  STOP  PROCESS 
TO  SHUT  OFF  PROCESS 
TO  RE-START  PROCESS 
ONE  AM 

START  PM  POISSON  PROCESS 

CONTINUE  PM  POISSON  PROCESS 
ROUTE  TO  OMS  QUEUE 
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4330=F:EG ,  13,1,1* 

4340=M0D, 16,12,13* 

4350=M0D, 15,13,12* 

4360=QUE,  14/FLTL.INE,  (10)  15* 
4370=RES, 1 /QMS, 38, 15* 

4 380= ALL, 15, FOR, 1,2,14/16* 
4390=REG , 16,1,1* 

E 

4400=ACT, 16, 17, CO, 1 . O, 5/ DEBRIEF* 
44 1 0=REG ,  17,1,1* 

4420=ACT, 17, 18, ER, 1,6/REMOVAL* 
4430=PAR , 1 , 0. 5, , , 4* 

4440=REG , 18,1,1* 

4450=ACT, 18, 19, CO, O. 5, 7/PUTAWAY* 
4460=FRE , 1 9 , D , 1 , 2 , 1 5* 

4470= ACT, 19, 23, CO, 0.5* 

E 

4480=S0U , 20 , 0 , 1 , D , M* 

4490=ACT, 20, 20, US, 1 , 8* 

4500= ACT , 20 ,21* 

45 1 0=REG ,21, 1,  1* 

4520=REG ,22,1,1* 


DEAD  END  TO  STOP  PROCESS 

TO  SHUT  OFF  PROCESS 

TO  RE-START  PROCESS 

QUEUE  FOR  QMS 

OMS  TECHNICIANS 

ALLOCATE  TWO  OMS  TECHNICIANS 

START  OMS  (FLIGHT  LINE)  MAINTENA 

DEBRIEF 

REMOVE  LRU 

MEAN  REMOVAL  TIME  IS  TWO  HOURS 

TAKE  LRU  TO  SHOP 

FREE  TWO  OMS  TECHNICIANS 

LOGIN  AND  SCHEDULE  SHOP  MAINTENA 

START  ATE  FAILURES 

TIME  BETWEEN  ATE  FAILURES 

FAILURES  FROM  EIGHT  AM  TO  MIDNIG 

NO  FAILURES  FROM  MIDNIGHT  TO  EIG 


4530= 
4540= 
4550= 
4560= 
4570= 
4580= 
4590= 
4600= 
46 1 0= 
4620= 
4630= 
4640= 
4650= 
4660= 
4670= 
4680= 
4690= 
4700= 
4710= 
4720= 
4730= 
4740= 
4750= 
4760= 
4770= 
4780= 
4790= 
4800= 
4810= 
4820= 


MOD, 
MOD, 
ACT, 
REG, 
ACT, 
ACT, 
ACT, 
ACT, 
QUE, 
RES, 
ALL, 
QUE , 
RES, 
ALL, 
QUE, 
RES, 
ALL, 
QUE, 
RES, 
ALL, 
REG, 
ACT, 
REG, 
ACT, 
ACT, 
ACT , 
ACT, 
ACT, 
QUE , 
RES, 


21,21,22* 

17,22,21* 

21 . 23,  CO, 0. 5* 

23, 1 , 1 , A* 

23.24,  (9) Al.EQ. 1* 


23,25, (9) Al.EQ. 


16, (9) Al.EQ. 3* 

23,27, ( 9 ) A 1 . EQ . 4* 
24/DIGITAL, (6) S/4, (10)28* 
2/DIGITAL. 1 , 66, 28* 

28,  POR, 2, 1,24/32* 
25/DIGAV, (6) S/4, (10)29* 
3/DIGAV, 1,67,29* 

29,  POR, 3, 1,25/32* 

26/RF, (6) S/4, (10)30* 

4/RF , 1 , 68,30* 

30,  POR, 4, 1,26/32* 
27/RADAREW, (6) S/4, (10)31* 
5/RADAREW , 1,69,31* 

31,  POR, 5, 1,27/32* 

,1,1* 

32, 33, CO, 0.0, 9* 

3,  1,  1,A* 


34, 

<9 

)  A2. 

EQ. 

1* 

35, 

(9 

)  A2 . 

EQ. 

2* 

36 , 

(9 

)  A2. 

EQ. 

3* 

37, 

(9 

)  A2 . 

EQ. 

4* 

38, 

<9 

)  A2. 

EQ. 

5* 

TO  SHUT  OFF  FAILURES 
TO  OPEN  UP  FAILURES 
LOGIN  AND  SCHEDULE  AT  SHOP 
SELECT  STATION  TYPE 
ROUTE  TO  CORRECT  QUEUE 


QUEUE  FOR  AVAILABLE  STATION 
ONE  DIGITAL  STATION 
ALLOCATE  ONE  STATION 


ONE  DIGITAL  ANALOG  VIDEO  STATION 


ONE  RF  STATION 


ONE  RADAR  EW  STATION 


DO  NOTHING 

SELECT  TECHNICIAN  TYPE 
ROUTE  TO  CORRECT  QUEUE 


34/COMMNAV, (10)39* 
6/COMMNAV, 7, 77, 39* 


QUEUE  FOR  AVAILABLE  TECHNICIAN 
SEVEN  COMM  NAV  TECHNICIANS 


ALLOCATE  ONE  TECHNICIAN 


4830=ALL, 39, FOR, 6, 1 , 34/44* 
4840=QUE, 35/ AFC I NS , < 10)40* 
4850=RES, 7/ AFC INS , 4 , 78,  40* 
iNS 

4860=ALL, 40, FOR, 7, 1 , 35/44* 
4870=QUE, 36/OFFAV, (10)41* 
4880=RES, 8/OFF AV, 16,79,41* 
4890=ALL ,41, FOR , 8 , 1,36/44* 
4900=QUE , 37/DEFAV , (10) 42* 

49 1 0=RES , 9/DEF AV ,10, 80 , 42* 

» 

4920= ALL, 42, FOR, 9, 1 , 37/44* 
4930=QUE, 38/PMEL, ( 10)43* 
4940=RES, 10/PMEL ,5,81, 43* 
4950=ALL, 43, FOR, 10, 1,38/44* 
4960=REG , 44 , 1,1* 

4970=ACT , 44 , 45 , AT , 6 , 1 O* 
4980=REG, 45, 1 , 1 , A* 

4990= ACT, 45, 47, (9)A8.EQ.2* 

5000= ACT, 45, 46, (9)A8.EQ. 1* 
5010=REG, 46, 1 , 1 , P* 

5020=ACT ,46,47, CO ,1.0, ( 8 ) O . 75* 
5030= ACT , 46 , 47 , CO , O . 25 , ( 8 ) 0 . 25* 


FOUR  AUTO  FLT  CNTLS  INST  TECHNIC 

SIXTEEN  OFF  AVIONICS  TECHNICIANS 

TEN  DEFENSIVE  AVIONICS  TECHNICIA 

FIVE  PMEL  TECHNICIANS 

SETUP  AND  FAULT  ISOLATION 

SRU  SPARE  NOT  NEEDED 

GET  SRU  SPARE  FOR  LRU  REPAIR 

GET  SRU  SPARE  FROM  BASE  SUPPLY- 
GET  SRU  FROM  FORWARD  SUPPLY  POIN 


5040=REG, 47, 1,1* 

5050=ACT, 47, 48* 

5060=ACT ,45,149, CO ,1.0, ( 9 ) AS . EO . O* 
5070=REG, 149, 1 , i , A* 

5080= ACT ,149,1 50 , E  X , 2 , 1 1 ,  (9) A4. EQ. 3* 
5090=PAR, 2,0. 25* 

5100=REG, 150, 1,1* 

5 1 1 0= VAS ,  1 50 , 6 , E  X , 2* 

5 120= ACT, 150,48* 

5 1 30= ACT , 149, 151,  <9>A4.EQ. 1* 

5 1 40=FRE , 151, D, 10, 1,43,81* 

5 1 50=ACT , 151, 152* 

5 1 60=REG , 152,1,1* 

5 1 70= VAS , 152, 6, EX , 2,8, CO, 1 . O* 

5 1 80= ACT ,152, 33 , CO ,192. O* 

5190=REG, 48, 1 , 1 , A* 

5200=ACT, 48, 49, (9IA2.EQ. 1* 

52 1 0=ACT , 48 , 50 , (9) A2.EQ. 2* 

5220=ACT ,48,51, ( 9 ) A2 . EQ . 3* 

5230=ACT ,48,52, ( 9 ) A2 . EQ . 4* 

5240=ACT, 48,53, (9) A2. EQ. 5* 

5250=FRE , 49 , D , 6 , 1 , 77 , 39* 

5260=ACT , 49 , 54* 

5270=FRE, 50, D, 7,1,78, 40* 

5280=ACT, 50,54* 

5290=FRE , 51, D, 8, 1,79,41* 

5300= ACT ,51, 54* 

53 1 0=FRE , 52 , D , 9 , 1 , 80 ,42* 

532  0=ACT  ■  52  n  54* 


NO  SPARE  AVAILABLE  ON  BASE- 

TEAR  DOWN  LRU 

FIFTEEN  MINUTE  TEAR  DOWN 

CHANGE  A  T  6  TO  SETUP 

LRU  WILL  GO  AWAP 

STATION  GOES  NMCS 

FREE  TECHNICIAN  BUT  NOT  STATION 


CHANGE  AT6  TO  SETUP 

EIGHT  DAY  ORDER  AND  SHIP  TIME 


FREE 

COMM 

NAV  TECHNICIAN 

FREE 

AF  C 1 NS 

FECHN 1 C 1 AN 

F  REE 

OFF 

AV 

TECHN I C I  AN 

FREE 

DEF 

AV 

TECHNICIAN 

5340=ACT, 53, 54* 
5350=REG , 54 , 1 , 1 , A* 
5360= ACT, 54, 55, (9) 
5370=FRE, 55, D, 2,1, 
5380=ACT, 55, 59* 
5390=ACT, 54, 56, (9) 
5400=FRE, 56, D, 3, 1 , 


,  55 ,  ( 9 ) A 1 . EQ .  1  * 
,  D,2,  1,66,28* 
,59* 

,56, ( 9 ) A 1 . EQ . 2* 
,  D,3,  1,67,29* 


5410= 

5420= 

5430= 

5440= 

5450= 

5460= 

5470= 

5480= 

5490= 

5500= 

5510= 


5540= 

5550= 

5560= 


5590= 
5600= 
5610= 
5620= 
5630= 
5640= 
5650= 
5660= 
5670= 
5680= 
5690= 
5700= 
57 1 0= 


5740= 
5750= 
5760= 
5770= 
5780= 
5790= 
5800= 
58 1 0= 
5820= 
5830= 
5840= 


AC  T  , 
ACT, 
FRE, 
ACT, 
ACT, 
FRE, 
ACT, 
REG, 
ACT, 
REG, 
VAS, 
ACT, 
ACT, 
REG, 
ACT, 
ACT, 
ACT, 
ACT, 
QUE, 
QUE, 
QUE, 
QUE, 
ALL, 
ALL, 
ALL, 
ALL, 
REG, 
ACT , 
REG, 
ACT, 
ACT, 
ACT, 
ACT, 
ACT, 
QUE, 
QUE, 
QUE, 
QUE, 
QUE, 
ALL, 
ALL, 
ALL, 
ALL, 
ALL, 
REG, 
ACT, 


FREE  DIGITAL  STATION 


FREE  DIGITAL  ANALOG  VIDEO  ST AT  10 


56,59* 

54.57,  (9) A1.EQ.3* 

57,  D, 4, 1,68,30* 

57,59* 

54.58,  (9) A1.EQ.4* 

58,  D, 5, 1,69,31* 

58,59* 

59. 1 . 1 ,  A* 

59.61,  (9) A8.EQ.0* 

61. 1 ,  1  * 

61 , 8, CO, 1 . O* 

61, 23, CO, 192.0* 

59,60, (9) A8.NE.0* 

60.1.1,  A* 

60.62,  (9) Al.EQ. 1* 

60.63,  (9) Al.EQ. 2* 

60.64,  (9) Al.EQ. 3* 

60.65,  (9) Al.EQ. 4* 
62/DIGITAL,  (6)L.,  (10)66* 
63/DIGAV, (6) L, (10)67* 
64/RF , ( 6 ) L , (10)68* 
65/REW, ( 6 ) L , (10)69* 

66,  FOR, 2, 1,62/70* 

67,  FOR, 3, 1,63/70* 

68,  FOR, 4, 1,64/70* 

69,  FOR, 5, 1,65/70* 

70, 1,1* 

70, 71, CO, 0.0, 12* 

71 . 1 . 1 ,  A* 


FREE  RF  STATION 


FREE  RADAR  EW  STATION 


AWAP  LRU 
ORDER  SRU  SPARE 

EIGHT  DAY  ORDER  AND  SHIP  TIME 

ON  TO  NEXT  REPAIR  PHASE 
ROUTE  TO  CORRECT  QUEUES 


WAIT  FOR  AVAILABLE  STATION 


ALLOCATE  AVAILABLE  STATION 


71.72,  (9) A: 

71.73,  (9) A; 

71.74,  (9) AS 

71 . 75,  (9) Ai 

71.76,  ( 9 ) A; 


,  EQ.  1* 
,  EQ. 2* 
,  EQ . 3* 
,  EQ. 4* 
,  EQ . 5* 


DO  NOTHING 


ROUTE  TO  CORRECT  QUEUES 


72/COMMNAV,  <6)L,  (10)7/* 
73/ AFC I NS , ( 6 ) L , (10)78* 
74/OFFAV, ( 6 ) L , (10) 79* 
75/DEFAV, ( 6 ) L , (10)80* 
76/PMEL, ( 6 ) L , (10)81* 

77,  FOR, 6, 1 , 72/82* 

78,  FOR, 7,1, 73/82* 

79,  FOR, 8, 1,74/82* 

80,  FOR, 9, 1,75/82* 

81,  FOR, 10, 1,76/82* 

82, 1,1* 

82, 83, AT, 7, 13* 


WAIT  FOR  AVAILABLE  TECHNICIAN 


ALLOCATE  AVAILABLE  TECHNICIAN 


REPAIR  AND  TEAR  DOWN 


'■.v.'-'vy.y 


ROUTE  TO  CORRECT  FREE  NODES 


5870= 
5880= 
5890= 
5900= 
5910= 
5920= 
5930= 
5940= 
5950= 
5960= 
5970= 
5980= 
5990= 
6000= 
6010= 
6020= 
6030= 
6040= 
6050= 
6060= 
6070= 
6080= 
6090= 
6100= 
61  10= 
6120= 
6130= 
6140= 
6150= 
6 1 60= 
6170= 


REG, 

ACT, 

ACT, 

ACT , 

ACT, 

FRE, 

ACT, 

FRE, 

ACT, 

FRE, 

ACT, 

FRE, 

ACT 

REG 

ACT, 

ACT, 

ACT, 

ACT, 

ACT, 

FRE, 

ACT 

FRE 

ACT 

FRE 

ACT 

FRE 

ACT 

FRE 

ACT 

REG 

ACT 


83. 1.1,  A* 

83.84,  (9) Al.EQ. 1* 

83.85,  (9) A1 . EO. 2* 

83.86,  (9) Al.EQ. 3* 

83.87,  <9>  A1 . EQ. 4* 

84,  D, 2, 1 , 28, 66* 
84,88* 

85,  D, 3, 1,29,67* 
85,88* 

86 ,  D , 4 , 1 , 30 , 68* 

86 , 88* 

87,  D, 5, 1,31,69* 
87,88* 

88. 1. 1, A* 

88.89,  <9) A2.EQ. 1* 

88.90,  (9) A2.EQ.2* 

88.91,  (9) A2.EQ.3* 

88.92,  ( 9 ) A2 . EQ . 4* 

88.93,  (9) A2. EQ. 5* 

89,  D, 6, 1,39,77* 
89,94* 

90,  D, 7, 1 , 40, 78* 
90,94* 

91,  D, 8, 1,41,79* 
91,94* 

92,  D, 9, 1,42,80* 
92, 94* 

93,  D, 10, 1,43,81* 
93, 94* 

94,1,1* 

94, 95, CO, 0.5* 


FREE  STATIONS 


ROUTE  TO  CORRECT  FREE  NODES 


FREE  TECHNICIANS 


REPAIR  COMPLETED 

RETURN  LRU  TO  BASE  SUPPLY  &  LOGO 


6 1 80=REG , 95 , 1 , 1 , F* 

6 1 90= ACT , 95 , 96 , (8) 1,A5.EQ.3* 
6200= ACT ,95,97, ( 8 ) 2 , A4 . EQ . 3* 
62 1 0= ACT , 95 , 98 , (8>3,A4.EQ. 1* 
6220=REG,96, 1,1* 

6230— REG, 97, 1,1, A* 

6240=ACT , 97 , 99 , (9) Al.EQ. 1* 
6250=STA, 99/DIGITAL, 1, 1,D, I* 
6260=ACT , 97 , 100, (9) Al.EQ. 2* 
6270=STA, 100/DIGAV, 1, 1,D, I* 
6280=ACT ,97,101,  <  9 ) A 1 . EQ . 3* 
6290=STA , 101 /RF, 1, 1,D, I* 
6300=ACT , 97 , 102, (9) Al.EQ. 4* 
6310=STA, 102/REW, 1, 1,D, I* 
6320=REG , 98 , 1 , 1 , A* 

6330=ACT ,98,1 03 ,  <  9 ) A 1 . EQ .  1  * 
6340=STA, 103 /DIGITAL, 1 , 1 , D, I* 
6350= ACT, 98, 104, (9) Al.EQ. 2* 
6360= ST A, 104/DIGAV, 1 , 1 , D, I* 
6370= ACT, 98, 105, (9) Al.EQ. 3* 
6380=ST A , 105/RF, 1 , 1 , D, I* 
6390=ACT ,98,1 06 ,  ( 9 ) A 1 . EQ . 4* 


CONDITION  TAKE  FIRST 
REMOVE  NRTS  LRUS  FROM  STATS 
ALL  OTHER  LRUS  (RTS,  R'TOK ) 

ATE  MAINTENANCE 

COLLECT  STATS  ON  LRUS  BY  STATION 


COLLECT  STATS  ON  ATE  MAINTENANCE 


6400= 

6410= 

FT 

6420= 
6430= 
GHT 
6440= 
6450= 
6460= 
6470= 
6480= 
6490= 
6500= 
65 1 0= 
6520= 
6530= 
6540= 
6550= 
6560= 
6570= 
6580 
6590 
6600 
66 1 0 
6620 
6630 
6640 
6650 
6660 
6670 
6680 
6690 
6700 
6710 
6720 
6730 
6740 
6750 
6760 
6770 
6780 
6790 
6800 
6810 
6820 
6830 
6840 
6850 
6860 
6870 
6880 
6890 
6900 
6910 


:STA ,  106/REW,  1 ,  1  ,  D,  1* 
SOU, 1 10,0, 1* 

:ACT ,  110,  1  11* 

:ALT ,  1  1  1 ,  D,  1 ,  -38,  15* 


ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
REG, 
ACT, 
REG, 
ACT, 
REG, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT , 
ALT, 
ACT, 
ALT, 
ACT, 
ALT, 
ACT, 
REG, 
ACT, 
REG, 


111, 

112, 

112, 

113, 

113, 

114, 

114, 

115, 

115, 

116, 
116, 
117, 

117, 

118, 
1  18, 
119, 

119, 

120, 
120, 
121 , 
121, 

L  x-x-  , I 
<  nn 

A  i 

123, 

123, 

124, 

124, 

125, 

125, 

126, 
126, 
127, 

127, 

128, 
128, 
129, 

129, 

1 30 , 

130, 

131, 

131, 

132, 

132, 

133, 

133, 

134, 

134, 

135, 


1 12* 

D,  2, — 1 , 66, 28* 
113* 

D,  3,  -1 , 67,  29* 
114* 

D , 4 , - 1 , 68 , 30* 

1  15* 

D ,  5 , - 1 , 69 , 3 1  * 
116* 

D, 6, -7,77, 39* 
117* 

D, 7,-4,78,40* 

1  18* 

D,  8,  -16, 79,41* 
119* 

D,  9, -10, 80, 42* 
120* 

D, 10, -5,81 , 43* 

121 ,  CO, 0.0, 21* 

1,1* 

122,  CO, 1.0, 15* 

1,1* 

1 23 ,  CO ,4.0, 1 6* 

1,1* 

124,  CO, 3.0, 17* 
D, 1 , 38, 15* 

125* 

D,6,  7,77,39* 
126* 

0,7,4,78, 40* 

12  7* 

D, 8, 16, 79,41* 
128* 

D, 9, 10,80,42* 
129* 

D, 10,5,81, 43* 
130, CO, 0.25, 18* 
D, 2, 1,66,28* 
131* 

D, 3, 1,67,29* 
132* 

D, 4, 1,68,30* 
133* 

D, 5, 1,69,31* 

134,  CO, 2. 75, 19* 
1,1* 

135,  CO, 4. 0,20* 

1,1* 


v.v-.v.’.v '.v  ■ 


BEGIN  CLOCK  AT  END  OF  SECOND  SHI 

WHICH  IS  AT  MIDNIGHT 

SHUT  DOWN  ALL  RESOURCES  AT  M1DNI 


MIDNIGHT  —  SEE  NODES  21,22 

ONE  AM  —  SEE  NODES  12,13 

FIVE  AM  —  SEE  NODES  12,13 

EIGHT  AM  —  SEE  NODES  21,22 
START  FIRST  SHIFT  AT  EIGHT  AM 


FIFTEEN  MINUTE  DELAY  FOR  STATION 
DAILY  CONFIDENCE  CHECK 
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ELEVEN  AM  —  SEE  NODES  9,10 
THREE  F’M  —  SEE  NODES  9,10 


.*  .*  /  *,•  '  •  * , 


6920= ACT , 1 35 , 1 1 0 , CO , 9 . O* 
6930=FIN* 

694Q=*EGR 


AND  BACK  TO  MI DN I GH1 
THAT'S  ALL,  FOLKS' 


Appendix  B:  Design  of  Simulation  Experiment 

Purpose  of  the  Experiment 

The  simulation  experiment  has  just  one  response 
(dependent.)  variable,  which  is  base  repair  cycle  time.  The 
factor  to  be  varied  is  test  station  quantities;  this  factor 
has  three  levels  (quantities  of  1,  2,  or  3  test  stations). 
The  purpose  of  the  experiment  is  to  determine  if  varying 
test  station  quantities  will  cause  a  statistically 
significant  difference  in  repair  cycle  time.  This 
experiment  is  performed  a  total  of  12  times,  once  for  each 
combination  of  the  three  aircraft  quantities  (16,  26,  or  32 
F'AA)  and  the  four  station  types  (Digital,  Digital  Analog 
Video,  RF,  or  R/EW) .  The  aircraft  quantities  and  station 
types  were  not  regarded  as  factors  since  a  separate 
management,  decision  (how  many  stations  to  buy)  must  be  made 
for  each  of  the  12  combinations. 


VisCLSDE®  Reduction 

To  make  the  simulation  experiment  as  powerful  as 
possible  (for  a  cjiven  sample  size),  it  is  desirable  to 
remove  as  much  variation  (caused  by  the  simulation  model)  as 
possible.  This  can  be  accomplished  by  the  variance 
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Through  use 


reduction  technique  of  common  random  numbers, 
of  a  FORTRAN  SAVE  statement,  it  was  possible  to  ensure  that 
each  level  (treatment)  -faced  the  same  random  numbers  which 
determined  LRU  inter arrival  times,  LRU  identification,  LRU 
type  of  maintenance  action  (RTS,  RTOK,  or  NRTS) ,  and  LRU 
maintenance  task  times.  The  FORTRAN  user  function  UF,  as 
modified  to  ensure  common  random  numbers,  is  attached  at  the 
end  of  this  appendix.  Since  the  test  station  quantities 
were  varied,  it  was  not  possible  to  ensure  common  random 
numbers  for  test  station  failure  interarrival  times  or 
repair  times.  However,  model  results  indicate  that  at  least 
95"/.  of  the  random  numbers  are  common. 

The  method  of  common  random  numbers  would  be  useless 
unless  the  simulation  experiment  can  isolate  variation  due 
to  different  random  number  streams  from  variation  due  to 
random  error.  This  isolation  can  be  accomplished  by 
blocking,  where  each  block  corresponds  to  a  random  number 
stream.  The  linear  model  for  such  an  experiment  would  then 
be: 


=  JK  + 

where 

=  dependent,  variable  resulting  from  block, 

j+H  level 

Z4 

=  over a 1 1  mean 

pi 

=  cfW 

block  effect 

% 

level  effect 

-  random  error  term  tor  block,  level 

V-  V  V 


fa 


In  a  blocking  experiment,  the  blocking  effect  is 


assumed  to  be  significant.  In  addition,  it  is  also  assumed 


that  there  is  no  block-level  interaction  of  effect, 


A NOVA  Layout 


In  this  experimental  design,  there  are  ten  replications 


performed,  which  is  the  maximum  permitted  in  a  single  Q-GERT 


run.  bach  replication  has  a  different  stream  of  random 


numbers  and  thus  corresponds  to  a  block  in  the  simulation 


experiment.  The  level  or  treatment  corresponds,  to  the  test 


station  quantities  considered.  The  layout  for  such  an 


experiment  would  be: 


BLOCK 


QTY~  1 


FACTOR 
QTY  =2 


QTY=3 


The  degrees  of  freedom  for  such  an  experiment  would  be 


as  foil ows : 


Treatment  or  Level  —  3-1  =  2 

Block  —  10-1  =  9 

Error  —  <3—1 >  x  (10-1)  =  18 

Total  —  30—1  =  29 

The  null  nypothesis  to  be  tested  is  that  varying  test 
station  quantities  does  not  cause  a  statistically 
significant  change  in  base  repair  cycle  time.  In 
statistical  terms,  this  is  represented  by: 

Ho  •  =  7i  *  2i  *  o 

ANOVA  Assumptions 

As  part  of  the  standard  (normal  distribution)  ANOVA,  i 
is  necessary  to  assume  that  the  error  terms,  ,  have  zero 

mean,  have  constant  variance,  are  independent,  and  follow  a 
normal  probability  distribution.  When  this  experiment  was 
performed,  it.  soon  became  quite  obvious  that  the  variance 
for  different  levels  was  not  constant.  It  was  therefore 
necessary  to  resort,  to  Friedman's  test,  which  is  the 
nonpar  ametr  i  c:  equivalent  of  the  one-way  ANOVA  with  blocking 
The  only  disadvantage  to  a  nonpar  ametr i  c  test,  is  that  it 
might  not  be  as  powerful  as  an  ANOVA.  As  will  be  seen,  thi 
is  not  a  concern  for  this  particular  situation. 


E'iQ§r.i!D®Qtalu  Results 

The  results  -for  the  12  simulation  experiments  are  shown 
in  Table  XVIII.  Significance  levels  for  each  of  the  12 
combinations  were  obtained  from  Friedman’s  test. 

Table  XVIII.  Significance  Levels 


Digital 

Digital 

An  Vi^deo 

RF 

R/EW 

16  PAA 

.  045 

.  000 

.  000 

.  000 

26  PAA 

.  002 

.  000 

.  000 

.  000 

32  PAA 

.  000 

.  000 

.  000 

.  000 

The  conclusion  that  can  be  drawn  is  that  varying  test 
station  quantities  causes  a  statistically  significant 
difference  in  base  repair  cycle  time  for  each  of  the  12 
combinations.  It  is  also  possible  to  perform  3  pairwise 
comparisons  for  each  of  the  12  experiments.  For  o*.  =  0.05, 
32  out  of  36  pairwise  comparisons  were  found  to  be 
statistically  significant.  For  the  remainder  of  this 
research  effort,  therefore,  it  was  assumed  that  varying  test 
station  quantities  (over  the  range  of  one  to  three  stations) 
always  caused  a  statistically  significant  difference  in  base 
repair  cycle  time. 

For  comparison  purposes,  a  one-way  ANOVA  with  blocking 
was  performed  for  each  station  type  for  the  16  PAA  case-.  In 
each  case,  the  null  hypothesis  was  rejected  for  »«.  =  0.5. 


L  I  ST ,  1 ooy 1 500 

1 00=LMR , CM 1 77000 , T600 , 10200 .  T830662 , ROARK , 4533 . 

1 10=ATTACH, IMSL, IMSL, I D=L I BRAR Y , SN= ASD . 

1 20=L I BR ARY , IMSL. 

1 30=ATT ACH , PROOF  I L , QGERTPROC ,  I D= AF I T . 

140=FTN5, ANSI=0. 

1 50=BEG I N , QGERT , PROCF I L , M=LGO , MODE=  X . 

160=*E0R 
1  70= 

180=  FUNCTION  UF(IFN) 

190= 

200=  COMMQN/QVAR/NDE , NFTBU ( 500) ,NREL(500) ,NRELP(500) , 

210=  1  NREL2  <  500 )  ,  NRUN ,  NRIJNS ,  NTC  ( 500 )  ,  PAR  AM  (  1 00 ,  4 )  ,  TBEG ,  TNOW 

220= 

230=  COMMON  /  UCOM 1  /DEMFH ,  FLYF’ER,  ACBASE,  FHACMO, 

240=  1  FACTDY (5)  , FACTTM (2) 

250= 

260=  COMMON /UC0M2/CUMPRB (200) , STNREQ (200) , AFSC (200) , FRRTS (200)  , 

270=  1  FRRTOK ( 200 ) , FRNRTS ( 200 ) , SETUP ( 200 ) , TSS I TA ( 200 ) , 

280=  2  WARMUP ( 200 ) , PERFT ( 200 ) , REPA I R ( 200 ) , TEARDN ( 200 ) 

290= 

300=  C0MMQN/UC0M3/DSEED, SRUAVL 

310= 

320=  DOUBLE  PRECISION  DSEED, XSEED 

330=  SAVE  XSEED 

340=  DSEED  =  DSEED  +  ( 100000. *NRUN> 

350=  IF (TNOW  .LE.  11.001)  XSEED  =  754321. DO  +  100000. *NRUN 

360= 

370=  ATI  =  GATRB(l) 

380=  I ATI  =  ATI 

390= 

400=  DEMSH  =  (DEMFH*FLYPER*4. 333) / (FACTDY (I ATI)* 

410=  1  FACTTM ( IFN) *ACBASE*FHACMO> 

420= 

430=  UN IF  =  GGUBFS( XSEED) 

440=  UF  =  — 1 . 0*DEMSH*AL0G ( UN I F ) 

450=  IF (UF  .GT.  16.0)  UF  =  16.001 

460= 

470=  UN IF  =  GGUBFS( XSEED) 

480=  I  =  1 

490=  IF (UN IF  .LT.  CUMPRB ( I ) )  GO  TO  110 

500=  100  1=1+1 

510=  J  =  I  -  1 

520=  I F ( ( CUMPRB ( J ) . LE . UN IF)  .AND.  (UNIF. LT. CUMPRB ( I )) )  GO  TO  110 

530=  GO  TO  100 

540=  110  CONTINUE 

550= 

560=  A 1  =  STNREQ ( I ) 

570=  CALL  PATRB ( A 1 , 1 ) 

580= 

590=  A2  =  AFSC ( I ) 

600=  CALL  PATRB <A2, 2) 


CALL  PATRB (A4, 4) 


630= 
640= 
650= 
660= 
670= 
680= 
690= 
700= 
7 1 0= 
720= 
730= 
740= 
750= 
760= 
770= 
780= 
790= 
800= 
810= 
820= 
830= 
840= 
850= 
860= 
870= 
880= 
890= 
900= 
910= 
920= 
930= 
940= 
950= 
960= 
970= 
980= 
990= 
1 000= 
1010= 
1 020= 
1 030= 
1 040= 
1 050= 
1 060= 
1070= 
1 080= 
1 090= 
1 1 00= 
1110= 
1  1 20= 
1 1 30= 
1  1 40= 
1 1 50= 
1  160= 


UN IF  =  GGUBFS ( XSEED) 

XLOW  =  FRRTS(I) 

XMID  =  FRRTS  ( I )  +  FRRTQK ( I ) 

XHIGH  =  XMID  +  FRNRTS ( I ) 

IF  <UNIF  .LT.  XLOW)  GO  TO  120 
I F  <  UN  I F  .LT.  XMID)  GO  TO  130 
IF (UN IF  .LE.  XHIGH)  GO  TO  140 

1 20  A5  =  1 . O 

CALL  PATRB <A5, 5) 

SUM1  =  SETUP  ( I  )  +  TSGITA  (  I  )  +  WARMUP  (  I )  +  O.  5*F'ERFT  (  I  ) 
SUM2  =  (0. 75) *REPAIR  < I )  +  (O. 25) *2. 0*REPAIR < I ) 

1  +  F’ERFT  (  I )  +  TEARDN  (  1  ) 

SUM  =  SUM 1  +  SUM2 

ERLANG  =  0.0 
DO  125  K  =  1,2 

UN IF  =  GGUBFS (XSEED) 

DRWEXP  =  -1 . 0* (SUM/2. O) *ALOG (UNIF) 

ERLANG  =  ERLANG  +  DRWEXP 
125  CONTINUE 

A6  =  ( SUM 1 /SUM ) *ERLANG 
CALL  PATRB (A6, 6) 

A 7  =  (SUM2/SUM) *ERLANG 
CALL  PATRB (A7, 7) 

UNIF  =  GGUBFS (DSEED) 

IF (UNIF  .LT.  SRUAVL )  A8  =  1.0 
IF (UNIF  .GE.  SRUAVL)  A8  =  0.0 
CALL  PATRB (AS, 8) 

RETURN 


130  A5  =  2.0 

CALL  PATRB (A5, 5) 

SUM1  =  SETUP ( I )  +  TSSITA(I)  +  WARMUP ( I ) 
SUM2  =  0.3*PERFT(I>  +  TEARDN ( I ) 

SUM  =  SUM1  +  SUM2 


PERFT ( I ) 


ERLANG  =  0.0 
DO  135  K  =  1,5 

UNIF  =  GGUBFS ( XSEED) 

DRWEXP  =  -1.0* (SUM/5. 0)*ALOG (UNIF) 
ERLANG  =  ERLANG  +  DRWEXP 
CONTINUE 

A6  =  (SUM1 /SUM) *ERLANG 
CALL  PATRB (A6, 6) 
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A 7  =  (SUM2/SUM) *ERLANG 
CALL  PATRB  < A7, 7) 


AS  =  2.0 

CALL  PATRB (A8, 8) 

RETURN 

140  A5  =  3.0 

CALL  PATRB (AS, 5) 

SUM  1  =  SETUP ( 1 )  +  TSSITA < I )  +  WARMUP ( I )  + 
SUM2  =  TEARDN ( I ) 

SUM  =  SUM 1  +  SUM'2 

ERLANG  =  0.0 
DO  145  K  =  1,5 

UN IF  =  GGUBFS  <  XSEED) 

DRWEXP  =  — 1 . O* (SUM/ 5. O) *ALOG (UNIF) 
ERLANG  =  ERLANG  +  DRWEXP 
145  CONTINUE 

AS  =  (SUM1/SUM) *ERLANG 
CALL  PATRB (AS, 6) 

A 7  =  (SUM2/SUM) *ERLANG 
CALL  PATRB (A7, 7) 

A8  =  2.0 
CALL  PATRB (A8, 8) 

RETURN 
END 


1  '30 


0. 5*PERFT (  I  ) 


Appendix  C 


Analytical  Queueing  Model  Program  Listing 


140 


1 00=LMR , CM  1 25000 , T 1 00 , 1 050 .  T830662 , ROARK , 4533 . 

1 1 0=ATT ACH , I MSL , I MSL5, I D=L I BRARY , SN= ASD . 

120=LIBRARY, I MSL. 

1 30=FTN5. 

1 40=LG0 . 

150=*EOR 

160=  PROGRAM  QUEUE4 

170= 

180=  COMMON/UCQM/ XMTBD ( 109) ,QPA(109> , STNREQ < 109) , SETUP (109) 

1 90=  1  TSS I TA ( 1 09 ) , WARMUP ( 1 09 ) , PERFT ( 1 09 ) , FRRTS ( 1 09 ) , 

200=  2  FRRTOK  <  1 09 )  ,  FRNRTS  <  1 09 )  ,  TEARDN  <  1 09 )  ,  REF'A  I R  (  1 09 )  , 

210=  3  ATEREL (4) , ATEQTY (4) , ATEREP (4) 

220= 

230=C  COMMON  BLOCK  UCOM  IS  USED  TO  TRANSFER  DATA  VALUES 

240=C  FROM  BLOCK  DATA  SUBROUTINE 

250= 

260=  DIMENSION  BMHAVG  < 4 ) , DEMFH ( 4 ) , DEMSH ( 4 ) 

270= 

280=C  NUMBER  OF  AIRCRAFT  FOR  THIS  BASE 

290=  DATA  ACBASE/ 16. 0/ 

300= 

3 1 0=C  FLIGHT  HOURS  PER  AIRCRAFT  PER  MONTH 

320=  DATA  FHACMO/29.2/ 

330= 

340=C  NUMBER  OF  MAINTENANCE  HOURS  PER  WORKDAY  —  TWO  SHIFTS 

350=  DATA  HRSDAY/ 16.0/ 

360= 

370=C  NUMBER  OF  DAYS  PER  WEEK  —  NO  MAINTENANCE  ON  WEEKENDS 

380=  DATA  DAYWK/5. 0/ 

390= 

400=C  4.333  WEEKS  PER  MONTH 

410=  DATA  WEEKMO/4. 333/ 

420= 

430=C  HOURS  PER  MAINTENANCE  DEBRIEFING 

440=  DATA  DBRIEF/ 1 . 0/ 

450= 

460=C  HOURS  TO  REMOVE  LRU  FROM  AIRCRAFT 

470=  DATA  REMOVE/2. 0/ 

480= 

490=C  HOURS  TO  TAKE  LRU  TO  AMS 

500=  DATA  PUT AW Y/ 0.5/ 

510= 

520=C  19  PAIRS  OF  OMS  TECHNICIANS 

530=  DATA  NUMOMS/ 19/ 

540= 

550=  DATA  BMHAVG/ 4*0. 0/ 

560=  DATA  RTSRT/O.O/ 

570=  DATA  SRUAVL/O. 9/ 

580= 

590=C 
600= 


TRAVEL  TIME  TO  OBTAIN  SRU  SPARE 
TRAVEL  =  (0.25) *0.25  +  (0.75) *1.0 


620=C 

630= 

640= 

650=C 

660= 

670= 

680=C 

690= 

700= 

710= 

720= 

730= 

740= 

750= 

760= 

770= 

780= 

790= 

800=C 

810= 

820= 

830=C 

840= 

850= 

860= 

870= 

880= 

890= 

900-C 

910= 

920= 

930= 

940= 

950=C 

960= 

970= 

980=C 

990= 

1 000= 
1010=C 
1020= 
1030= 
1040=C 
1050= 
1060= 
1070= 
1080=C 
1090= 

1 100= 

1 1 10=C 
1 120= 

1 1 30= 

1 140=C 


EIGHT  DAY  ORDER  AND  SHIP  TIME 
HRSQST  =  8.0*16.0 


109  LRUS 
NLRU  =  109 

COMPUTE  TOTAL  LRU  ARRIVAL  RATE  BY  STATION  TYPE 

DO  110  J  =  1,4 

REALJ  =  J 
SUM  =  0.0 
DO  100  I  =  1 , NLRU 

IF (STNREQ ( I )  .EQ.  REALJ)  THEN 
SUM  =  SUM  +  QPA(I) /XMTBD(I) 

END  IF 
100  CONTINUE 

DEMFH  IS  AGGREGATE  MTBD  FOR  JTH  STATION  TYPE  LRUS 
DEMFH(J)  =  1.0/SUM 

CONVERT  FLYING  HOURS  TO  MAINTENANCE  SHIFT  HOURS 
DEMSH(J)  =  DEMFH < J ) * (HRSDAY*DAYWK*WEEKMO> / 

1  <  ACBASE*FHACMQ ) 

110  CONTINUE 

SUM  =  0.0 

COMPUTE  TOTAL  LRU  ARRIVAL  RATE  —  ALL  STATIONS  COMBINED 
DO  120  J  =  1,4 

SUM  =  SUM  +  ( 1 . O/DEMSH  < J ) ) 

120  CONTINUE 

COMBINED  AVIONICS  RELIABILITY 
XMTBF  =  1.0/ SUM 

O-LEVEL  REMOVAL  TIME  (MEAN) 

XMTTR  =  DBRIEF  +  REMOVE  +  PUTAWY 

SUBROUTINE  FOR  M/M/S  QUEUE 

CALL  MMS( XMTBF, XMTTR, NUMOMS, RHOOMS, WOMS) 

TOTAL  TIME  FOR  OMS  QUEUE;  CONVERT  16  HOUR  DAY  TO  REAL  TIME 
WOMS  =  (24. O/HRSDAY) *WOMS 
PRINT  *,  •'  WOMS  =  ’  ,  WOMS 

FOR  ALL  LRUS  —  COMPUTE  WEIGHTED  AVEREAGES 
DO  130  I  =  1 , NLRU 

7.  REPAIR  THIS  STATION 

RTSRT  =  RTSRT  +  (QF’A  (  I )  *  XMTBF  /  XMTBD  (  I  )  )  *FRRTS  ( I  ) 


MEAN  REPAIR  TIME  FOR  RTS 


1150=  BMH1  =  SETUP ( I )  +  TSSITA(I)  +  WARMUP ( I ) 

1160=  1  +  0.5*PERFT(I)  +  TRAVEL  +0. 75*REPAIR ( I ) 

1170=  2  +  0.25*2.0*REF'AIR(I)  +  PERFT ( I )  +  TEARDN ( I ) 

1180= 

1 190=C  MEAN  TASK  TIME  FOR  RTOK 

1200=  BMH2  =  SETUP  < I >  +  TSSITA(I)  +  WARMUP ( I ) 

1210=  1  +  PERFT ( I )  +  0.3*PERFT(I>  +  TEARDN < I ) 

1220= 

1 230=C  MEAN  TASK  TIME  FOR  NRTS 

1240=  BMH3  =  SETUP ( I >  +  TSSITA(I)  +  WARMUP < I ) 

1250=  1  +  O. 5*PERFT ( I )  +  TEARDN < I ) 

1260= 

1270=C  MEAN  TASK  TIME  OVERALL 

1280=  BMHLRU  =  FRRTS ( I > *BMH1  +  FRRTOK ( I ) *BMH2  +  FRNRTS ( I ) *BMH3 

1290= 

1 300=C  PARTITION  AVERAGES  BY  STATION  TYPE 

1310=  DO  140  J  =  1,4 

1320=  REALJ  =  J 

1330=  IF (STNREQ  < I >  .EQ.  REALJ)  THEN 

1340=  HOLD  =  <  QPA ( I ) *DEMFH  <  J )  / XMTBD ( I ) ) *BMHLRU 

1350=  BMHAVG(J)  =  BMHAVG(J)  +  HOLD 

1360=  END  IF 

1370=  140  CONTINUE 

1380=  130  CONTINUE 

1390= 

1 400=C  POST  —  7.  OF  DEMANDS  WHICH  CAUSE  LRU  AWAP 

1410=  POST  =  RTSRT* < 1 . 0— SRUAVL) 

1420= 

1430=C  AVERAGE  TIME  FOR  ORDER  AND  SHIP  TIME 

1440=  WOST  =  P0ST*8. 0*24.0 

1450=  PRINT  WOST  =  ’,WOST 

1460= 

1470=  DO  160  J  =  1,4 

1480=  IF (J  .EQ.  1)  PRINT  ’ DIGITAL  STATION’ 

1490=  IF (J  .EQ.  2)  PRINT  *, ’DIGITAL  ANALOG  VIDEO  STATION’ 

1500=  IF ( J  .EQ.  3)  PRINT  *, ’RADIO  FREQUENCY  STATION’ 

1510=  IF (J  .EQ.  4)  PRINT  RADAR  EW  STATION’ 

1520=  PRINT  *, ’  DEMSH  =  ’.DEMSHCJ) 

1530=  PRINT  *, ’  BMHAVG  =  ’ , BMHAVG ( J ) 

1540= 

1550=C  MEAN  TIME  BETWEEN  STATION  NMCS 

1560=  ATEREL(J)  =  ATEREL ( J ) / ( 1 . O  -  SRUAVL) 

1570= 

1 580=C  MEAN  ORDER  AND  SHIP  TIME 

1590=  ATEREP(J)  =  HRSOST 


1600= 

1610=C 

EITHER  ONE  <K  = 

1  )  OR 

TWO  <K  =  2)  STATIONS 

1620= 

1 630= 

DO  150  K  =  1,2 

1 640= 

IF (K  .EQ.  1) 

THEN 

1650=C 

SUBROUTINE 

FOR  ONE 

STATION  CASE 

1 660= 

CALL  ME21 (DEMSH (J) 

, ATEREL ( J )  , BMHAVG ( J )  , 

1670= 

1  ATEREP(J), 

RHO . W ) 

1680= 

END  IF 
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1700= 

1 710=C 

1720= 

1730= 

1740= 

1750= 

1 760=C 
1 770=C 
1 780=C 
1 790=C 
1 800= 
1810= 
1820= 

1 830= 
1840= 
1850= 


I F ( K  . EQ ■  2 )  THEN 

subroutine  for  two  station  case 

CALL  ME 2 2 <  DEMSH ( J )  , ATEREL ( J ) , BMHAVG ( J ) , 
ATEREP(J) , RHO, W) 

END  IF 


RHO  IS  STATION  UTILIZATION 

W  IS  TOTAL  TIME  IN  STATION  QUEUE  +  SERVICE  TIME 
W  IS  CONVERTED  FROM  16  HOUR  TO  24  HOUR  DAYS 
WTOT  IS  TOTAL  BASE  REPAIR  CYCLE  TIME 


1950= 
1960= 
1970= 
1 980= 
1 990= 
2000= 
2010= 
2020= 
2030= 
2040= 
2050= 
2060= 
2070= 
2080= 
2090= 
2100= 
2110= 
2 1 20= 
2 1 30= 
2140= 
2150= 
2160= 
2170= 
2180= 
2190= 


or? 


00= 


2210= 

2220= 


W  =  (24. O/HRSDAY) *W 
WTOT  =  WOMS  +  W  +  WOST 


1 860= 

I F  ( RHO 

.  LT. 

1.0) 

THEN 

1870= 

PRINT 

*,  ’ 

K  = 

’  ,  K ,  ' 

RHO  =  ’ , RHO 

1 880= 

PRINT 

*5  ' 

W  = 

7  ,W,  ’ 

WTOT  =  ? , WTOT 

1890= 

END  IF 

1 900= 

1910= 

IF  (RHO 

.GE. 

1.0) 

THEN 

1920= 

PRINT 

*,  ’ 

K  = 

7  ,K,  ’ 

RHO  =  •’  ,  RHO 

1 930= 

END  IF 

1940=  150 

CONTINUE 

160  CONTINUE 
STOP 
END 


SUBROUTINE  MMS ( XMTBF , XMTTR, NUMSRV, RHO, W> 
REAL  NUMER , LQ 


FR  =  1.0 /XMTBF 
SR  =  1.0/ XMTTR 
S  =  NUMSRV 


RHO  =  (FR)/(S*SR) 


SUM  =  0.0 

LIMIT  =  NUMSRV  -  1 


DO  400  N  =  O, LIMIT 


XN  =  N 

DENOM  =  FACT ( XN) 
NUMER  =  (FR/SR) **N 
HOLD  =  NUMER /DENOM 
SUM  =  SUM  +  HOLD 
400  CONTINUE 


NUMER 

DENOM 


(FR/SR) **S 
FACT (S) 


DENOM 


DENOM* ( 1 . O-RHO ) 


2230= 

2240= 

2250=  HOLD  =  NUMER/DENOM 

2260= 

2270=  SUM  =  SUM  +  HOLD 

2280= 

2290=  F’O  =  1.0  /SUM 

2300= 

2310=  NUMER  =  PO*NUMER*RHO 

2320=  DENOM  =  DENOM* < 1 . O-RHO) 

2330= 

2340=  LQ  =  NUMER/DENOM 

2350= 

2360=  WQ  =  LQ/FR 

2370= 

2380=  W  =  WQ  +  <1.0 /SR) 

2390= 

2400=  RETURN 

2410=  END 

2420= 

2430=  FUNCT I  ON  FAC  T' t  X  ■ 

2440=  FACT  =  1.0 

2450=  IX  =  X  +  0.5 

2460=  IF  < I X  .GT.  1)  THEN 

2470=  DO  500  I  =  2, IX 

2480=  FACT  =  FACT* I 

2490=  500  CONTINUE 

2500=  END  IF 

2510= 

2520=  RETURN 

2530= 

2540=  END 

2550= 

2560=  SUBROUTINE  ME2 1 < XMTBF 1 , XMTBF2, XMTTR 1 , XMTTR2, RHO, WAIT ) 

2570= 

2580=  D I MENS ION  A  < 1 00 , 1 00 ) , B  < 1 00 , 1 ) , WK ARE A < 1 00 ) 

2590= 

2600=  DO  290  I  =  1,100 

2610=  DO  280  J  =  1,100 

2620= 

2630=  A  <  I ,  J  )  =0.0 

2640= 

2650=  280  CONTINUE 

2660= 

2670=  B  < 1 , 1 )  =0.0 

2680= 

2690=  290  CONTINUE 

2700= 

2710=  M  =  1 

2720=  NMAX  =  8 

2730=  N  =  4*NMA X  +  2 

2740=  I A  =  100 

2750=  IDGT  =  3 

2760= 
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2770= 

FR1  =  1.0/XMTBF1 

2780= 

FR2  =  1.0/ XMTBF2 

2790= 

SRI  =  1 . 0/XMTTR1 

2800= 

2810= 

SR2  =  1 . 0/ XMTTR2 

2820= 

A <1,1)  =  -1.0* (FR1  +  FR2) 

2830= 

A (1,2)  =  SR2 

2840= 

A  < 1 , 3)  =  2. 0*SR1 

2850= 

A  <2, 1 )  =  FR2 

2860= 

A <2, 2)  =  -1.0* (FR1  +  SR2) 

2870= 

A  <3, 3)  =  — 1 . 0* (FR1  +  FR2 

+  2. 

2880= 

A  <3, 4)  =  SR2 

2890= 

A <3, 5)  =  2. 0*SR1 

2900= 

A <4, 3)  =  FR2 

2910= 

2920= 

A  <4, 4)  =  — 1 . O* (FR1  +  SR2) 

2930= 

2940= 

LIMIT  =  NMAX  -  2 

2950= 

DO  300  K  =  1, LIMIT 

2960= 

11=  4*K  +  1 

2970= 

12  =  4*K  +  2 

2980= 

13  =  4*K  +  3 

2990= 

3000= 

14  =  4*K  +  4 

3010= 

A< I 1 , I 1-4)  =  FR1 

3020= 

A<I1,  ID  =  —  1 . 0*  <  FR 1  + 

FR2 

3030= 

A(I1, I 1+1 )  =  SR2 

3040= 

Adi,  11+2)  =  2. 0*SR1 

3050= 

A  < 12, 12-4)  =  FR1 

3060= 

A (I 2 , 12-1)  =  FR2 

3070= 

A <12, 12)  =  — 1 . 0*  < FR1  + 

SR2) 

3080= 

A< 13, 13-4)  =  FR 1 

3090= 

A <13, 13)  =  — 1 . 0*  <FR1  + 

FR2 

3 1 00= 

A  < 13, 13+1)  =  SR2 

31 10= 

A <13, 13+2)  =  2. 0*SR1 

3120= 

A ( I 4 , 14-4)  =  FR1 

3130= 

A  < 1 4 , 1 4- 1 )  =  FR2 

3140= 

3 1 50= 

A< 14,  14)  =  — 1 . 0*  <FR1  + 

SR2) 

3 1 60= 
3170= 

300  CONTINUE 

3180= 

3190= 

K  =  LIMIT  +  1 

3200= 

11  =  4*K  +  1 

3210= 

12  =  4*K  +  2 

3220= 

13  =  4*K  +  3 

3230= 

3240= 

14  =  4*K  +  4 

3250= 

A  <  1 1 , 11-4)  =  FR1 

3260= 

A(I1,I1)  =  — 1 . 0*  <FR1  +  FR 

2  +  : 

3270= 

A  <  1 1 , 1 1  + 1 )  =  SR2 

3280= 

A <11, 11+2)  =  2. 0*SR1 

3290= 

A  < 12, 12-4)  =  FR1 

3300= 

A  <12, 12-1)  =  FR2 

0*SR1 ) 


+  2.0*SR1) 


+  2. 0*SR1 ) 


2 . 0*SR 1 ) 
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A  (12, 12)  =  -1.0* (FR1  +  SR2) 

A ( 13, 13-4)  =  FR1 

A (13, 13)  =  -1.0* (FR2  +  2 . 0*SR 1 ) 

A ( 1 3,  13+1)  =  SR2 

A (13, 13+2)  =  2. 0*SR 1 

A(  14,  14-4)  =  FR1 

A  ( 14,  14-1)  =  FR2 

A ( 14, 14)  =  -1. 0*SR2 

11  =  II  +  4 

12  =  12  +  4 

A  ( I  1 ,  I  1 -4 )  =  FR 1 
A  ( 1 1 ,  II)  =  -1.0* (FR2  +  2. 0*SR1 ) 
A  ( 1 1 , 1 1  + 1 )  =  SR2 

DO  310  K  =  1,12 

A  <  1 2 ,  K )  =  1.0 


510  CONTINUE 

B(I2, 1)  =  1.0 

CALL  LEQT1F ( A, M, N, IA, B, IDGT, WKAREA, IER) 

PRINT  *, '  ME21  IER  =  ’ , IER 

SRAVG  =  (FR1*SR1  +  FR2*SR2) / (FR1  +  FR2) 

RHO  =  (FR1  +  FR2) /SRAVG 

XL  =  0.0 

DO  320  K  =  1 , NMAX 

11  =  4* (K— 1 )  +  3 

12  =  4* (K— 1 )  +  4 

13  =  4* ( K— 1 )  +  5 

14  =  4* (K— 1 )  +  6 
REALK  =  K 

HOLD  =  REALK* (B(I1, 1)  +  B(I2,1)  +  B(I3,1)  +  B(I4,1)) 

XL  =  XL  +  HOLD 
520  CONTINUE 

WAIT  =  XL/FR1 

RETURN 


SUBROUTINE  ME22 ( XMTBF 1 , XMTBF2, XMTTR1 , XMTTR2, RHO, WAIT) 
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D I MENS  ION  A(1 00 ,  1 00 )  , B ( 1 00 ,  1 )  , WKAREA < 1 00 ) 

DO  395  I  =  1,100 
DO  385  J  =1,1 00 
A (I ,  J )  =  0.0 
85  CONTINUE 

B<  I  ,  1  )  =  0.0 
95  CONTINUE 

M  =  1 
NMAX  =  8 
N  =  9*NMAX 
I A  =  100 
I DGT  =  3 

FR1  =  1.0/XMTBF1 
FR2  =  1 . 0/ XMTBF2 
SRI  =  l.O/XMTTRl 
SR2  =  1.0/ XMTTR2 


Ad,  1) 

=r 

-1.0* (FRl 

+ 

2. 0*FR2 ) 

Ad, 2) 

= 

SR2 

A  <  1 , 4  ) 

= 

2 . 0*SR  1 

A (2,  1  ) 

= 

2. 0*FR2 

A  (2, 2) 

= 

-1.0* (FRl 

+ 

FR2  +  SR2 ) 

A  (2,  5) 

= 

2. 0*SR1 

A  ( 2 , 3 ) 

= 

2 . 0*SR2 

A  <3, 2) 

= 

FR2 

A  <3, 3) 

= 

-1.0* (FRl 

+ 

2. 0*SR2) 

A  (4,  4) 

-1.0* (FRl 

+ 

2. 0*FR2  +  2. 0*SR1 ) 

A  (4, 5) 

= 

SR2 

A ( 4 , 7) 

= 

2.  0*SR 1 

A  ( 4 ,  1 0 ) 

=  4 . 0*SR 1 

A  ( 5 , 4 ) 

= 

2. 0*FR2 

A  ( 5 , 5 ) 

= 

-1.0* (FRl 

+ 

FR2  +  2 . 0*SR 1  +  SR2) 

A  (5,  6) 

= 

2.  0*SR2 

A  <5, 8) 

= 

2. 0*SR 1 

A (5, 11) 

=  2. 0*SR1 

A (6, 5) 

=s 

FR2 

A (6, 6) 

= 

-1.0* (FRl 

+ 

2. 0*SR2> 

A (7,  1 ) 

= 

FRl 

A  <7, 7) 

= 

-1 . 0* (FRl 

+ 

2. 0*FR2  +  2. 0*SR 1 ) 

A  (7, 8> 

- 

SR2 

A ( 7 , 13) 

=  2. 0*SR1 

A  (8,  2) 

= 

FRl 

A  <  8 , 7 ) 

= 

2. 0*FR2 

A  ( 8,  8 > 

-1.0* (FRl 

4 

FR2  +  2. 0*SR1  +  SR2) 

A  <8, 9) 

= 

2. 0*SR2 

A (8, 14) 

=  2.  0*SR 1 

A  <  9 , 3 ) 

= 

FRl 

A  <9, 8) 

=r 

FR2 

A  (9, 9) 

= 

-1.0* (FRl 

•4 

2. 0*SR2 ) 

4390= 

A  <10,  10)  =  -l.OMFRl  + 

2 . 0*FR2  +  4 . 0*SR 1 ) 

4400= 

A  < 10,  1 1 )  =  SR2 

4410= 

A (10, 13)  =  2. 0*SR1 

4420= 

A(ll, 10)  =  2 . 0*FR2 

4430= 

A  < 1 1 ,  1 1 )  =  — 1 . O*  <  FR 1  + 

FR2  +  2.0*SR1  +  SR2 ) 

4440= 

Adi,  12)  =  2. 0*SR2 

4450= 

A <12, 11)  =  FR2 

4460= 

A ( 1 2 ,  12)  =  - 1 . 0*  <  FR 1  + 

2 . 0*SR2 ) 

4470= 

A (1 3 , 4 )  =  FR 1 

4480= 

A  d  3 , 1 3 )  =  - 1 . 0*  <  FR 1  + 

2. 0*FR2  +  4 . 0*SR 1 ) 

4490= 

A <13, 14)  =  SR2 

4500= 

A <13, 16)  =  4 . 0*SR 1 

4510= 

A <13, 19)  =  4. 0*SR1 

4520= 

A <14, 5)  =  FR1 

4530= 

A <14, 13)  =  2. 0*FR2 

4540= 

A <14,  14)  =  — 1 . 0*  <FR1  + 

FR2  +  2.  O+SR 1  +  SR 2! ) 

4550= 

A <14, 15)  =  2 . 0*SR2 

4560= 

A  d  4 ,  1 7 )  =  2 . 0*SR 1 

4570= 

A <14, 20)  =  2 . 0*SR 1 

4580= 

A  < 1 5 , 6 )  =  FR 1 

4590= 

A  < 1 5,  14)  =  FR2 

4600= 

4610= 

A (15,  15)  =  - 1 , 0* ( FR 1  + 

2 . 0*SR2 ) 

4620= 

4630= 

LIMIT  =  NMAX  -  2 

4640= 

4650= 

DO  400  K  =  2, LIMIT 

4660= 

11  =  9*<K-1)  +  7 

4670= 

12  =  11  +  1 

4680= 

13  =12+1 

4690= 

14  =  13  +  1 

4700= 

15  =  14  +  1 

4710= 

16  =  15  +  1 

4720= 

17  =  16  +  1 

4730= 

18  =  17  +  1 

4740= 

4750= 

19  =  18  +  1 

4760= 

A(I1, 11-9)  =  FR1 

4770= 

A<I1,  ID  =  —  1 . 0*  (FR1 

+  2.0+FR2  +  4 . 0*SR 1 ) 

4780= 

Adi  ,  12)  =  SR2 

4790= 

Adi,  17)  =  2. 0*SR1 

4800= 

A <12, 12-9)  =  FR1 

4810= 

A  <12,  ID  =  2.0+FR2 

4820= 

A  < 12, 12)  =  —  1 . 0*  < FR 1 

+  FR2  +  2 . 0*SR 1  +  SR2) 

4830= 

A <12, 13)  =  2. 0*SR2 

4840= 

A <12, 18)  =  2. 0*SR1 

4850= 

A <13, 13-9)  =  FR 1 

4860= 

A <13, 12)  =  FR2 

4870= 

A <13,  13)  =  —  1 . 0*  < FR 1 

+  2.0+SR2) 

4880= 

A  < 1 4 ,  14-9)  =  FR 1 

4890= 

A  (  I  4,  14)  =  —  1  .  O*  <FR  1 

+  2.0+FR2  +  4 . 0*SR 1 ) 

4900= 

A  <1 4 ,  1 5 )  =  SR 2 

4910= 

A (I  4,  17)  =  2-0*SRl 

4920= 

A  <15,  15-9)  =  FR 1 

4930= 

A  (  I  5 , 

14)  = 

2 .  O  *  1  R  2 

4940= 

A  (  15, 

15)  = 

- 1 . 0* ( F  R 1 

+ 

F  R  2  + 

4950= 

A  (  15, 

16)  = 

2. 0*SR2 

4960= 

A  (16, 

16-9) 

=  FR1 

4970= 

A  (  16, 

15)  = 

FR2 

4980= 

A  <16, 

16)  = 

- 1 . 0* ( F  R 1 

+ 

2. 0*SR 

4990= 

A  (  I  7  , 

17-9) 

=  FR1 

5000= 

A  <  I  7 , 

17)  = 

- 1 . 0* ( F  R 1 

+ 

2, 0*FR 

50 1 0= 

A  (  I  7 , 

18)  = 

SR2 

5020= 

A  (  I  7 , 

1  1+9) 

=  4 . 0*SR 1 

5030= 

A  <  I  7 , 

14  +  9) 

=  4.0*SR1 

5040= 

A  (  18, 

18-9) 

=  FR1 

5050= 

A  (  18, 

17)  = 

2. 0*FR2 

5060= 

A  <  18, 

18)  = 

-1 . 0* (FR1 

+ 

FR2  + 

5070= 

A  <  1 8 , 

19)  = 

2. 0*SR2 

5080= 

A  (18, 

12+9) 

=  2. 0*SR1 

5090= 

A  (  18, 

15+9) 

=  2. 0*SR 1 

5 1 00= 

A  (  19, 

19-9) 

=  FR1 

51  10= 

A(  19, 

18)  = 

FR2 

5120= 

A  (  19, 

19)  = 

-1 . 0*  <FR1 

+ 

2 . U*SR 

i’ .  U*SR  1 


5 1 30= 
5140= 
5150= 


0*SR 1 ) 


+  2. U*SR 1  +  SR2 ) 


400  CONTINUE 


5 1 60= 

11  =  11 

+  9 

5170= 

12  =  12 

+  9 

5180= 

13  =  13 

+  9 

5190= 

14  =  14 

+  9 

5200= 

15  =  15 

+  9 

5210= 

16  =  16 

+  9 

5220= 

17  =  17 

+  9 

5230= 

18  =  18 

+  9 

5240= 

19  =  19 

+  9 

5250= 

5260= 

Adi,  11- 

9)  =  FR1 

5270= 

A(U,  ID 

=  —  1 . 0*  (FR'l 

5280= 

A ( I 1 , 12) 

=  SR2 

5290= 

A ( I 1 , 17 ) 

=  2.0* SRI 

5300= 

A  <  1 2 ,  1 2- 

9)  =  FR'l 

53 1 0= 

A ( 12, 11) 

=  2. 0*FR2 

5320= 

A ( 12, 12) 

=  - 1 . 0* ( FR 1 

5330= 

A< 12, 13) 

=  2. 0*SR2 

5340= 

A< 12, 18) 

=  2 . 0*SR 1 

5350= 

A< 13, 13- 

9 )  =  FR 1 

5360= 

Ad  3  ,  12) 

=  FR2 

5370= 

A ( 13, 13) 

=  — 1 . 0*  <  FR 1 

5380= 

A ( I 4 , 14- 

9)  =  FR 1 

5390= 

A( 14, 14) 

=  -1.0* (2.0 

5400= 

A ( 14, 15) 

=  SR2 

5410= 

A( 14, 17) 

=  2. 0*SR1 

5420= 

A (15, 15- 

9)  =  FR'l 

5430= 

A( 15, 14) 

=  2. 0*FR2 

5440= 

A ( 15, 15) 

=  — 1 . 0* (FR2 

5450= 

A ( 15, 16) 

=  2. 0*SR2 

5470= 

A (16, 15)  =  FR2 

5480= 

A <16, 16)  =  -1 . 0* <2. 0*SR2) 

5490= 

A  a 7,  17-9)  =  FR1 

5500= 

A ( I 7 , 1 7 )  =  - 1 . 0* ( 2 . 0*FR2  +  4 . 0*SR 1 ) 

55 1 0= 

A (17, 18)  =  SR2 

5520= 

A (17,11+9)  =  4.0*SR1 

5530= 

A ( 18, 18-9)  =  FR1 

5540= 

A (18, 17)  =  2. 0*FR2 

5550= 

A (18, 18)  =  -1.0* (FR2  +  2. 0*SR 1  +  SR2 ) 

5560= 

A (18, 19)  =  2 . 0*SR2 

5570= 

A (18, 12+9)  =  2 . 0*SR 1 

5580= 

A( 19, 19-9)  =  FR1 

5590= 

A ( 19, 18)  =  FR2 

5600= 
56 1 0= 

A ( I 9 , 1 9 )  =  - 1 . 0* ( 2 . 0*SR2 ) 

5620= 

11  =  11  +9 

5630= 

12  =  12  +  9 

5640= 

13  =  13  +  9 

5650= 

A  ( I  1 ,  I 1-9)  =  FR1 

5660= 

A  < I 1 ,  1 1 )  =  - 1 . 0* ( 2 . 0*FR2  +  4 . 0*SR 1 ) 

5670= 

A  (  1 1 ,  1 2 )  =  SR2 

5680= 

A (12, 12-9)  =  FR1 

5690= 

A(I2,  ID  =  2 . 0*FR2 

5700= 

A  (12,  12)  =  -1.0* (FR2  +  2. 0*SR 1  +  SR2 ) 

5710= 

5720= 

A (12, 13)  =  2. 0*SR2 

5730= 

DO  410  K  =  1,N 

5740= 

A  ( 1 3 ,  K )  =  1.0 

5750= 

410 

CONTINUE 

5760= 

5770= 

B  < 13, 1 )  =  1.0 

5780= 

5790= 

CALL  LFQT1F (A, M, N, IA,B, IDST, WKAREA, IER) 

5800= 

5810= 

PRINT  *,’  ME 22  IER  =  '  ,  IER 

5820= 

SRAVG  =  ( FR 1 *SR 1  +  FR2*SR2) / (FR1  +  FR2) 

5830= 

5840= 

RHO  =  <  FR 1  +  FR2) / (2. 0*SRAVG> 

5850= 

XL  =  0.0 

5860= 

XL  =  B  ( 4 ,  1 )  +  B  (5,  1  )  +  B  ( 6 ,  1  ) 

5870= 

5880= 

XL  =  XL  +  B  ( 7,  1 )  +  B  (8,  1  )  +  B(9,l) 

5890= 

DO  420  K  =  2 , NMAX 

5900= 

REALK  =  K 

59 1 0= 

11  =  9  *  <  K  —  1 )  +  1 

5920= 

SUM  =  B ( I 1 , 1 )  +  B ( I 1 + 1 , 1 )  +  B ( I 1 +2 , 1 ) 

5930= 

SUM  =  SUM  +  B( I 1+3,1)  +  B< I 1+4,1)  +  B( 11+5,1) 

5940= 

SUM  =  SUM  +  B< 1 1+6,1)  +  B< I 1+7,1)  +  B(  I 1+8,1) 

5950= 

5960= 

XL  =  XL  +  SUM*REALK 

5970= 

5980= 

420 

CONTINUE 

5990= 

6000= 

WAIT  =  XL/FR1 

RETURN 

END 


60  1  0= 
6020= 
6030= 
6040= 
6050= 
6060= 
6070= 

6080= 
6090= 
6 1 00= 

61 10= 

6120= 
6 1 30= 

6140= 

6 1 50= 
6 1  60= 

6 1  70= 
6 1 80= 

6 1 90= 
6200= 
62 1 0= 
6220= 
6230= 
6240= 
6250= 
6260= 
6270= 
6280= 
6290= 
6300= 
6310= 
6320= 
6330= 
6340= 
6350= 
6360= 
6370= 
6380= 
6390= 
6400= 
6410= 
6420= 
6430= 
6440= 
6450= 


BLOCK  DATA 

COMMON /LJCOM/  XMTBD  (  1 09 )  ,  ORA  (  1 09 )  ,  STNREQ  <  1 09 )  ,  SETUP  ( 1  09 )  , 

1  TSS I T A ( 1 09 )  , WARMUP ( 1 09 )  , PERFT < 1 09 )  , FRETS  < 1 09 )  , FEE  T  Ok  < 1 0 


FRNRT'S  ( 1 09 )  ,  TEARDN  (  1 09 )  ,  REP  A I R  (  1 09 ) 

, ATEREL <  4 )  , ATEQTY <  4  >  , ATEREP < 4 ) 

DATA  XMTBD/ 16458. , 1 199. , 707. , 1488. , 2024. , 1576. , 952. , 100O00 


1 923. , 11731. , 151.  , 560. , 590 . , 780 . , 2260 . 


j690  . 


10.  ,  101 


620. .  740. ,  1800. , 460. , 232.  , 272. , 3803. ,4218. , 424. ,  1  ISO. 

1923. .  2085. . 446. .2175. . 1931 . . 3653. . 3538. . 372. . 2019. . 3 

1 3600 . .  1 00000 . . 1 530 . .41. . 14223. . 8509 ..1391.. 2085 . . 86 1 


7 

8 

9 

A 

B 

1 


9667 . .  5856 . . 967 . . 5 1 86 . . 1 2 1 7 1 . . 5877 . . 6385 ..3176.. 5037 . 

2203. .  36000. . 3030. . 198. . 289. . 478. . 282. .418. . 1391. . 437 

245. .  1597.  ,  1597.  , 221 .  , 7838.  , 847.  ,  1380.  , 3030.  , 3030.  , 

159..  225 . . 208 . . 486 ..191.. 269 . . 225 . . 1 56 . . 1 09 . . 344 ..185 

157. .  223. . 136. .111.. 265. . 198. . 136. . 159. . 58. . 449. . 449. 

422. . 422. .422. .216. .216. . 175. . 175. .50.  ,2906. ,82. , 

706. . 4480. / 

DAT  A  QPA/3 1 .  , 3*2 . ,  1 . , 3*2 . , 5.  ,  1 .  , 3. , 8.  ,  1 .  , 3.  , 2.  ,  1 .  ,  1 .  , 3. , 


8*1 . , 2*1 . , 2. , 1 . , 2*2. , 1 . 
2*1. ,2. ,5*1. ,2. ,5*1. ,2. 
1 . , 2 . , 3 . , 3*3. ,7*1. , 2*3 . 
1.  ,  1.  ,2.  ,  1  .  ,2.  ,2.  ,  1.  ,2. 


3 . ,2*1. ,2. ,4*1. ,2. , 1. ,3. , 
7*1 . , 4*2. , 10. , 1 . , 2. , 1 . , 3 . , 
1  .  ,  2.  , 2*1  .  ,  3.  ,  1  .  , 2.  ,  1  .  , 

2.  ,  1 .  ,  1  .  / 


DATA  STNREQ/41*2. 0,34*1.0, 18*4.0, 16*3.0/ 

DA  TA  SET  UP/2*.  23,  2*.  33,  .  28,  .  18,  .  33,  .18,  2*.  28, 2*.  33,  3*.  28, 

. 33 ,.18,. 33 , 2* . 1 8 , . 28 , . 33 , . 28 , . 33 , . 28 , . 1 8 , . 28 , . 33 , 
3*. 28, 3*. 33, . 28, . 23, . 33, . 18, . 33,2*. 18,2*. 33,3*. 23, 

. 18, . 23, . 18, 2*. 23, 2*. 18, 2*. 18, . 23, . 28, . 28, . 18, . 18, 

. 33 , . 33 , . 28 , . 33 , . 28 , . 23 , . 23 , .28, . 18, . 23 , . 33 , -  33 , 

. 23 , 5* . 18, . 23 , 5* . 23 , . 28 , . 18,5*. 28, 3*. 33 , 2* . 28 , 

7*. 18,2*. 23, . 23,2*. 18, . 28, . 18/ 

DATA  TSS I TA/ 1 09*0 .17/ 

DATA  WARMUP / 1 09*0 . 0/ 

DATA  PERFT / . 20 ,.59,1.44,. 37 , 1 . 08 , .29, 1.31, . 16, . 75 ,.61,. 05 
. 83, . 55, . 55, . 75, . 74, . 95, . 95, . 89, 1 . 25, .51, . 65, 1 . 63, 


1 


wJ 

6) 


i 


4 

5 

6 

7 


. 57 , . 63 , . 55 , . 64 , . 68 , . 55 
. 77, . 23, . 13, . 27, . 27, . 71 
. 29 , . 25 , . 32 , . 23 , . 36 , . 36 
.24, . 19, . 28 , . 24 , . 09 , . 1 3 


1  .  15,  .97,  .31,  .  5  /  ,  .  5 5 .  6 5 , 

.  32 ,  .  12,  .25,  .  12,  1  .<  '9  ,  .  72 , 

.  56 ,  .  72 ,  . 28 ,  . 68 ,  . 2 7 ,1.1 6 , 

. 85, 1 . 35, . 43, .12,1. 7 1,1. O- 
. 99, . 89, 1.16,1.16,. 73, . 81 , . 31 , 1 . 04, . 21 , . 08, 1 . 52, 

. 38,  1 . 29,  1.11,  1.16,1. 59,  . 2u,  . 59,  . 99,  1  .  33,  . 53, .61, 


6430= 
4  9  O — 
o5o0= 
6510= 
65.20= 
6530= 
654 0= 
6550= 
6560= 
6570= 
6580= 


DATA 


1  HR  I  fcj / .8, . 4b , . 4b .  . 5 ,  . 6b ,  . 45 
.  75,  8,  .65,  .8,  .  5,  .  85,  .  5,  .  45 


nc  -yc,  ~/c r  "7  c  'tc  / 

03  \_J  «(  ■  /  ij  ||  •  /  1-1  n  ■  /  \-2  q  a  /  \_1  ^  n  u-J  a  a  O  «|  a 

75 ,  85 ,  .  8 ,  .  75 ,  .  55 ,  .  5 ,  .  6 ,  .5 

45,  ..  45,  .  55,  .  60,  .  75,  .  75,  .  B, 


DATA 


8  A ,  .  8 ,  .  8 ,  .  8 ,  .  8 ,  ..  8 ,  .  7  ,  .  7  ,  .  7 
RRTGK/ . 1 , . 45, .  45,  . 40, . 25, . 
15,-1,. 25 ,.1,.4,. 05 , . 4 ,  .  45 

,  -  ■■■■'  ■■■  I  ,  ■  ,  ■  jL.  ,  «  1  ... 


6590= 

-:r 

a  1  5  a  l  5  ,  a  1  5  j  a  15,  . 

4,  .  3,  .1,-1 

6600= 

4 

. 1 , . 15, .35, .4, .3, 

.  35,  -1,-2, 

66  1 0= 

cr 

x.J 

a  -3  (|  a  1  \J  ||  8  1.  «.J  1  ■  1  1  •  1  f 

.  1,  -2,  .  1  ,  . 

6620= 

.6 

■  j  ■*•*—*»  «»  -fc-  y  •  *  j  » 

»-•  n  ■ 

6630= 

y 

a  if  a  rf  a  In  M  1  f  a  U  x- j  ,  • 

05/ 

6640= 

DATA 

FRNR  1  8/  109*0.  1.0/ 

6650= 

DATA 

1  LARON/ 2*  a  20  ,  au’* . 

t  n  *  2-3  ,  a  1  /  , 

6660= 

1 

a  •.!«  y  ■  *!.  *.J»  ,  a  ail  j  ,  ••  a^.  O  y  a 

27,  .  17,  .  1  7 

667 O — 

•~i 

r  >  ■/  /-5  y  r  y  ' 

"  ■<-  /  ,  ■  a_  /  ,  a  .a_  •_*  ,  -  aU.  ,  » 

aU.  ^  a  .aL  /  ,  a  a6-  / 

6680= 

•:rt 

.  17, . 17, . 27 , -  27 ,  . 

20 , . 20 , . 20 

6690= 

4 

■j  “7  •->.*  i  7  :•  '  :>  “T 

■  J.  /  ,  ”  ■  *•  t  <i  ■  a’...  ,  •«  •A  •..  *  , 

. 23, -17.-1 

6700= 

5 

■  ,  a  a_-  '  H  ■  •*-  ,  •  a...  / 

H  m  XL  /  n  "  A  .  ,  a  1 

67  1 0= 

6 

-t  ?t  4-u.  '.r-r 

"I,*-  ■;>-/  ';y-r 

n  "  a  a-.  /  ,  "  ■»•••  *-• 

6720= 

7 

•  l/  y«17ya  -a-  -3  ,  .  1/  / 

6730= 

DATA 

REPAIR/ . 37, . 37, . 50, . 50, . 43, 

6740= 

1 

2*. 50, 3*. 43, . 50, . 

37 , .50, . 30 

6750= 

1;* 

. 43, . 50,  . 43, -  23, . 

43, . 50, 3*. 

6760= 

. 37, . 50, . 30, . 50, . 

30, . 30, . 43 

6770= 

4 

.  37, . 30, . 37, . 30, . 

37, . 37, 2*. 

6780= 

5 

.  43,  .  43, 2*. 30, 2*. 

50, . 43, . 50 

6790= 

6 

. 43, . 30, . 37, 2*. 50 

, . 37 , 5* . 30 

6800= 

7 

.  43 ,  .  30 ,  5* .  4 3 3* . 

50 , 2* . 43 , 7 

6810= 

8 

•  -3  7  ,  a  -3  J  ,  ■  -30  ,  a  4  -3  ,  . 

30/ 

6820= 

DATA 

ATEQTY/4*! . 0/ 

6830= 

DATA 

A  TEREP/4* 1 . 0/ 

6840= 

DA  1  A 

A  1  LKtiL./  .  ,  1  4trj  -  , 

167. , 130. / 

6850= 

686O=*E0R 

END 

Bi  b  1  :t  ography 


Morse,  Philip  li.  Queues^  inventories  and  Maintenance 
New  York:  John  Wiley  &  Sons,  Inc:.,  1958. 

Sauer,  Charles  H.  and  K.  Mam.  Chandy.  "Approximate 
Solution  of  Queueing  Models,  "  Computer:.  25-31  (Apr i 
1980) . 

Ross,  Sheldon  M.  introducti  on  to  Probabi  l  i  ty  Models.. 
New  York:  Academic  Press,  Inc.,  1980. 

Cooper,  Robert  B.  introduction  to  Queueing  Ihegry. 
New  York:  North  Holland,  1981. 

Kleinrock,  Leonard.  Queueing  Systems,  Vgiyme  1.1 
Theory,.  New  York:  John  Wiley  &  Sons,  Inc.,  1975. 

kleinrock,  Leonard.  Queueing  Systems,.  Volume  iii 
Computer  Applications,  New  York:  John  Wiley  &  Sons, 
Inc.,  1976. 

Basket,  Forest,  et  al .  "Open,  Closed  and  Mixed 
Networks  of  Queues  with  Different  Classes  of 
Customers,  "  Journal,  of  the  Association  of  Computing 
Machinery,.  22: 248-260  (April  1975). 

Cox,  Michael  H.  Myi.ti.-Cl.ass  Analytical  Models  of  the 
D E L s y s te nr-1 0  Job-Swapping  behavior,  MS  thesis. 

Wr i qh t —Patter son  AFB,  Ohio:  Air  Force  Institute  of 
Technology,  December  1981.  (G0R/MA/81D— 4)  . 

Shanthi kumar ,  J.  G.  "M/G/l  Queues  with  Scheduling 
within  Generations  and  Removable  Server,"  Gperatfons 
Resear'h,.  29:  1010—1018  (September-October  1981). 

Carroll ,  J.  L. ,  et  al .  "Solutions  of  M/G/l //N-type 
Loops  with  Extensions  to  M/G/l  and  GI/M/1  queues," 

Q P §ra 1 1  o n s  Research x  30:490-514  (May- June  1982). 

Neuts,  Marcel  F.  "Explicit  Steady-state  Solutions  to 
Some  Elementary  Queueing  Models,"  Operations 
Research,  30: 480-489  (May-June  1982) . 

Clark,  Thomas  D.  and  Mary  K.  Allen.  "A  Logistical 
Analysis  of  the  F— 111  Automatic  Test  Equipment:  A 

Network  Simulation  Approach,"  AFALD  Presentation, 
Wright— Patterson  Air  Force  Base,  Ohio. 


General  Dynamics,  F o r t  W orth  Division.  S  t  a t i_ o n 

LQ§.di,ng  Model  User's  Manual ^  General  Dynamics 
document  16PR1345  (12  June  1980). 

Bryson,  Gary  C.  ,  et.  al  .  A  Decision  Suggort  System  of 
Accjui  si  t  i  on  of  F-16  Avi  oni  cs  Intermedi  ate  Shop  Test 
Sets  Using  the  System  Science  Paradigm  and  Q—GERT^  MS 
Thesis.  Wr i ght-Patter son  APB,  Ohio:  Air  Force 
Institute  of  Technology,  1982.  (USSR  11-82) . 

Strategic  Air  Command.  Maintenance  Portion  of  the 
System  Operational  Concept  for  the  B-1B^  HQ  SAC 
Planning  Document.  Offutt  APB,  Nebraska:  HQ  SAC, 

1983. 

Deputy  for  B-1B.  B— IB  integrated  Logistics  Support 
Plan.  B-1B  SPO  Planning  Document.  Wr i ght-Patter son 
APB,  Ohio:  ASD/B1,  January  1983. 

Hermes,  Phillip  H.  SAC  Base  Operational  and 
Maintenance  Data.  Informal  Report..  Wright— Patterson 
APB,  Ohio*.  ASD/YYEE,  1983. 

Biffin,  Walter  C.  Queueing..  Columbus:  Grid,  Inc., 

1978. 

1MSL,  Inc.  IMSL  Library  Reference  Manual j_  Vgl_ume  2._ 
Houston:  IMSL,  Inc.  1982. 

Concepts  and  Analysis  Directorate,  Acquisition 
Logistics  Division,  Air  Force  Logistics  Command. 
Logistics  Suggort  Cost  Mgdei  Uteris  Handbook..  Version 
1.1.  AFALD/XRSC,  Wr i ght-Patter son  APB,  Ohio.  January 

1979. 

Church,  Capt  Ted,  et  al .  Personal  interviews  with  many 
personnel  from  HQ  SAC/LGXB  and  ASD/B1L.  June  through 
Mo vemb  er  1 982 . 

Shannon,  Robert.  E.  Systems  Si  myiat  i  on ._  Englewood 

Cliffs,  New  Jersey:  Prent i ce— Hal  1  ,  Inc:.,  1975. 


I 


Captain  Lance  M.  Roark  was  born  on  August.  24,  1955  in 

Burbank,  California.  He  graduated  from  high  school  in 
Miraleste,  California  in  1972  and  attended  Harvey  liudd 
College  in  Claremont,  California  from  which  he  received  the 
degree  of  Bachelor  .of  Science  in  Mathematics  in  June  1976. 
Upon  graduation,  he  was  employed  as  a  member  of  the 
technical  staff  for  Rockwell  International  Corporation,  B-l 
Division.  In  October  1978,  he  received  a  commission  in  the 
USAF  as  a  distinguished  graduate  of  the  Officer  Training 
School.  He  then  served  as  an  operations  research  analyst  in 
the  F-16  System  Program  Office  until  entering  the  School  of 
Engineering,  Air  Force  Institute  of  Technology,  in  June 
1 982 . 


Permanent  address:  6645  W.  Coley  Avenue 

Las  Vegas,  Nevada  89102 


156 


